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Abstract: The antenna subtraction method developed originally for the computation of 
higher order corrections to jet observables from a colourless initial state is extended for 
hadron collider processes involving a pair of massive particles and jets in the final state 
at the next-to-leading order (NLO) level. Due to the presence of coloured initial states, 
the subtraction terms need to be divided into three categories (final-final, initial-final and 
initial- initial). In this paper, we outline their construction and derive the necessary ingredi- 
ents: phase space factorisation, antenna functions and also integrated antennae, including 
the effects of massive final states in all of those building parts. As a first application, 
we explicitly construct the colour-ordered real radiation and the corresponding antenna 
subtraction terms required at NLO for the production of a top quark pair and for the 
production of a top quark pair in association with a hard jet. The latter constitutes an 
essential ingredient for the computation of the hadronic production of a top-antitop pair 
at NNLO. 
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1. Introduction 

At LHC, physics beyond the Standard Model will almost inevitably manifest itself by the 
creation of massive particles which decay instantaneously into multiparticle final states. 
Searches for supersymmetric particles often involve final states with four or more jets. Top 
quarks [1, 2] are measured through their decay into a bottom quark and a subsequently 
decaying IF-boson, yielding up to six-jet final states for top quark pair production. Mean- 
ingful searches for these signals require not only a very good anticipation of the expected 
signal, but also of all standard model backgrounds yielding identical final state signatures. 
Since leading-order calculations are affected by large uncertainties in their normalisation 
and their kinematical dependence, it appears almost mandatory to include NLO correc- 
tions. For a long time, these corrections were available only for at most three final state 
particles. Up to recently, the principle obstacle to NLO calculations of final states with 
higher multiplicities were the one-loop virtual corrections to multi-parton scattering ampli- 
tudes. Application of standard quantum field theory methods to those produced extremely 
large and numerically unstable expressions. In the recent past, several new algorithms 
have emerged to circumvent these numerical instabilities [3]. The recently released pack- 
ages CutTools [4], BlackHat [5], Golem [6], Rocket [7] and Samurai [8] provide automated 
implementations of these new methods. They were already applied to in a number of 
pioneering calculations [9-17]. 

With a mass m t = 173 ± 1.3 GeV, the top quark is the heaviest particle produced 
at colliders and due to its very large mass it decays before it hadronises. By studying its 
properties in detail, it is hoped to elucidate the origin of particle masses and the mechanism 
of electroweak symmetry breaking. Since its discovery at the Fermilab Tevatron, a number 
of its properties (mass, couplings) have been determined to an accuracy of ten to twenty 
per cent. With the large number of top quark pairs expected to be produced at the LHC, 
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the study of its properties will become precision physics. This large production rate will 
allow precise measurements of their properties and their production cross sections with an 
expected experimental accuracy of five per cent. 

Current theoretical predictions for the top quark pair production cross section include 
NLO corrections [16-18] and next-to-leading- logarithmic resummation (NLL) [19]. More 
recently even the NNLL resummation effects have been completed in [20] . These predictions 
lead to a theoretical uncertainty of the order of ten per cent. The same precision is available 
for single top quark production [21], top-pair-plus-jets production [15,22] and for top-pair- 
plus-bottom-pair production [13,14]. 

The top quark appears as virtual particle in hadron collider processes and due to the 
small ratio between the top quark width and its mass, it is possible to factor the cross 
section of processes involving top quarks into the product of the production cross section 
for on-shell top quarks and the top quark decay width. Most of the calculations mentioned 
above are performed for on-shell top quark pair production. Only most recently, the decay 
of the top quark has been included in NLO calculations [16, 17, 23] leading to a similar 
theoretical accuracy. 

Even for on-shell top quark pair production, a full fixed order calculation of the total 
top-antitop rate at NNLO, required to match the experimental accuracy, is missing. NNLO 
calculations involving massive quarks require the same ingredients as their massless coun- 
terparts. Three classes of contributions enter: double real, mixed real-virtual and two-loop 
type virtual contributions. However, the quark mass introduces one additional scale into 
the calculation. Especially the two-loop virtual corrections become therefore more involved 
than in the massless case. On the other hand, the treatment of the real radiation n + 2 
parton processes is expected to be easier, since the heavy quark mass acts as an extra 
infrared regulator, thus eliminating part of the singularity structure. 

Recent progresses has been accomplished concerning the two-loop contributions. Part 
of these two-loop virtual corrections are built with products of one-loop virtual amplitudes. 
Those corrections have been computed in [24]. Concerning the two- loop virtual corrections 
build with product of two-loop and tree-level amplitudes, the situation is different. The 
two-loop virtual corrections for the processes qq — > ti and gg — > ti are not fully available 
at present. A purely numerical evaluation of the quark-initiated process [25] could be 
partly confirmed by analytical results [26], which were most recently extended also to the 
gluon-induced subprocess. 

Real radiation corrections involving heavy quarks have up to now been investigated 
only at NLO. Most recently a semi-numerical scheme to evaluate the double real radiation 
corrections at NNLO has been introduced and applied to top quark pair production [27]. 

In this paper, towards calculating the double real contributions to top-antitop pro- 
duction in hadronic collisions at NNLO, we extend the antenna subtraction formalism to 
be applicable to the production of massive heavy quark pairs in the presence of coloured 
initial states and present results for the real and subtraction contributions for ti and ti 
+jet production at NLO. 

Generally, for hard scattering observables, the inclusive cross section with two incoming 
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hadrons H±, H2 can be written as 

dcJ = E/ ^^/a/i(^)A/2(6)da a6 (Ci J ffi,6^2) , (1.1) 

a,b 

where £1 and £2 are the momentum fractions of the partons of species a and b in both incom- 
ing hadrons, / being the corresponding parton distribution functions and d<7 a b(£i-ffi, £2-^2) 
is the parton-level scattering cross section for incoming partons a and b. 

The partonic cross section da a b has a perturbative expansion in the strong coupling a s 
such that theoretical predictions for a hadronic process at a given order in a s are obtained 
when all partonic channels contributing to that order to the partonic cross section are 
summed and convoluted with the appropriate parton distribution functions as in eq.(l.l). 

In general, beyond the leading order, each partonic channel contains both ultraviolet 
and infrared (soft and collinear) divergences. The ultraviolet poles are removed by renor- 
malisation in each channel. Collinear poles originating from the radiation of initial state 
partons are cancelled by mass factorisation counterterms and absorbed in the parton distri- 
bution functions. The remaining soft and collinear poles cancel among each other when all 
partonic contributions are summed over [28] . As these observables depend in a non trivial 
manner on the experimental criteria needed to define them, they can only be calculated 
numerically. The computation of hadronic observables including higher order corrections 
therefore requires a systematic procedure to cancel infrared singularities among different 
partonic channels before any numerical computation of the observable can be performed. 

For the task of next-to- leading order (NLO) calculations [29], the infrared divergen- 
cies present in real radiation contributions can be systematically extracted by process- 
independent procedures, called subtraction methods. 

More specifically, let us consider the hadronic production of m-jets at NLO. A theoret- 
ical prediction for this observable is obtained by summing the following massless partonic 
contributions: At LO, the tree- level contribution contains m partons in the final state 
which build m jets. 

At NLO, the differential cross section for the production of m-jets, d&NLO m &y sym- 
bolically be written as, 

daNLO = I da* LO 4T+ 1 ) + / {da v NLO + d&%[ ) 3^ (1.2) 

where J d $ corresponds to the integration over the N parton phase space. The cross 
section is built with the real radiation cross section contribution da R which has (m + 1) 
massless partons in the final state, the one-loop cross section da v and the mass factorisation 
counterterm da^[ Q which have both m partons in the final state. A jet algorithm is applied 
separately on each of these contributions to ensure that out of n partons, m jets are built 

(n) 

in the final state. Symbolically this recombination procedure is denoted by Jm ■ 

The purpose of any subtraction method is to provide a subtraction term da^ LO which 
has the same singular behaviour as the real radiation squared matrix element and is suf- 
ficiently simple to be integrated analytically over a factorised form of the (m + l)-phase 
space. 
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Using a subtraction method, the NLO partonic cross section given in eq. (1.2) may 
then be written as, 



With this, the first integral is finite and can be integrated numerically in four dimen- 
sions. The integrated form of the subtraction term has the same number of final-state 
partons as the virtual contributions and the mass factorisation counterterms. It can there- 
fore be combined with those thereby canceling analytically the explicit infrared divergences. 
The second integral in eq. (1.3) is therefore finite as well. 

The actual form of the subtraction term dafj LO depends on the subtraction formalism 
used. The approximations of the matrix-element in the unresolved limits being non unique, 
several successful subtraction formulations have been proposed in the literature [30-34]. 
The dipole formalism of Catani and Seymour [30] and the FKS [31] of Frixione, Kunszt and 
Signer have been implemented in an automated way in [35-40]. In its original formulation, 
the formalism of Catani-Seymour [30] deals with massless partons in final and/or initial 
state at NLO. 

An alternative subtraction formalism is given by the the antenna subtraction method 
. This formalism [41,42] was originally derived for processes involving only (massless) 
final state partons in e + e~ collisions. It has been applied in the computation of NNLO 
corrections to three-jet production in electron-positron annihilation [43-46] and related 
event shapes [47-51], which were subsequently used in precision determinations of the 
strong coupling constant [52-56]. 

For processes with initial-state partons and massless final states, the antenna sub- 
traction formalism has been so far fully worked out only to NLO in [41,57]. It has been 
extended to NNLO for processes involving one initial state parton relevant for electron- 
proton scattering in [58] while an extension of the formalism to include two initial state 
hadrons at NNLO is under construction [59,60]. 

Subtraction formalisms which deal with massive final state particles have been so far 
only developed up to the NLO level [61-64]. The kinematics is more involved due to the 
finite value of the parton masses. QCD radiation from massive particles can lead to soft 
divergencies but cannot lead to strict collinear divergencies, since the mass is acting as an 
infrared regulator. In a calculation of observables involving massive final state fermions, 
logarithmic terms of the form ln(Q 2 /M 2 ), where M is the parton mass and Q, the typical 
scale of the hard scattering process can be generated. In kinematical configurations where 
Q » M, these logarithmically enhanced contributions can become large and can spoil the 
numerical convergence of the calculation. The cross section calculation of ti at LHC is an 
example where such enhanced logarithmic terms arise. These terms are related to a process- 
independent behaviour of the matrix elements; its singular behaviour in the massless limit 
(M —7- 0). This singular behaviour is related to the quasi- collinear [61] limit of the matrix 
element. In the presence of massive particles, the factorisation properties of matrix element 




(1.3) 
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and phase space in collinear and soft limits need to be generalised to take the mass effects 
into account. 

The dipole formalism of Catani and Seymour [30] and the FKS formalism [31] have 
been extended to deal with massive particles in [61] for colourless and for coloured initial 
states up to the NLO. The subtraction terms constructed within this formalism account 
for the quasi-collinear limits of the matrix-element squared. 

The antenna subtraction method has so far only been extended to deal with the produc- 
tion of massive fermions from a colourless initial state in [64] . It is the purpose of this paper 
to present an extension of this method to include radiation off final state massive fermions 
produced in hadronic processes. In this paper we aim to derive all necessary ingredients, 
massive antennae, phase space factorisation and finally integrated massive antennae for 
this extension. As a first application, we construct the subtraction terms daff LO required 
for the hadronic production of a top quark pair in association with no or one hadronic jet. 

These two processes, pp — > ti and pp — > ti +jet have been calculated in [19, 22] using 
the dipole subtraction method [61]. Our aim here is however not to redo this calculation 
using another subtraction scheme. Instead, we construct here the subtraction terms for 
ti +jet production in a colour-ordered form which is essential for the computation of top 
quark pair production without any jets at NNLO within the framework of the antenna 
formalism. Those subtraction terms can be used to capture all single unresolved radiation 
from the double real radiation contribution for the ti pair production at NNLO. 

The plan of this paper is as follows. In Section 2 we outline the construction of the 
subtraction terms for the hadronic production of a heavy quark pair in association with 
(m — 2)-jets at NLO. We present the form of the subtraction terms in all kinematical con- 
figurations with particular emphasis on the changes caused by the presence of massive final 
state particles in the expressions of the subtraction terms compared to those expressions 
in the massless case. Section 3 contains a list of all massive antenna functions required. In 
Section 4 we tabulate all non- vanishing single unresolved limits of those massive antennae 
while in Section 5 results for the integrated antennae are given. Section 6 presents a check 
on one of the integrated antenna. In Section 7, for all partonic process involved, we present 
the colour ordered real contributions and their subtraction terms required to evaluate the 
hadronic production cross section of a tt pair and of tt pair and a jet at NLO. Finally 
Section 8 contains our conclusions. 

2. Antenna subtraction with massive final states 

In this section, we present the general formalism necessary to evaluate the hadronic pro- 
duction of a pair of heavy quarks QQ in association with (m — 2) jets at the next-to- leading 
order (NLO) in perturbative QCD. 

2.1 Real radiation contributions to heavy quark pair production 
in association with jets 

The leading order (LO) m-parton contribution to the hadronic production of a pair of 
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heavy quarks QQ in association with (m — 2) jets may be written as, 
da LO {pi,p 2 ) = ^2 d$ m (kQ,kQ, ki, . . . , k m - 2 ;pi,p 2 ) 

m-2 

m—2\Pl i P2 

)| 2 J^\k Q ,k Q MT--,k m - 2 ). (2.1) 

The momenta pi and p 2 are the momenta of the initial state partons, the massive partons 
Q and Q have momenta feg and fcg, while the momenta of the remaining (m — 2) massless 
final state partons are labelled k± . . . k m - 2 . S m - 2 is a symmetry factor for identical massless 
partons in the final state. Jm l \kQ, kg, k±, . . . k m -2) is the jet function. It ensures that out 
of (m — 2) massless partons and a pair of heavy quarks Q and Q present in the final 
state at parton level, an observable with a pair of heavy quark jets in association with 
(m — 2) jets is built. At this order each massless or massive parton forms a jet on its 
own. The normalization factor N includes all QCD-independent factors as well as the 
dependence on the renormalised QCD coupling constant a s . Ylm-2 denotes the sum over 
all configurations with (m — 2) massless partons. d^m, is the phase space for an m-parton 
final state containing (m — 2) massless and two massive partons with total four-momentum 
Pi +P2- In <i = 4 — 2e space-time dimensions, this phase space takes the form: 

d^kn d^kn 
d* m (k Q , k Q M,-.., k m - 2 ; Pl , P2 ) = 2gq(2?r)d _ 1 2E^(2^ 

In eq.(2.1) lA^m) 2 denotes a colour-ordered tree-level m-parton matrix element squared for 
m partons out of which two are massive. These terms only account for the leading colour 
contributions to the squared matrix elements. On the other hand, colour subleading con- 
tributions are, in general, given by the interference between two colour-ordered n-parton 
amplitudes. However, to keep the notation simpler we denote these interference contri- 
butions also as |.M m | 2 . Related to these interference terms, it is here worth noting the 
following: As soon as more than five coloured partons are present in a given partonic pro- 
cess, the subtraction of infrared singularities present in interference terms is more involved 
than for colour ordered squared matrix elements. This particular issue will be treated in 
Section 7.2.1. 

At, NLO the real radiation partonic contribution to the production of the heavy quark 
pair production in association with (m — 2) jets involves (m + l)-final state partons with 
two of them being massive. It may be written as, 

d-°NLo(PiiP2) = ^2 d<f> m+1 (k Q , kg, ki, . . . , k rn ^ l ;p 1 ,p 2 ) 

m+l 

x-^— \M m +i(k Q , k Q , h, ....... , km-i,;Pi,P2)\ 2 Jin +1) ( k Q, kn,ki,..., km-!). (2.3) 

(tyi I 1^ — 

The jet function J m ensures that out of (m — l)-massless partons and a QQ pair, an 
observable with a pair of heavy quark jets in addition to (m — 2) jets, are built. In other 
words, an m-jet observable is formed. 
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In this contribution, when the real matrix element squared |A4 m +i| 2 is integrated over 
the phase space, it develops singularities when one parton in the final state is unresolved. 
In the presence of massive partons in the final state, a parton is called unresolved, either 
when it becomes soft or collinear to another massless parton or when it is quasi- collinear 
to a massive parton. In this latter case, it leads to finite logarithmic terms involving the 
mass of the massive parton. The notion of quasi- collinear limit will be explicitly presented 
in Section 4. To extract the unresolved behaviours of the real matrix element, subtraction 
terms which take both the massless and massive effects need to be considered. 

At the next-to-leading order, the subtraction terms derived in the antenna formalism 
[41, 42, 57] are constructed solely with tree-level three-parton antenna functions. Those 
functions encapsulate all singular limits due to the emission of one unresolved parton 
between two colour-connected hard partons, called radiators. 

Depending where the two radiators are located, in the initial or in the final state, we 
distinguish three types of configurations: final-final, initial-final and initial-initial. In any 
of those configurations, the radiated parton is always in the final state. 

The subtraction terms in a given configuration are constructed from products of the 
corresponding antenna functions with reduced matrix elements. Those can be integrated 
over a phase space which is factorised into an antenna phase space (involving all unresolved 
partons and the two radiators) multiplied by a reduced phase space, where the momenta of 
radiators and unresolved radiation are replaced by two redefined momenta. These redefined 
momenta can be in the initial or in the final state depending on where the corresponding 
radiator momenta are and are defined by appropriate mappings. The full subtraction 
term is then obtained by summing over all antennae required in one configuration and by 
summing over all configurations needed for the problem under consideration. 

The antenna subtraction terms do not provide a strictly local subtraction of collinear 
singularities in the case of a gluon splitting to two gluons or to a quark-antiquark pair. In 
these, the antenna subtraction term accounts for the singular behavior only after integration 
over the azimuthal angle of the two parton system with respect to the collinear direction. 
As a consequence, the numerical integration of the difference of matrix element and antenna 
subtraction term is potentially unstable. By an appropriate partitioning of the final state 
phase space [43,65], this azimuthal integration variable can be separated off for each limit. 
Once this variable is separated, the angular terms can be averaged out by a smooth one- 
dimensional integration, or by combining different phase space points. 

The massless and massive three parton final-final antenna functions, besides being 
fundamental entities of the antenna subtraction formalism developed for colourless ini- 
tial states and for massless partons in [42] or for massive partons in [64] have another 
fundamental role. Those can be used as basic ingredients in parton showers. The event 
generator VINCIA uses these antenna functions as evolution kernels. In its present formula- 
tion [66,67], VINCIA describes the evolution of timelike showers based on the massless [42] 
and the massive [64] antenna functions. A study concerning the importance of quark mass 
effects is currently ongoing. The results presented in this paper concerning the massive 
initial-final antennae will become relevant for initial-state parton showers for observables 
with massive final states. 
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For the NLO corrections to QQ + jets production in hadronic collisions, we will need 
all three types of subtraction terms and therefore all three types of antenna functions. 
Since one massive radiator is always in the final state, the subtraction terms will involve 
final-final and initial-final massive antennae but no initial-initial antennae with massive 
particles. Antennae involving only massless partons though, will be needed in all three 
configurations. Those have been derived and integrated in [42,57,60]. All required massive 
antenna functions will be presented in Section 3 and their integrated forms will be given 
in Section 5. 

In the following, we shall present the general form of the subtraction terms needed 
in each of the three configurations (final-final, initial-final and initial-initial) to account 
for single unresolved radiation in processes involving a heavy quark pair in association 
with jets in the final state. We will in particular focus on the changes introduced in the 
subtraction terms due to the presence of massive final state particles compared to those 
when only massless partons are involved [42,57,60]. 

2.2 Subtraction terms for final-final configurations 

In the final-final configuration, the subtraction term related to the real contributions to 
the partonic process yielding a heavy quark pair in association with (m — 2) jets given in 
eq.(2.3) has to take into account the presence of an unresolved parton j of momentum kj 
emitted between two hard final-state radiators i and k of momenta ki and k^ respectively. 
It reads, 



This subtraction term involves the phase space for the production of m + 1 partons, d<3? m +i, 
the massive final-final three-parton antenna function X^- k , the reduced m-parton amplitude 

squared |.M m | 2 and the jet function The jet function ensures that out of (m — 2) 

massless partons and a pair of massive partons, (m — 2) jets and a QQ jet pair is build. 
The jet function and the reduced m-parton amplitude do not depend on the individual 
momenta ki, kj and kk, but will only depend on the redefined momenta kj and kx which 
are linear combinations of the original momenta ki, kj, k^. 

Two cases are implicitly considered here. Either i and k are massive hard final state 
radiators in which case the redefined partons I and K are massive or, i is massless and k is 
massive and the redefined partons / and K are massless and massive respectively. In this 
case, one of the parton momenta k a with a ^ i, k is massive in order to obtain a reduced 
matrix element with two massive final state partons. 

Eq.(2.4) holds strictly fot the subtraction of singularities of colour ordered matrix 
elements squared. Furthermore, as mentioned before, also interferences between partial 
amplitudes with different colour orderings appear in the subleading colour pieces. It will 




x 



Yl X ijk IMmih, ...,K!,K K ,.. . ,k m+1 ; Pl ,p 2 )\ 2 



j 



xJ^\k 1 ,...,K I ,K K ,...,k m+l ). 



(2.4) 
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be seen in section 7.2.1 that the subtraction of infrared singularities appearing in these 
interferences needs a special treatment. However, to keep equations as brief and clear as 
possible, we still write our subtraction terms in the final-final configuration symbolically 
as in eq.(2.4). 

Most of the massive final-final three-parton antenna functions X^- k needed to evaluate 
the NLO corrections to QQ + jets production have been derived in [64]. Those will be 
listed in Section 3. Solely, the new flavour-violating massive antennae, which are related 
to flavour violating vertices, will be explicitly derived in that section. 

The phase space d<I> m+ i can be factorised as follows, 

d$ m+ i(fci, . . .,ki,kj,k k , . . . , k m+1 ;p 1 ,p 2 ) = 
d$ m (fci,..., if/, K K ,...,k m+1 ;p 1 ,p 2 ) ■ d<f> Xijk (ki, k j: k k \Ki + K K ) . (2.5) 

d<E> m is the ti-dimensional phase space for m outgoing particles with momenta k±, ■ ■ ■ , k m +i 
with two of those momenta being massive, d&x ijk is the NLO final-final antenna phase 
space. It is proportional to a massive three-particle phase space relevant to a 1 — > 3 
decay [64]. Depending whether the two radiators are of equal masses or whether one of 
them is massless, different parametrisations of this antenna phase space are obtained. The 
parametrisations necessary to integrate the final-final massive antenna functions will be 
given in Section 5. 

Appropriate final-final phase space mappings are furthermore required to define the 
final-state momenta fcj and kx in the relevant kinematical configurations. Those are how- 
ever not unique. A possible mapping can be found in [61]. As in the massless case, the 
mapping given in [61] is not symmetric under the exchange of i and k. A possible sym- 
metric version for a massive final-final definition of the mapped momenta will be given 
elsewhere. 

For the analytic integration, we can use the phase space factorisation formula given in 
eq.(2.5) to rewrite each of the subtraction terms in eq.(2.4) in the form 

\M m (k 1 ,...,K I ,K K ,...,k m+l ; Pl ,p 2 )\ 2 jt ) d<l>m J d$ Xi . k X G ijk . (2.6) 

The integrated massive final-final antennae, normalised appropriately are defined by 
analogy to the massless case by 

= cji) / d **« A ( 2 - 7 ) 

with, the normalisation factor 

C(e) = (^T^-r- (2-8) 

The integration is performed analytically in d-dimensions such that the integrated subtrac- 
tion terms can be combined with the one-loop m-parton contribution. 
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2.3 Subtraction terms for initial-final configurations 

In the initial-final configuration, the subtraction term related to the real contributions given 
in eq.(2.3) has to take into account the presence of an unresolved parton j of momentum kj, 
(which can be massive), emitted between a massless initial state radiator i of momentum 
Pi and a massive final state radiator k of momentum kk and mass m k . It is given by, 



d& NLO = N ^2 d ^m+i{ki,---,k j ,k k ,...,k m+1 ;p i ,p 2 )-^— 
m+l * m+1 

Yl X lik \M m (k u . . . , K K , . . . , k m ; Pi ,p 2 )\ 2 J^\h, ...,K K ,..., k m ). 

(2.9) 

This subtraction term involves the phase space for the production of (to + 1) par- 
tons, d$ m +i, the massive three parton initial-final antenna function denoted by X^- k , the 

reduced m-parton amplitude squared |.M m | 2 and the jet function jffl. The reduced m- 
parton matrix element squared |.M m | 2 does not contain any explicit dependence on the 
original final state momenta kj and k k , but only depends on them through the redefined 
momentum Kk- The same holds for the jet function jffl. The three parton massive 
antenna functions Xfj k depend only on the momenta pi, kj and k k and on the masses of 
the final state partons rrij and m k . These will be presented in Section 3. 

In this configuration, the initial state radiator i is always massless while the final state 
radiator k is always massive. The unresolved parton j can be either massless or massive 
such that two situations are covered in eq.(2.9). In case it is massless, one of the partons 
with k a where a / j,k has to be taken massive in order to build an event with a pair of 
massive partons in addition to (to — 2) jets. Depending whether the unresolved parton 
is massless or massive, different phase space parametrisations will be needed. Those will 
be presented in Section 5 when the massive initial-final antennae are integrated over the 
unresolved phase space. 

As it was mentioned for the final-final case, the subtraction terms as presented in 
eq.(2.9) are strictly valid for the subtraction of infrared singularities of colour ordered 
matrix elements squared. Interferences of different partial amplitudes which also generate 
infrared singularities need a special treatment which we shall present in section 7.2.1. To 
keep our equations as brief and clear as possible, we still write our subtraction terms in 
the initial-final configuration symbolically as in eq.(2.9). 

2.3.1 Phase space factorisation and mappings 

For the subtraction term to be integrated and then added to the virtual contributions and 
mass factorisation counterterms, we need the phase space to factorise adequately. As the 
presence of massive final states lead to differences in this factorisation with respect to the 
massless case considered in [57] we derive this phase space factorisation explicitly below. 

We start from the 2 — > (to + l)-particle phase space with a priori two massive final 
state particles of momenta kj and k k and the remaining (m, — 1) particles staying massless. 
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The momenta of the initial state partons are named pi and p 2 - 

d$ m+ i (k u ...,k m+1 ;pi,p 2 ) = (2vr) d S^ Pt +p 2 -Y j k^j 11 [dk] , (2.10) 

where the phase space measure for massless and massive partons are respectively given 
by: [dk{\ = d d k S + (kf) / (27r) d_1 for I ^ j,k and [dfc,] = d% 5+ (kf - mf) / (27r) d_1 for 
/ = j,k. We insert 

1 = J d d q8(q + Pi -kj-k k ), (2.11) 
and, to take into account the masses of the final state particles j and k we insert 

1 = Q2 \f K I ^ J[dK K ](2ir) d 8( q + Pl -K K ) (2.12) 

with 

Q Q (c^ ~\~ Tfi jy- 

Q = -q x i = ^, and Pl=XiPi- (2.13) 

2pi • q 

We find that the original phase space for (to + 1) partons given in eq.(2.10) can be written 
as the product of an m-parton (d<I> m ) and a two-to-two particle phase space (d<&2) as 
follows, 

d$ m+ i(fci, k rn+1 ;pi,p 2 ) = d$ m (fei, K K , k m+1 ; Xip h p 2 ) 

x (Q 2 + O d$2(fe h . Phq) ^i. (2 . 14) 

Replacing this phase space in the subtraction term given in eq.(2.9) we can explicitly carry 
out the integration of the antenna functions over the two-to-two particle phase space and get 
the integrated expression for the subtraction term to be added to the virtual contributions 
and the mass factorisation counterterms. For this purpose it is convenient to define the 
integrated massive initial-final antenna as, 

x 1 frit (Q 2 + Q Y r91 ^ 

" jk C(e) J 2vr ' jk ' ' ' 

where C(e) is given in eq.(2.8) and the initial-final massive antenna phase space denoted 
by d$ Xi , jk is given by 

d^ ]k =d ^ Q2+ 2 f K \ (2.16) 

Due to the presence of massive particles, the phase space mapping for initial-final 
configurations derived in the massless case in [57] needs to be generalized as well: In a 
process of the form pi + q — > kj + k^ with (q 2 < 0) , to find a mapping from the original 
momenta pi in the initial state, kj and kk in the final state, to the redefined momenta 
pi and Kk that ensures phase space factorisation, the following conditions need to be 
fulfilled: The remapped final state momentum Kk must be on-shell, and momentum must 
be conserved. This implies, 

Pi - kj - k k = XiPi - K K =pi - K K (2.17) 
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for the phase space to factorise as in eq.(2.14). 

In addition, concerning the single unresolved behaviour of the remapped momenta, the 
following requirements have to be satisfied: 

XiPi —s- Pi Kk — > kk when j becomes soft 

XiPi —s- Pi Kk — > kj + kk when j becomes collinear or quasi-collinear with k 

XiPi Pi — kj Kk — > kk when j becomes collinear or quasi-collinear with i. 

(2.18) 

It can be seen that when parton j becomes soft, collinear or quasi-collinear to parton k we 
have Xi — >■ 1, while for the case in which partons i and j are collinear or quasi-collinear, 
Xi — > 1 — z, with z being the fraction of the initial state momentum pi carried by the 
unresolved momentum j. The notion of quasi-collinear limit can be viewed as an extension 
of the collinear limit when at least one of the partons becoming collinear to another is 
massive. Its explicit definition in this initial-final configuration will be given in Section 4. 
All these conditions given above are interpolated by, 

Kk = kj + k k - (1 - xi)pi 

Xl = *i (2.19) 

&ij Sik 

With this choice, the mass relation for momenta before and after mapping reads, 

m\ = m 2 k + m). (2.20) 

Note also that in eq.(2.19), s a b stand for 2p a • p^. as everywhere in this paper. In terms of 
those, massless [57] and massive definitions of Xi are the same. The massless phase space 
factorisation and mappings given in [57] can be recovered by setting the masses of the final 
state particles to zero. 

2.4 Subtraction terms for initial-initial configurations 

Additional divergent contributions may finally also occur in the real matrix element squared 
given in eq.(2.3) when a massless final state parton becomes unresolved with respect to 
two initial state partons. In this case, the subtraction terms are constructed exclusively 
with massless three parton initial-initial antennae. In those, the initial state partons are 
the hard radiators and this situation was studied in detail in [57]. 

The subtraction term associated to an unresolved massless parton j and two hard 
initial state radiators i and k with momenta pi and pk in the partonic process Pi + Pk 
kj + kQ + kq + (m — 2) partons takes the form 

da NLO = N d$ ™+l( fc Q' fc Q' fc l' ■ , km-^PuPk)-^— 

m+l ' m+1 

x ^2 X ik,j \ M m(kQ, kg, h, . . . , kj-^kj+i k m -i;xiPi,x k Pk)\ 2 
j 
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x 4i \^Qi kQ,kl, ... , kj—l, • • • , • • • , k m —l). 

(2.21) 

This subtraction term involves the phase space for the production of (m + 1) partons, 
d3> m +i, the massless three parton initial- initial antenna function denoted by Xf k j, the 

reduced m-parton amplitude squared |A4 m | 2 and the jet function jffl. As mentioned for 
the final-final and initial-final cases, the subtraction terms as presented in eq.(2.21) are 
strictly valid for the subtraction of infrared singularities of colour ordered matrix elements 
squared only. Interferences of different partial amplitudes need a special treatment which 
we shall present in Section 7.2.1. However, to keep our equations as brief and clear as 
possible, we still write our subtraction terms in the initial-initial configuration symbolically 
as in eq.(2.21). 

All momenta (massless and massive) in the arguments of the reduced matrix elements 
|.M m | 2 and the jet function jffl have to be redefined. They are denoted with tildes in eq. 
(2.21). The two hard radiators are simply rescaled by factors x% and Xk respectively. The 
other momenta present in the reduced matrix element squared are boosted by a Lorentz 
transformation onto the new set of momenta {ki, I ^ j} as described for massless partons 
in [57]. We have checked that the same boost [30] required to redefine the momenta 
of the massless partons can be used to redefine the momenta of the massive partons. 
Consequently, the presence of massive partons in the final state does not influence the way 
the phase space factorises and how the mapping is defined so that the massless factorisation 
[57] can be used in this context. 

The phase space d$ m +i factorises into the convolution of a massive m-particle phase 
space involving only redefined momenta and a massless initial-initial antenna phase space 
related to the phase space of the unresolved parton j, as in the massless case. 

All required initial-initial massless antennae X® k ■ needed for the construction of the 
subtraction terms for ti production in association with jets have been derived and integrated 
in [57,60]. 

3. Massive antenna functions 

In this section, we aim to list all massive antenna functions which are needed to con- 
struct the subtraction terms for the hadronic production of a heavy quark pair and for the 
production of a heavy quark pair in association with one additional jet at NLO. 

3.1 General features of antenna functions 

In Section 2, we saw that the subtraction terms defined in the three configurations (final- 
final, initial-final and initial-initial) are built with products of the corresponding antenna 
functions (denoted by X) with reduced matrix element squared. At this order (NLO), 
only tree level three-particle antenna functions are required. Those describe all unresolved 
(soft, collinear and quasi-collinear) radiation emitted between a pair of colour-ordered hard 
partons, the radiators. Originally, in [42], the antenna functions were derived for massless 
final-final configurations. Those are defined by the pair of hard partons they collapse to 
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in the unresolved limits and in all cases, the antenna functions are derived from physical 
matrix elements. Generally, the quark-antiquark antenna functions are obtained from 
7 — ^ QQ + (partons), the quark-gluon antenna functions from \ ~ > 9 + (partons) [68] and 
the gluon-gluon antenna functions from H — > (partons) [69]. 

The three parton final-final antenna functions were obtained by normalising the colour- 
ordered three-parton tree level matrix elements squared to the matrix element squared for 
the basic two-parton process, omitting all couplings and colour factors. As such the tree- 
level three parton massless or massive final-final antenna functions are scalars in colour 
space, have mass dimension (—2) and are defined by [42,64] as, 



S denotes the symmetry factor associated to the antenna, which accounts both for potential 
identical particle symmetries and for the presence of more than one antenna in the basic 
two-parton process. It is chosen such that the antenna function reproduces the unresolved 
limits of a matrix element with identified particles. 

At NLO the existing three parton tree level massless final-final X^- k antennae are: 

• Quark-antiquark: the only antenna functions of this kind at NLO are the A- Type 
antennae, and they are obtained from the ratio \M(^* — > qqg)\ 2 /|A4(7* — > qq)\ 2 . 
Since the quark and the antiquark are of the same flavour, in the following these 
antennae will be referred as flavour conserving A- Type antennae. 

• Quark-gluon: there are two different antenna functions of this kind: D-Type and 
E-Type. The D-Type antennae are obtained from the ratio |A4(x — > ggg)\ 2 /\M(x — > 
gg)\ 2 , while the E-Type are computed from \M(x — > gqq)\ 2 /\M(x - > <7fi0| 2 [68]. 

• Gluon-gluon: there are also two different antennae of this kind: F-Type and G-Type. 
The F-Type antenna functions are obtained from \M(H — > ggg)\ 2 /\M(H — > gg)\ 2 , 
while the G-Type are calculated from the ratio \M(H ->■ gqq)\ 2 /\M(H -> gg)\ 2 [69]. 

The initial-final antennae and the initial-initial antennae denoted respectively by Xfj k 
and Xf k j are in principle defined by crossing one or two massless partons from the final 
to the initial state in the final-final antennae X^- k . However, this crossing may not be 
unambiguous for initial-final configurations [57]. 

All these types of antenna functions are needed for the construction of the subtraction 
terms for heavy quark pair production except the G-Type antennae. Those antennae are 
required in all three (final-final, initial-final or initial-initial) configurations and in their 
massless or massive versions. In the latter case, the final state quark (q) or gluino (g) 
is taken massive. Note that if both of those partons are present in the final state, the 
corresponding massive antenna is obtained by taking only the gluino g to be massive. The 
antenna obtained taking both of these partons massive is not required for the construction 
of the subtraction terms for heavy quark pair production +jets at NLO, considered here. 



X 



-o 

ijk 




(3.1) 
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In addition to these flavour conserving antennae, we will also need flavour violating 
massive quark-antiquark A- Type antennae in both final-final and initial-final configura- 
tions. Those involve flavour violating vertices and will be defined below. 

In all these antennae denoted by in the final-final configuration and crossings of 
those, the labels i,j, k can either stand for massless or massive partons. The particle with 
label j always denotes the unresolved parton. The massive antennae will depend on the 
invariants Sij , Sik and Sjk (with s a b = 2paPb throughout this paper) and on the masses of 
the final states heavy partons. 

In the following, when presenting the expressions of the required massive antennae we 
will use a clear labeling of those antennae. We shall specify which partons in them are 
taken massless, those will be indexed with q, or taken massive and then indexed with Q. 
For conciseness, the 0(e) in the expression of the unintegrated antennae will be omitted. 

3.1.1 Flavour- violating antennae 

For the construction of the subtraction terms for ti +jet production at NLO, we encounter 
two types of flavour violating antennae. Both are A- Type quark-antiquark antennae and 
are labelled : A Q q an d A° qg Q. I n those, q represents a massless quark and Q and Q 
represent a massive quark or antiquark of different flavour than the massless quark q. Note 
that for symmetry reasons of the A- Type antenna, the role of q and Q in A^q or of q 
and Q in A® g Q are interchangeable resulting in the same two flavour-violating antennae. 
We encounter these two types of flavour violating antennae in final-final and initial-final 
configurations. In the latter case, a massless quark is always in the initial state and plays 
the role of the initial state radiator. It is worth noting that the limits covered by flavour 
violating antennae with a gluon in the initial state can be covered by flavour conserving 
initial-final antennae, such that these gluon initiated flavour violating antennae are not 
required here. 

The massive final-final flavour violating antennae have their massless counterparts 
given by A qgq > and A qgq i where q' and q' are massless quarks or antiquarks of different 
flavour than q. These massless counterparts have been used as essential ingredients for 
the construction of the subtraction terms required for the computation of e + e~ — > 3 jets 
at NNLO [43]. Those massless A- Type flavour violating antennae have exactly the same 
properties and singular structure as a massless flavour-conserving A- Type quark-antiquark 
antenna given by A qgq [42]. Those have also the same unintegrated form as this antenna. 

In the massive case, the final-final flavour violating A- Type antennae A^q and A qg Q, 
having one massive particle in the final state, have the same properties and singular struc- 
ture as a massive flavour-conserving A- Type antenna of the form Aq^q where the massive 
quark Q or the massive antiquark Q (by symmetry arguments) present in this flavour 
conserving antenna is taken massless. 

By construction, within the antenna formalism, in the unresolved (collinear, soft or 
quasi-collinear) limits, the three parton antennae yield massless or massive universal un- 
resolved factors given either by splitting functions or eikonal factors in their massless or 
massive forms. The spin properties of the final states are determinant to obtain these 
unresolved factors. Consequently, all flavour-violating antennae (massless or massive) can 
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be generated by processes resulting from the decay of a charged W-boson into massless or 
massive final state fermions. When the particle in the final state is a (massless or massive) 
antiquark such that the flavour violating antenna contains a quark and an antiquark of 
different flavours, those final states can be Dirac fermions. When the particles in the an- 
tennae are two quarks of different flavours, those final states can be produced in extended 
versions of the Standard Model [70] and will involve Dirac and Majorana fermions. 

3.2 Massive NLO final-final antenna functions 

For the construction of the subtraction terms relevant for heavy quark pair production 
at hadron colliders at NLO, we will need three-parton massive final-final flavour conserv- 
ing and flavour violating A- Type antennae and flavour conserving quark-gluon antenna 
functions. All flavour-conserving massive final-final antennae required have been derived 
in [64]. They are given below for completeness. The flavour violating massive final-final 
antennae are new. 

The generic process considered to define the massive three-parton final-final antenna 
X®j k is given by q — > k{ + kj + kk where q is the virtuality of the colourless initial state, 
with q 2 > 0. In the following, the center of mass energy of the decaying particle E cm will 
be used instead of q with q 2 = E 2 m The final state radiator partons i and k are either both 
massive or only one of them is such. Parton j defines always the unresolved massless final 
state parton. 

The massive three-parton final-final antenna will depend on on the invariants sij , 
Sjk and Sik, on the masses ra; and of the final state partons i and k and on E 2 m . 

3.2.1 Quark- Antiquark antennae 

The flavour conserving quark-antiquark massive antenna function denoted by ^Q g Q nas a 
massive quark Q and a massive antiquark Q as radiators and reads, 

Al(l Q ,3 9 ,2 Q ) = /^ + ^ + 2^ + £ 23 +£l 3 

4^E 2 m + 2m 2 Q j \S13S23 Sis S 23 S 13 S23 
,2 



( 8si2 _ 2s 12 _ 2s 12 _ 2s 23 2 2 2s 13 ^ ( 3 2 ) 



+ m O I 9 9 

\S13S23 Sf 3 Sfo Sf 3 Sl 3 S 2 3 % 



\ ^23 b 13/ / 

This function is normalised to the two-particle process 7* — > QQ, whose matrix element 
squared (omitting couplings) is given by 

A° 2 (1 Q ,2 Q ) = 4 [(1 - e) E 2 m + 2m 2 Q ]. (3.3) 

3.2.2 Quark-Gluon antennae 

For the subtraction terms required for the production of heavy particles in addition to jets, 
only massive quark-gluon antennae with one massive radiator in the final state are needed. 

The quark-gluon massive antennae with either two gluons or a massless quark-antiquark 
pair in the final state are normalised by the two-particle matrix element squared relevant 
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for the process x ~~ ^ 99> with the gluino g being massive with mass rriQ and the gluon g 
being massless. This two-particle matrix element squared omitting couplings reads, 

X 2 °(l Q ,2 g ) = 4(l-e) {E 2 cm -ml)\ (3.4) 

where X can stand for the E or D-Type antennae here. 

The quark-gluon massive antenna Eq^-, with a pair of massless quarks q'q' and a 
massive radiator quark Q in the final state reads, 

E° 3 (1 Q , 3 q ,, 4 ? = 1 4 (s 13 + Sl4 + £13 + £m _ s^mg) + O (e) . (3.5) 

The quark-gluon massive antenna with two gluons and a massive radiator Q in 

the final state reads, 

4s?, 4s?, 4s?, 2sf 3 



^3°(lQ,3 9 ,4 9 ) = — - 1 4 ( ( 9 s 13 + 9s 14 + Iffa + l£M + Iffs + 



4(£'2 m - m 2 Q ) 2 V V «14 si3 s 34 S14S34 

6S13S14 4sf 4 2sf 4 3si 3 s 34 , 3s i4 s 34 s§ 4 s§ 4 
H i 1 h os 3 4 H 1 1 1 

«34 S34 S i3 S 34 S M S13 S13 S14 

m 2 f fi , 2s ?3 , 4s i3 , 4si4 2sf 4 6s 34 4si 3 s 34 

— m O D H o 1 1 1 5 1 r 



Q I u "I 2 1 1 1 2 

V s 14 s 14 s 13 Sf 3 S13 S| 4 

6s 34 4si 4 s 34 2s| 4 , 2s| 4 , 2s 2 



^ 1 34 , ^ J 34 , ^ 34 

"I 1 9 1 9 1 9 1 

S14 Sf 3 S| 3 Sf 4 Sl 3 Sl 4/ 

„ p ^ 3 2s 34 4 2s 3 4 \ 

+2£ cm m Q - E cm m Q h m Q + O (e) . 

V Sl 3 Si4 V S13S14/ 

(3.6) 

This tree-level antenna function D 3 (1q, 3 S , 4 S ) contains two antennae, corresponding 
to the following configurations: gluon (3 g ) radiated between the massive quark and gluon 
(4 g ) and also the configuration where gluon (4 S ) is radiated between the quark and gluon 
(3 S ). The separation between these two configurations is not free from ambiguity, since 
the collinear limit of the two gluons has to be split between the two configurations. We 
consider the following decomposition 

D 3 (l Q ,3 g ,4 g ) = 4(l Q ,3 g ,4 g ) + 4(l Q ,4 9 ,3 g ) (3.7) 
where the sub-antenna denoted by is given by 

9s i3 + 9si4 + 6s 34 4sf 4 4sf 4 2sf 4 3si 3 si 4 3si 4 s 34 s^ 
2 si 3 s 34 Si 3 S 34 S 34 S13 Si 3 

-m 2 (s + ^ + lfM + ^ + ^ + lflpi + ^fi + ^M (3.8) 

V s 13 s 13 «13 S13 Sf 3 Sf 3 Si 3 Si 4 / 

+E cm m Q - (E cm - m Q )m% 534 + 0(e) 

V Si 3 Si4_ 

Both of these sub-antennae will be needed individually to construct the subtraction 
term for tt +jet at NLO. Those will not need to be integrated separately though. 
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3.2.3 Flavour violating antennae 

In addition to the massive flavour-conserving final-final antennae given above, for the con- 
struction of the subtraction terms for ti +jet production at NLO, we also need two types 
of massive flavour violating final-final antennae. We will require an antenna involving as 
radiators: one massive quark Q and one massless antiquark q , denoted by A Q Q g ^ and one 
antenna with a massive quark Q and a massless quark q denoted by Aq 9 ~ . The expressions 
of these two antennae are related by the exchange of a massless quark versus an antiquark 
in the final state. Therefore, we here present only one of those namely A°Q g -. It is given 
by 



4j(l ,3 ff ,2,-) 



2S12 2si2 2s\ 2 £13 £23 
(E 2 m + m? Q ) [ Sl3 S 23 S13S23 S23 S13 



2 / 2si2 2s 23 . 2 \ 
-m n I —5- H 5- H 

s 13 s 13 s 13/ 



+ 0(e) (3.9) 



This antenna function is normalised with the following two-parton matrix element 
squared: 

A° 2 (l Q , 2 q ) = 4 (1 - e) [£ c 2 m + m£] . (3.10) 

It is worth noting that the expression of this final-final flavour violating antenna Aq g - given 
in eq.(3.9) differs significantly from the flavour conserving antenna A^ g Q given above in 
eq.(3.2) as expected. 

3.3 Massive NLO initial-final antenna functions 

To construct our subtraction terms for heavy quark pair production in association with 
jets at NLO, we will also need massive initial-final antennae of different types. Flavour- 
conserving and flavour-violating massive quark-antiquark A- Type antennae and flavour 
conserving massive quark-gluon antennae are required. 

In principle, the massive initial-final massive antennae can be obtained from the cor- 
responding expressions for the massive final-final antennae given above by appropriately 
crossing one massless parton from the final to the initial state. By this crossing procedure, 
the presence of an overall uneven number of fermions crossed to define the initial-final 
antenna leads to an overall minus sign in the definition of the antennae. Furthermore, this 
crossing procedure is not non ambiguous for the quark-gluon D-Type antennae initiated 
by a gluon. 

In general, the initial-final antennae Xijk are normalised to the reduced two-parton 
matrix element squared K \ 2 to which the three-parton matrix element squared - fc | 2 
tends in the limits. For the D-Type antennae initiated by a gluon, depending which limit 
is considered, the reduced two-parton matrix element to which the three-parton matrix 
element collapses to, can be different. In particular, the nature of the initial state parton 
in the three and two-parton matrix elements squared needed as ratio to define the antennae 
may change. 

The generic process necessary to define the massive initial-final three parton antenna 
Xfj k is given by q + pi — > kj + kk with q 2 < and q 2 = —Q 2 . In this process, pi is 
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the momentum of the initial state radiator i whereas kj and k^ are the momenta of the 
unresolved parton j and the final state radiator k respectively. Depending on the situation 
considered, the one massive parton present in the final state can either be the unresolved 
parton or the final state radiator. 

In any case, the initial-final antennae -Xjjfc will depend on the invariants Sjj, sjk and 
Sik, on the masses rrij and of the final state partons j and k and on Q 2 . 

3.3.1 Quark- Antiquark antennae 

By crossing a gluon to the initial state in the expression of the final-final massive quark- 
antiquark antenna ^Q g Q one obtains the initial-final massive quark-antiquark antenna 
AP qq. In it, the gluon plays the role of the initial state radiator, the final state radi- 
ator is a massive quark or an antiquark (by symmetry arguments). The unresolved parton 
is correspondingly either a massive antiquark or a quark which can become quasi-collinear 
to the initial state gluon. This antenna is given by, 



4°tt -1 9 ^ 1 ( H2 ,2512, 
A 3 (d g , 1q,2q) = -— — 1 h 

v [Q 2 - 2m^\ V S13S23 S13 



2s\ 2 2si2 , 2S12 S13 «23 



«23 «23 «13 



2 / 8si2 2s 12 2s 12 2s 13 2s 2 3 , 2 2 \ 
-m 5 7> 1 5— H 5— H 1 (3-11; 

V V«13S23 Sf 3 S* 3 S% 3 Sf 3 S13 S 2 3 / 



This function has been normalised to the two-particle matrix element related to the 
process j*Q —> Q. It corresponds to the matrix element in which the process 7*5 — > QQ 
reduces to in all its limits. Omitting couplings, this normalisation factor is given by 

A° 2 (l Q ;2 Q ) = 4[(l-e)Q 2 -2m 2 Q }. (3.12) 

Note that the resulting antenna has an overall minus sign made explicit in eq.(3.11) 
due to the uneven number of fermions crossed to define it. 

3.3.2 Quark-gluon antennae 

As in the massless case, the quark-gluon final-final massive antennae are separated into 
two categories depending if the final state radiator gluon splits into a quark-antiquark pair 
(E-Type) or into two gluons (D-Type). These initial-final massive antennae depend fur- 
thermore on the mass m x of the decaying neutralino with momentum q. Following our 
definition of Q 2 , it is given by m x = \J —Q 2 . 

A)E-Type antennae 

Only one case of E-Type initial-final massive antenna functions is required: E 3 (4 9 ; 3 q , 1q) 
which has a massless initial state radiator 4 q , a massive final state radiator denoted by 1q 
and an unresolved parton which is the final state quark 3 q of the same flavour as the initial 
state quark \ q . This antenna accounts for the massless initial state collinear behaviour and 
is given by, 
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E° 3 (4 q ; 3„ 1q) = -2 f-si4 + S13 - — - — - 2m Q m x ) + 0(e). (3.13) 

(Q 2 + m 2 Q ) V *34 «34 7 

It is normalised to a gluon initiated process, namely X9 ~~ ^ Q- The corresponding two- 
particle matrix element squared omitting couplings reads, 

£ 2 °(2 9 ;l Q ) = 4(l-e) {Q 2 + m 2 Q f . (3.14) 

As a consequence of the crossing procedure, the antenna defined in eq.(3.13) has an 
overall minus sign made explicit in that expression. 

B)D-Type antennae 

The massive final-final D-Type antenna denoted by D®(1q, 3 g , 4 g ) has only one massive 
final state particle, the quark denoted by 1q. The massive initial-final D-Type antenna 
is in principle obtained by crossing one of the two gluons in this function to the inital 
state. However, since the initial state gluon can split either into a quark-antiquark pair 
or into two gluons, two possible reduced matrix elements can serve as normalization for 
this antenna. A simple crossing of a gluon from the final-final quark-gluon antenna is 
not sufficient to define it unambiguously. This ambiguity requires the decomposition of the 
gluon-initiated D-type antenna function into sub-antennae according to the reduced matrix 
elements it factorises to in the different limits. According to the limit considered, those 
reduced matrix elements are related to the process x + 9 — ► Q or to the process \ + Q ~ > 9- 
These two different reduced matrix elements define the two different limiting behaviours 
in which the antenna needs to be decomposed. The decomposition can be achieved by 
separating the terms in the crossed function according to their contributions to a given 
limiting behaviour. For those terms which give contributions to more than one limiting 
behaviour, partial fractioning is applied. The antenna corresponding to a reduced matrix 
element initiated by a quark will be denoted by D® Q g and the antenna corresponding to a 
reduced matrix element initiated by a gluon by D®. g Q. 

For D®.Q g , the final state hard radiator is a gluon and the parton becoming unresolved 
is a massive quark Q which can become quasi-collinear to the initial state gluon. This 
antenna is given by, 
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(3.16) 



This antenna is normalised to the matrix element associated to the process \Q ~^ 9 with 
a massive quark in the initial state. Omitting couplings it reads, 

D 2 °(l Q ;2 g ) = 4(l- e ) {Q 2 + m 2 Q )\ (3.17) 

As this D-Type antenna is initiated by a gluon but normalised by a quark-initiated 
process, it has an overall minus sign made explicit in eq.(3.16). 

For Dg g g , the final state hard radiator is a massive quark and the parton becoming 
unresolved is a gluon which can become soft, collinear to the initial state gluon or quasi- 
collinear to the final-state massive quark. This three-parton antenna reads, 
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This function is normalised to the matrix element associated to the gluon-initiated 
process \9 ~^ Q with the corresponding two-particle matrix element squared omitting 
couplings reading, 

^2 



£>°(l 5 ;2 g )=4(l-e) (Q 2 + m 2 Q Y 



(3.20) 



Being initiated and normalised by a gluon, this antenna has no overall minus sign. 



3.3.3 Flavour violating antennae 

Finally, as for the final-final case, in addition to the flavour conserving antennae given 
above, for the construction of the subtraction terms for QQ +jets production at NLO, we 
also need two types of massive flavour violating initial-final antennae. Those are such that 
in both cases a massless quark is playing the role of the initial state radiator; the final state 
radiators can either be a massive quark or a massive antiquark. The expressions of these 
two antennae are related by the exchange of a massive quark versus a massive antiquark 
in the final state such that only one expression is given below. The unintegrated form of 
the antenna having a massive quark and a gluon in the final state given by A® g q is, 
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+ C(e). (3.21) 



This antenna function is normalised with the following two-parton matrix element 
squared: 

A° 2 (2 q ;l Q )=4(l-e)[Q 2 -m 2 Q ] (3.22) 
and has therefore no overall minus sign. 

4. Singular limits of the massive antennae 

The factorisation properties of tree-level QCD squared matrix elements for massless partons 
are well known [71]. An (m + l)-parton squared amplitude factorises into a product of a 
reduced m-parton squared matrix element and a soft eikonal factor in the soft limit, or an 
Altarelli-Parisi splitting function in the collinear case. If massive partons are involved the 
factorisation still takes place but the collinear and soft behaviour have to be generalised 
to take the mass effects into account. Those generalised limits will be described both in 
final-final and initial-final configurations below. 

In the following, we shall here first recall all unresolved factors (massless and massive) 
in the final-final and initial-final configurations before tabulating the limits of all antenna 
functions encountered in Section 3. 

4.1 Single unresolved massless factors 

When only massless partons are involved, when a gluon j emitted between two massless 
hard radiators i and k becomes soft, the squared matrix element factorises and the eikonal 
factor that factorises off the squared matrix element is 

S ijk = (4.1) 

When two massless partons become collinear, the matrix element factorises yielding spe- 
cific Altarelli-Parisi splitting functions corresponding to a particular parton-parton split- 
ting. Those functions depend on z, the fractional momentum carried by the unresolved 
parton. Depending whether the unresolved parton is collinear to an initial or to a final 
state parton, the definition of z will be different. For two final state particles i and j of 
momenta pi and pj becoming collinear, we have, in the limit, 

Pi zpij, pj (1 - z)pij, s ik zs ijk , s jk (1 - z)s ijk , (4.2) 

whereas for a final state particle j of momentum pj becoming collinear with an initial state 
parton i of momentum p, L we have 

Pj ->• zpu Pij^{ l ~ z )Ph Sik ->• y~~;> s i k ~* y~; ■ ( 4 - 3 ) 

The splitting functions denoted by Pij^^(z) corresponding to the collinear limit of 
two final state partons i and j are given in [72] by, 

r \ 1 + {1 ~ z) 2 - ez 2 
P qg ^ Q {z) = ^ -!■ (4.4) 



- 23 - 



jw*) - " 2+( ;:f' e <4.5) 



2 1 — 2 

+ + z(l - z) 



1 - Z 



(4.6) 



When the collinearity arises between an initial % and a final state parton j, the splitting 
functions denoted by Pij^^(z) are given in [72] by, 
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The additional factors (1 — e) and 1/(1 — e) account for the different number of polarizations 
of quark and gluons in the cases in which the particle entering the hard processes changes 
its type. 

In all splitting functions defined above, the label q can stand for a massless quark or 
an antiquark since charge conjugation implies that P qg ->Q = Pqg^Q an d Pqg^Q = Pqg*-Q ■ 
The labels Q and G in those denote the parent parton of the two collinear partons, which 
is massless. 



4.2 Single unresolved massive factors 

When the final state partons are massive, the emission of extra radiation from those can 
still lead to soft divergences, but not to collinear singularities since the mass of the final 
state parton regulates those. The relation between matrix element squared and the splitting 
functions needs to be extended from massless to massive. Similar factorisation formulae for 
matrix elements as in the massless case hold provided the collinear limit is generalized [61] 
to the quasi- collinear limit. 

4.2.1 Quasi-collinear limit in final- final and initial- final configurations 

In the final-final configuration, two final state massive partons can become quasi-collinear 
to each other resulting in a parent parton which is massive. The limit when a massive 
parton (ij) of momentum puj) and mass muj) decays quasi-collinearly into two massive 
partons i and j of masses mi and mj is defined by, 

(1 - (4.H) 
Ph = m h (4-12) 

with the constraints 

Pi -pj, mi, mj, m ( y) ->• 0, (4.13) 
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at fixed ratios 



< ™) m h (414) 



Pi ■ Pj Pi ■ Pj Pi ■ Pj 

The difference obtained between taking the quasi- collinear limit between two final state 
particles or taking this limit when one initial and one final state particles are involved is 
closely related to the difference obtained in these two situations for the massless collinear 
limit. When a massive parton of momentum pj becomes quasi-collinear to an initial state 
massless parton pi we have, 

Pj -> m, P(ij) -> (1 - z)Pi, ( 4 -!5) 

with the constraints, 

Pi • Pj, rrij, m^j) ->• 0, (4.16) 

at fixed ratios 



(4.17) 
Pi ■ Pj Pi ■ Pj 

Only the fractional momentum z carried by the unresolved parton needs to be defined 
accordingly in both final-final and initial-final situations. For the quasi-collinear limits, 
it is defined exactly as for the collinear limits. For the final-final case, z is defined as in 
eq.(4.2) whereas in the initial-final case it is defined by eq.(4.15). 

The key difference between the massless collinear limit and the quasi-collinear limit is 
given by the constraint that the on-shell masses squared of the final state partons have to 
be kept of the same order as the invariant Sij = 2pi • pj, with the latter becoming small. 

4.2.2 Factorisation in the quasi-collinear limits 

In these quasi-collinear limits (final-final or initial-final), the (m + l)-parton matrix element 
squared factorises into a reduced m-parton matrix element and unresolved massive factors. 
These single unresolved massive factors are generalizations of the massless unresolved fac- 
tors defined above. 

The generalized soft eikonal factor Sijk(rrii, rrik) for a massless gluon j emitted between 
two massive partons i and k depends on the invariants S[ m = 2pi -p m built with the partons 
i, j and k but also on the masses m 8 and of partons i and k. It is given by [61,64] 

C ( \ 2Sik 2m ^ 2m k (A 1 

S ijk {mi, m k ) = ^- 2"*. ( 4 - 18 ) 

SijSjk Sj k 

The massive splitting functions will depend on z, the fractional momentum carried by 
the unresolved parton j and on the masses and rrij of the partons i and j becoming 
quasi-collinear. All the mass dependence can be parametrized by fj,^ given by, 

mf + m 2 - 
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The massive splitting functions denoted by P^^^j^z, fifj) for a massive parton (ij) 
which splits into partons i and j, both being in the final state, have been given in the 
appendix of [62] and in [64]. Those read, 

= 1 + (1 "5" £Z2 - 2/4 (4.20) 

^2 , f-i „\2 _ i ,,2 



z 2 + (1 - zY - e + ul n 

p q ^ G (z,^) = — — — > ; '- 2| ; 

where 



= — and (4-22) 

Naturally, the gluon-gluon splitting function P gg ^c( z ) is left unchanged. The splitting 
function P qq ^c{z, l&q) given in eq.(4.21) and related to the quasi-collinear limit of two 
massive partons does not correspond to a limiting behaviour of the antenna functions 
required here and given in Section 3. It is given here for completeness. 

So far we have treated all generalisations of soft and collinear massless factors needed 
to treat final-final configurations involving massive partons. For the initial-final situations, 
since all the initial state partons are taken massless the only initial-final splitting function 
that changes when we allow massive partons in the final state is P qq <^G{z, H qg ) given by, 

P qq ^ G (z, fa) = " 2 + ( 1 1 :" )2 ~ £ + 2/4. (4.23) 

The definition of the momentum fraction z present in this formula will be the same as in 
the massless case given by eq.(4.3). 

4.3 Singular limits of the massive antenna functions 

In this subsection we list all the non- vanishing soft, collinear and quasi-collinear limits of the 
massive final-final and initial-final antenna functions given in Section 3. The limits of the 
massless antennae also needed to construct the subtraction terms for QQ+ jet production 
can be found in [42,57,60]. 

4.3.1 Final- final antenna functions 

The limits of the massive flavour conserving quark-antiquark and quark-gluon antenna 
functions have been derived in [64]. We give them here for completeness. The limits of the 
flavour violating antennae are new. 

The singular limits of the massive quark-antiquark antenna A Q g Q are 

3 o ->-0 



a °0-Q^ 2 q) — — ► Si32(m Q ,m Q ), 



A l( l Q^9i 2 Q) Q) — P qg^Q( z -,V 2 qg), (4.24) 

^(1q,3 9 ,2q) -^-P qg ^ Q (z,fa). 
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The only non-vanishing singular limit of the quark-gluon E-Type antenna Eq g ,-, is the 
collinear massless limit of the massless quark-antiquark pair, 

£ 3 °(1q, 3^,4,-0 — P OT -_> G (z), (4.25) 

•S34 

while for the D-type Dg gg and d-type d°Q gg antennae we have, 
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For the massive A- Type flavour violating antenna Aq ^ we have 

A° 3 (l Q ,3 g ,2 g ) ^ S 132 (m Q ,0) (4.27) 
A° 3 (l Q ,3 g ,2 g ) ^? — P m -q(*,/4) (4-28) 

Sl3 

^(l Q ,3 fl ,2 ? -) ^ -^-Vq(z). (4.29) 

S23 

Again, for symmetry reasons the role of Q and q are interchangeable in these formulae. 
4.3.2 Initial-final antenna functions 

For the quark-antiquark initial-final massive antenna A Q qq, only quasi-collinear limits are 
present. Those are: 

A° 3 (3 g ;l Q ,2 Q ) ^? — P M -«_ G (z,^) (4.30) 
S13 

A° 3 (3 g ;l Q ,2 Q ) ^ ±P qq <-G(z,v? qg ) (4.31) 
S23 
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For the E-Type antenna E q . q Q, the only non- vanishing singular limit is, 

E° 3 (4 q ;3 q ,l Q ) — P OT <_ (z). (4.32) 

The limits of the two quark-gluon D-Type antennae initiated by a gluon D g ^Q g and 
Dg :9 Q are given by : 

D° 3 (4 g ;l Q ,3 g ) * ^-P q ^ G (z,n 2 Qg ), (4.33) 

which has only a quasi-collinear limit when the unresolved quark 1q is collinear to the 
initial-state gluon. The limits of D 3 (4 g ; 3 g , 1q) with an unresolved gluon are 

D° 3 (4 g ; 3 g , 1 Q ) ^? —P qg ^ Q (z, ii\ g ) (4.34) 

D° 3 (4 g ; 3 g , 1 Q ) ^ —P gg ^a(z) (4.35) 

D° 3 (4 g ; 3 g , l Q ) S 134 (m Q , 0). (4.36) 
Finally, the initial-final flavour violating antenna A q g Q has its limits given by, 

A° 3 (2 q ;3 g ,l Q ) ^? —P^Qiz,^) (4.37) 

A° 3 (2 q ;3 g ,l Q ) ^ -Vq(^) (4-38) 
S23 

A° 3 (2 q ; 3 g , l Q ) ^ S l32 (m Q , 0). (4.39) 
5. Integrated massive antenna functions 

To combine the antenna subtraction terms with the virtual corrections and the mass fac- 
torization counterterms in a given kinematical configuration, the antenna functions have 
to be integrated over the appropriate factorised antenna phase space. After integration, 
the implicit soft and collinear singularities present in the antenna functions turn into ex- 
plicit poles in the dimensional regularization parameter e, and the remaining phase space 
corresponds to the same n-particle kinematics as the virtual contributions or the mass 
factorisation counterterms. In this section, we derive the integrated forms of the massive 
antenna functions defined in Section 3. Only full antennae, denoted by capital letters 
(or crossing of those), need to be integrated, while partial antennae denoted Xijk sum up 
to X^k prior to integration. 

5.1 Properties of the integrated massive antennae 

The results will be presented in two forms, in expanded and unexpanded forms in the di- 
mensional regularization parameter e. The unexpanded forms of the integrated antennae 
are functions of a few master integrals obtained after standard reduction techniques [73-77] 
have been applied. The masters will be given here analytically to all orders in e. Addition- 
ally, the integrated antennae will be presented after an e-expansion has been performed 
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on these all-order results up to finite order in e. In this expanded form, the poles in e 
become explicit and can be related to process independent infrared singularity operators 
and splitting kernels. 

Final-final integrated antennae have their pole part entirely related to colour ordered 
massive infrared singularity operators lQ\ which will be defined below. For initial-final 
antennae, those infrared operators are not sufficient to capture all singularities present in 
the integrated antennae. Additional pole terms can arise due to the presence of massless 
intial-final parton-parton collinear singularities which cancel against the mass factorisation 
counterterms. These pole terms are proportional to universal and process- independent 
splitting kernels p^\x). Those are defined for example in [57] and will be given below. 

As a check on our results for the expanded forms of the integrated antennae we consider 
the following: At NLO, a particular antenna can be regarded as the sum of two particular 
dipoles in the dipole formalism of [30,61]. The two radiators present in an given antenna 
play then both the role of emitter and spectator in the corresponding two dipoles. Up 
to terms which do not give rise to singularities when integrated over the phase space, the 
sum of these two dipoles and the antenna are the same (up to coupling factors). The pole 
parts of a given integrated antenna and those of the corresponding integrated dipoles can 
therefore be related. For each of the integrated antenna we will specify how this comparison 
is performed. 

5.1.1 Infrared singularity structure 

Owing to the universal factorisation properties of QCD amplitudes in infrared singular 
limits, it is possible to describe the infrared pole structure of virtual loop corrections and, 
consequently, also of integrated subtraction terms, by the product of infrared singularity 
operators with tree level amplitudes. These infrared singularity operators are a priori 
tensors in colour space. In a colour ordered framework, they decompose into different 
colour-ordered infrared singularity operators. In massless QCD, there is solid evidence to 
assume that the infrared singularity operators consist only of combinations of two-particle 
correlations [78-80] at all orders in perturbation theory. Explicit forms of the massless 
operators are known to three-loop order [81]. In massive QCD, only the one-loop infrared 
singularity operator is made up entirely of two-parton correlations [62], while multi-particle 
correlations can contribute at higher loop order [82]. The explicit form of the massive 
infrared singularity operators is known to two loops [82] and was used to predict the pole 
structure of the two-loop matrix elements for qq — > tt and gg — > tt. 

We are concerned only with the infrared singularity structure at NLO in the present 
study. Consequently, the integrated massive antenna functions will contain a pole structure 
in terms of the massive one-loop infrared singularity operators. Containing only two-parton 
correlations irrespective of the particle masses involved, these can be expressed straightfor- 
wardly in a colour-ordered form. We introduce the mass-dependent colour-ordered NLO 
real radiation singularity operator 

lf)(e,Sij,mi,mj, Xfi 2 ) , (5.1) 
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which describes the unresolved real radiation between partons i and j. It is a function 
of the invariant mass of the parton pair, of the masses of the partons and of a product 
of kinematical parameters X/j, 2 , which determines the logarithmic pole coefficient. In this 
form, /j, 2 is dependent only on the mass combination of the radiator partons, while A takes 
account of the nature of the kinematically allowed endpoint in the different kinematical 
configurations. A = 1 for all final-final antenna functions and for initial-final antenna 
functions with equal masses, while A = Xq with 

Q 2 



Xq 



Q 2 + m 2 Q 



(5.2) 



for initial-final antenna functions with one massless and one massive radiator. 

In their most general form those infrared operators have their functional dependence 
given by, 

I QQ( e ' s 0f' m 0' m Q'^ 2 )' O e ' s Q9' m Q' m Q' V 2 ) and I^ ^e, Sij,m Q ,0, fi 2 ). 

The third operator describes contributions arising from the splitting of a gluon into a 
quark-antiquark pair which are proportional to the number of light quark flavours Np. 

The massive infrared operators which have a non-trivial mass dependence, required in 
both final-final and intial-final configurations are, 
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(5.3) 
(5.4) 



The antiquark-gluon operators are obtained by charge conjugation: 



r(l) 



Iq{z, s g Q,m Q , 0, A/x 2 ) = lJ2(e, s gQ , m Q , 0, A//), 

^gQ,F^ s <?<2> m Q' °' ^) = ^gQ,F( e i s gQi m Q, °> ^) ■ 
In addition to these infrared operators which have their massless counterparts defined 
in [42], we also have an infrared operator associated with the flavour- violating antennae 
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AQgq and AQ gq denned in Section 3. As those antennae have the same unintegrated form 
both of them are related to one infrared operator denoted by 1^, where Q can stand 

for a massive quark or antiquark. As the quark-gluon operator Iq^, this operator has 
a mass-dependent logarithmic term proportional to X/j 2 with A = Xq in the initial-final 
configuration. It is given by, 
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2T(l-e) 



s q Q + m % 
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5.2 Integrated massive final-final antennae 

The integrated massive final-final flavour conserving antennae have been derived as function 
of a few master integrals in [64]. In this section, all final-final integrated antennae will be 
presented in expanded form such that their pole part can be related to the colour-ordered 
massive singularity operators 1^ defined above. Before we shall derive the integrated 
massive flavour-violating antenna and give its expanded and unexpanded forms. 

In general, the integrated final-final antennae denoted by Xij k are given as the integra- 
tion over the final-final antenna phase space d$x ijfe of the unintegrated final-final antennae 
X^k as given in eq.(2.7). Those will depend on the masses of the final state particles and 
on E 2 m . 

We start by giving the NLO massive final-final antenna phase space d<&x ijk necessary to 
evaluate these integrated final-final antennae. In the most general case, where the involved 
partons i, j, k have three different masses rrii, rrij and m k , the massive antenna phase space 
d^' k mj ' mk) is given by [63,64], 
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3-d 
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dsij ds jk ds ik 5{E 2 cm - m 2 - m 2 - m 2 k - Sij - s jk - s ik ) 



[4:A 3 ( Pi , Pj ,p k )] 2 e(A 3 (p i ,p j ,p k )). 

The masses mj and nix appearing in this equation are combinations of the masses mj, mj 
and m k . The function A 3 (pi,pj,p k ) is the Gram determinant for massive particles of 
momenta Pi,Pj,p k given in terms of invariants = 2pi ■ pj and masses m,, rrij, m k by, 

A 3 (pi, Pj ,p k ) = - (sijS ik s jk - m 2 s 2 k - m\s 2 j - m 2 s 2 k + Am 2 m 2 rn 2 k ) . (5.7) 

To be able to evaluate the integrated flavour-violating massive final-final antennae we 
need to consider the case where only one particle in the final state is massive while the 
other two are taken massless. We consider, mi = rrij = and m k = iuq, in which case 
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the masses of the remapped momenta pi and px have their masses given by mj = nii = 0, 
rax = TUk- We use the following parametrisation of the final- final massive antenna phase 
space with one massive parton d$^'°' m ' ) given by, 

ijk 
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d **«* - r(l -e) ^ ( o) 
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s ife = ^cm u o(l - ■")• (5.9) 

Integrating the flavour violating antenna A®(1q, 3 s , 2g) given in eq.(3.9) over this phase 
space d$x' '™^ we obtain, 

,o n o 2 v _ 8^ 2 (4vr)-e^ 

2 - E cm^o( 1 - n o)e 2 (l-2e) 

x|^ m [12(1 - u ) 2 - 2e(29 - 54u + 24u§) + e 2 (98 - 176u + 73u 2 ) 

+e 3 (-68 + 120u - 47uq) + 2e 4 (8 - 12u + 3ujj)] jf'°' m) (5.10) 

+3(1 - e) [4(1 - n ) + 2e(-7 + 6u ) + e 2 (14 - llu ) + 2e 3 (-2 + u )] 4°'°* m) j . 

In this expression, the master integrals required are I^°'°' m ^ and /^ ' 0,m ) were derived in [64]. 
The master integral /j ' 0,m ) corresponding to the integrated 1 — > 3 phase space measure 
with one massive final state is given by, 

7 (0,0,m) = /" d$ (0,0,m) 
1 / Xijk 



/ p 2 ^-e„,2-2e -4+2e_-2+e r ( 2 2g ) r (l 6 ) p /-, r, 9 . . s 

\&cm) u 2 * r ^ _ 4 ^ 2^1(1-6,2-26,4-46,^0), 



(5.11) 



while I^°'°' m ^ i s given by, 



r (0,0,m) 
J 2 



_ / p 2 ^2-e„,3-2eo-5+2e_-2+e ^(2 26 ) r (l <Q 

~ \^cm) u 2 n _ 4 ^ 2^1(1-6,2-26,5-46,^0), 

(5.12) 
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Expanding the integrated flavour- violating A- Type antenna -4,g 9 g given in eq. (5.10) in 
powers of e up to finite order, we obtain 



AQgqO-QiZgilq) = ~2l£J I €, SQgq, m Q , 0, 



m 



Q 



-) 



m Q + S Q 9 qJ 



1 19 5vr 2 

i-/x 2 ~ T + "12" 



l 



(73^ - i) HV 2 ) + ^V) + Li 2 (l - 



where 



2_ m Q 



EL' 



+ 0(e), (5.13) 



(5.14) 



All other integrated flavour conserving final-final antennae given in terms of masters 
in [64] can be written in terms of the real infrared singularity operators 1^ defined above 
up to finite order as follows, 

_ _ 9T (i) f, „ _ ^ „ 1 " 

40 ln(2) (r - 3) (1 + r) tanh" 1 ( y/r) 



A^Q^g^q) = -^-qq ^ s QgQi m Qi m QiYT^ 
1 [ 1 

+2^-39 + 17r - 161n(2)(r - 3)) 



+(-33 + lOr - r 2 + 161n(2)(r - 3)(1 + r)) In 



1 + y/r 
1-a/F 



+4(1 + r 



-16(r - 3) tanh _1 (l - 2r) ^-2^ + (1 + r) In Q + ^ 
+ 10(r - 3)(1 + r) (In 2 (l + y/7) - In 2 (l - y/f)) 

i + Jr 



(5.15) 



5Lio 



5Li 2 ( X —^- ) + 3Li 2 (r) - 12l.i, ( N <7) 



m 



Q 



m Q + S Q?Y 



6(1 -/x 2 ) 3 
-2// 3 (-3-3// + /x 3 )m(// 2 



6 + 3/i - 14^ 2 + 14/i 4 - 3^ 5 - 14,u 2 + 14/x 4 - 3,u 5 - 6^ 6 



+ 0(e), 



(5.16) 



P 3 °(1q,3 9 ,4 9 ) = -4lg> le,s Qgg ,m Q A-^r— + 



2^3 



- 2(1 - ^ 2 )(4(-14 + vr 2 ) - 3/i + /i 2 (126 - 8vr 2 ) + 9/x 3 + /i 4 (-74 + 4vr 2 )) 
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-2(-3 + 33/T + 6// - 51^ 4 + 17/i b ) ln(//) + 3(1 - /u) 2 (-2 - 4/i - ^ + /T) ln 2 (^) 
+12(1 - ^) 2 (-2 - 4/x - /x 2 + M 3 )Li2(l - ^ 2 ) + 0(e), 

(5-17) 

with r and /x denned as in eq. (5.4) and eq. (5.14) respectively. 

Let us notice that the last term present in the expansion of Aq g Q arises through the 
expansion at finite order of 

+2e 



1 /(l + r)\ . fl-y/r 
In 



2m Q 

with E 2 m = Sq 9 q + 2mg. This term arises since we have chosen to factor 



s 2 - 
QQ 



in ^qq as unexpanded overall factor whereas the integrated antenna -4,q 9 q is naturally 



proportional to 



*QgQ + 2m Q 



We have chosen to define the 1^ operator such, in 



order to have similar overall factors in all I-j^-type operators. 

From these results on the expanded integrated final-final antennae, we see that the 
pole parts of all of those can be captured solely with poles present in the massive colour- 
orderered I-j^ operators defined above. 

Finally, for all integrated final-final antennae, we have compared the pure pole parts of 
those with the pure pole parts of the corresponding sums of the integrated dipoles presented 
in [63]. Setting couplings and colour factors to one we found full agreement. 

More explicitly, the poles of Aq 9 q can be compared with the poles of two dipoles having 
both two massive final states, a spectator Q and an emitter Q. The poles of Sq^^ can 
be compared with those of two dipoles with a massless final state emitter q and a massive 
final-state spectator Q while the poles of T^Q gg have to be compared with the sum of two 
dipoles where the massive quark and the gluon in the final state can both be emitter or 
spectator. 

5.3 Integrated massive initial-final antennae 

In order to obtain the integrated massive initial- final antennae Xijk, the massive unin- 
tegrated initial-final antennae -Xjjfc defined in Section 3 need to be integrated over the 
initial-final massive antenna phase space d&x ijk as given in eq.(2.16). As a result, the 
integrated antennae will all depend on Q 2 = —q 2 , on x, the momentum fraction carried by 
the initial state parton pi , and on the masses of the final states present in a given antenna. 
As we saw in Section 2, the momentum fraction x depends on the number of massive par- 
ticles present in the final state in a given antenna. For two massive final states j and k of 
masses rrij and x is given by, 

Q 2 + m 2 + ml 

X = . 

2pi • q 
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Generally, the pole parts of the integrated massive initial-final antennae are related to 
the massive 1^ operators defined above and to the x-dependent colour ordered splitting 

kernels p^ (x) . The splitting kernels describe the initial-final massless collinear singulari- 
ties and are given by [57], 

V { S{x) = 1 6(1 - x) + 2V (x) -1-x, 
p$(x) = l-2x + 2x 2 , 

yH x 

P { aKx) = —S(l-x) + 2V (x) + --4 + 2x-2x 2 , 
MM b x 

pSf(*) = 4«(1-x), (5.18) 



99,*" ~> 3 

with the distributions 

VJx 



ln n (l - x) 



1 — x 

In the following, we shall specify the phase space parametrisations for the initial-final 
antenna phase space d<S>Xi jk needed to integrate the different antenna types before giving 
the results for the integrated initial-final massive antennae in unexpanded and expanded 
forms. 

5.3.1 Phase space parametrisations for initial- final configurations 

In the most general case, the massive two-to-two phase space d<3?2 to which the initial- 
final antenna phase space d&x tjk defined in eq. (2.16) is proportional to, is related to the 
process q+Pi — > Pj +Pk with pj and pk the momenta of the final state partons with masses 
Tfij and m k . It can be written as 

(Air) 6 r i 1 ~ 2e 

d<j) 2 {m j: m k ) = ^ _ ^ E 2e ~ 2 (E 2 cm - m 2 - m 2 k f - Am 2 m 2 k 2 dy y^{l - y)~ e , 

(5.19) 

where y runs from to 1. In this parametrisation the invariants take the following form 

Q 2 + m 2 + ml 
2x& X 

(e 2 ^ + m 2 - m\ - (2y - 1)^ (iP^ - m 2 - m{f - 4m]m^ (5.20) 
Q 2 + m 2 + ml 

\E 2 m - roj +ml + (2 v -\)^ (e?„ - rof - roj) * - 4m?mA , (5.21) 
and the partonic center of mass energy reads 



; Q 2 (l - x) +m 2 + ml 

E cm = V¥+W=\— 1 -■ (5-22) 
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In this case also x is defined as, 

Q 2 + m 2 + m\ 



x = 



In our application to ti +jet production, we are only interested in the cases where 
rrij = rrik = itiq or in the case where rrij / and = 0. For the first case, the 
parametrisation given above simplifies to 

dUm Q ,m Q ) = ffl^ cm" 2 [{ E L - 1rn 2 Q f - 4m 4 Q ] ~ dy y^{\ - y)~\ (5.23) 

while the invariants read 

Q 2 + 2m 2 n ( „ H ~V2 ~\ 

2pi " Pj = 2xEl m [ Ecm ~ (% ~ 1} V ~ Q ) ~ Q ) (5 " 24) 

and the center of mass energy is 



Q 2 (l-x) + 2m 2 Q 2 + 2rn 2 

E cm = \— ^ with x = — 2 5.26 

V x 2pi.q 

In the case in which rrik = and rrij = uiq the phase space reduces to 

dMm Q ,0) = ^Z_E 2 ^ 2 (E 2 cm - mlf-^dyy-^l - y)~\ (5.27) 



the invariants are 

2 ^ • PJ = Q T tr Q i E cm ~ V ( E cm ~ m 2 Q )] 
^2 i ,^2 



and the center of mass energy is 



(5.28) 



2 K • Pk = [V { E L ~ m 2 Q )] (5.29) 



Q 2 (l-x)+m 2 Q 2 + m 2 

E cm = \— ^ with x = — 5.30 

V x 2 Pi -q 

5.3.2 Integrated forms of the initial-final massive antennae 

The integrated forms of the antennae can be obtained using reduction techniques, using 
the extension of the integration by parts method [73, 74] in [75, 76] to reduce the phase 
space integrals to master integrals. In this task, we express all the invariants in the an- 
tenna functions as massive propagators and write the three on-shell conditions p 2 = m 2 
(a = j, k) as cut propagators. Since the invariants s a b appearing in the antennae are not all 
independent from each other, the integrated antennae are written as a one-loop diagrams 
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with on-shell conditions and not as two-loop diagrams as in the final-final case. The reduc- 
tion to master integrals is therefore easier. It was done using the Laporta algorithm [77] 
with two independent implementations: the Mathematica package FIRE [83], and an in- 
house implementation in FORM [84]. For all integrated initial-final antennae required and 
defined in Section 3, we find four master integrals. Those can be evaluated analytically in 
terms of gamma functions and hyper geometric functions. The all order expressions as well 
as the expanded expressions of all initial-final massive antenna functions will be presented 
below. 



A) Quark- Antiquark antennae 

The integrated form of the massive three-parton initial-final quark antiquark antenna given 
by Ag QQ is obtained by integrating its unintegrated form A q q defined in eq.(3.11) over 
the initial-final massive antenna phase d&x ijk given in eq.(2.16) using the parametrisation 
of the two-to-two parton phase space d<&2( m Q, m o) given in eq.(5.23). It depends on the 
virtuality of the incoming boson Q 2 , the masses of the final states both being uiq and on 

Q 2 +2m 2 Q 

x given by x = — 9 „ „ • It reads, 



^3(3 s ;1q,2q) 



+- 



_ 2(4vr) 1 - e e e T B 
~ Q 2 + 2m 2 Q 

2™ 2 Q (2-e) 
Q 2 (l - e) - 2m 2 Q 

4(2 - e) 2 m% 
Q 2 {l-e)-2m 2 Q 



+ 2x{\ - 2e) 1 



Q 2 x 



Q 2 + 2m 2 Q ^ 



+ 2m 2 Q (2- e + 2x 2 e) 



+(1 - 2x + 2x 2 )(l - e){Q 2 + 2m 2 



Q> 



I 3 (mQ,m Q ,s 13 ) 



h{m Q ,m Q ) 



(5.31) 



h(mQ, tuq) is the integrated phase space measure with two massive final states. It is given 
by 



I 2 (mQ,m Q ) = (4tt) 



-i r(l-e) 



x 



2r(2-2e) \Q 2 (\-x) + 2m 2 Q 
while the remaining master integral reads 



AmqX 



l-2e 
2 



Q 2 (l-x) + 2m^ 



(5.32) 



( d<i> 2 (m Q ,m Q ) — 

J Sij 



= (47T 



T(l-e) 



x 



T(2 - 2e) \Q 2 + 2m 2 Q J \Q 2 (1 - x) + 2m 2 Q 



,l-2e 



1 + V 



F 1 l,l-e,2-2e; 



2v 
1 + v 



(5.33) 
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where, 

-y-^i) -('- gsa-xH^ J ■ (5 - 34) 

The integrated form of the gluon-initiated quark-antiquark massive antenna ^.qq is 
finite, since it has only quasi-collinear limits. Expanding the all order result given in 
eq.(5.31) in powers of e up to finite order, we find 

A(3 9 ; lQ, 2q) = {Q2 - 2m 2 Q y {Q 2 _ 2m 2 Q) 

x I [Q A (1 -2x + 2x 2 ) + Am 2 Q Q 2 + m^(4 + 8x - 8x 2 )} In ( (5-35) 
+v [Q 4 (l -2x + 2x 2 ) + 4m Q Q 2 (l -x) + m 4 Q (4- 8x)] j + 0(e). 
B)Quark-Gluon antennae 

As only one massive parton is present in the final state, the fractional momentum x carried 
by the initial state momentum pi is defined as x = 2p . For the initial-final antenna 
phase space d$x ijk over which the antennae are integrated, we use the parametrisation 
of the two-by-two parton phase space d<l>2(mQ,0) given in eq.(5.27). After reduction, we 
found that two master integrals are needed. 

The master integral l2(m,Q,0) corresponds to the integrated phase space in eq.(5.27) 
and reads, 



7 2 (m Q ,0) = (47rr 1 LTl 



m^x 



l-2e 



(5.36) 



2r(2-2e) I Q 2 (l-x) + m 2 Q I I Q 2 (l-x)+m 2 Q 
while the remaining integral is given by 

h(mQ,0,s ik ) = I d$ 2 (^Q,0) — 

J Sik 

l r(l — e) / x \ ( x \ e 

= ^ ]1 '\v{2-2e) [WTm 2 , ) \Q 2 {l-x)+ m 2 Q ) ^ 

x u l ~ 2e 2 Fi(l,l-e,2-2e;n) 



with u being given by 

u = 1 



Q 2 (l -x)+m 2 Q 



B.l) E-Type antennae 

The integrated form of E® q Q is only proportional to the integrated phase space measure 
hijnQi^) and reads, 

(47r) 1_e e e7 - B 
f s (4,;3 s ,lg)= [Q 2 + w 2 q] 
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Jo (Q 2 (l-x)+m^)(l + (l-x) 2 ) 2 /l + x\ .1 . . 

x j 2m m x ex{l-x) + Q \l^x~) ~ Q ( " } | 2 ^ mQ ' 0) 

(5.38) 

Its e expansion is given by, 

^ 3„ 1q) = f Q- 7j [Q 2 + m 2 Q ] ~ e x | - - „£>(*) 

+ 2~(2 - 2x + x 2 ) [-2 + 21n(l - x) - ln(x) - ln(l - x x)] (5.39) 

Q 2 (-l + 3x - 2x 2 ) + m? Q (l + x) + 2m Q m x (l - x) 1 
+ 2(Q^l-x) + ml) + ° {e) ) ' 

with xq being given in eq.(5.2), and where we have kept the natural phase space factor for 
initial-final quark-gluon integrated antennae, 

unexpanded. Since this antenna has only one massless initial-final collinear limit, the pole 
part in its integrated form is only related to the x-dependent splitting kernel pfg (x) . This 
expression can be directly compared with the expression obtained in [63] for the integrated 
dipole involving a massless initial emitter q and a massive final-state spectator Q and full 
agreement is found. 

B.2) D-Type antennae 

The integrated form of the gluon-initiated initial-final antenna D® q g , which is such that 
the reduced matrix element is induced by a quark and which has a massive quark Q as 
unresolved parton, is 

^3°(4 9 ;1q,3 9 ) 



(Q 2 + m 2 Q ) (l-e)x(l-x) 
1 



2(Q 2 + m 2 Q )(Q 2 (l - x) + m 2 Q ) 

x Q 6 (l - x) 2 [-3 + 6x - 4x 2 + 4x 3 + 2e(l - 4x + 6x 2 - 6x 3 ) 

+e 2 (l + 2x-8x 2 + 8x 3 )] 
-Q 4 m 2 Q (l - x) [9 - 25x + 32x 2 - 16x 3 + 4x 4 - 2e(3 - 14x + 30x 2 

-24x 3 + 6x 4 ) - e 2 (3 + 5x - 24x 2 + 32x 3 - 8x 4 )] 
+Q 2 m i Q [-9 + 32x - 57x 2 + 38x 3 - 8x 4 + e(6 - 32x + 82x 2 - 76x 3 + 24x 4 ) 

+e 2 (3 + 4x - 25x 2 + 38x 3 - 16x 4 )] 
+m 6 Q [-3 + lOx - 19x 2 + 4x 3 + 2e(l - 5x + 12x 2 - 6x 3 ) 

+e 2 (l + 2x-7x 2 + 8x 3 )] (5.40) 
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+2e(l - e)Q 4 m Q m x (-l - x + 2x 2 ) 
—AQ 2 m^m x x [— x(l — x) + e(l — e)(l + x — x 2 



+2m^m x x [2x - e(l + x) + e 2 (l + x)] 



(1 - e) 2 Q 4 x(l - 3x + 4x 2 - 2x 3 ) 



h(m Q ,0) 



+2Q 2 m 2 Q [1 - 4x + 4x 2 - 2x 3 + e(-2 + 7x - 6x 2 + 2x 3 ) + e 2 (l - 2x + 2x 2 )] 
+m Q [1 - 7x + 4x 2 - 2x 3 + e(-2 + 8x - 4x 2 ) + e 2 (l - x + 2x 3 )] 



— 2e(l — e)Q 2 mQm x x + 2(1 — e + e 2 )m^m x x 
Its expansion up to finite order in e reads, 

^3°(4 9 ;i Q ,3 9 : 



4(Q 2 + m 2 Q f (1 - x 
(Q 2 + m 2 ? )(l-x) ^ 6/i ^ 2 
(Q 2 (l - x) +m 2 Q ) 2 x 

+32x 2 - 16x 3 + 4x 4 ) + Q 2 m Q (-9 + 32x - 57x 2 + 38x 3 - 8x 4 ) 



- x) 2 (-3 + 6x - 4x 2 + 4x 3 ) - Q 4 m 2 Q (l - x)(9 - 25x 



+m Q (-3 + lOx - 19x z + 4x d ) + kQ 2 u? Q m x x l {\ - x) + 4to q 



m v x 



-2 In 



m Q X 



Q 2 (l -x) + 



Q 4 (l - x)(l - 2x + 2x 2 ) - 2Q^m Q (-l + 4x - 4x 2 + 2x 3 ) 



-mg(-l + 7x - 4x 2 + x 3 ) + 2m Q m x 



x 



+ 0(e) 



(5.41) 



This integrated antenna is finite since its unintegrated form has only quasi-collinear limits. 

Finally, the integrated form of the initial-final massive D-Type quark-gluon antenna 
D® g Q which is initiated by a gluon, which has a reduced matrix element induced by a 
gluon and where the final state gluon is the unresolved parton, is given by, 



^3°(4 s ;3 3 ,1q) 



(47T) 



Q 2 (l - e) 2 + 2m Q m x e{\ - e) + m Q (3 - 4e - e 2 
(l-e)[Q 2 (l-x) + m 2 3 ] 

4(1 - 2e)[Q 2 (l - x) + m 2 Q ] (l - x + x 2 ) 2 
ex(l-x) 2 (Q 2 + m 2 3 ) 
2 



+ (l- g )(l-x) 2 (Q 2 + m Q ) h Qmx(1 " 6 + 6 )(1 " XY 
+Q 2 {1 - x)[-5 + 3x + e(8 - 6x) - e 2 (l - x)] 

+m 2 Q [-9 + 16x - 9x 2 + 2e(4 - 7x + 4x 2 ) + e 2 (l - x) 2 ] 



(5.42) 



>I 2 (m Q ,0) 
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The e expansion of this integrated antenna reads 

P 3 °(4 9 ;3 9 ,1q) = -2lg> (e,Q 2 ,m Q ,0,x 2 Q2 m + Q m2 ^ j 5(1 - x) - ^(x) 
3 / ^ 



--+\l--)S(l-x)- V (x) + 2V 1 (x) 



+ 



Q 4 (l - x) - 2Q 2 m 2 Q {\ - x){2 - 9x) - 5m J (1 - 3x) 

(5.43) 



4(Q 2 (l-x)+m 2 ) 2 

+4Q 2 m Q m x (l - x) 2 + 4niQm x (l - 



+2 ^-2 + i + x - x 2 ^) ln(l - x) + Q<5(1 - x) + y-^-j- - ^o(x)j M 1 ~ x o) 

+ ^(1 " x) ln(x 2 ) + ln(x 2 ) ln(l - x ) + \ ln 2 (x 2 ) + 1 ln 2 (l - x ) 

(1-x + x 2 ) 2 . . , . N1 1 mi/ n, /xn(Q 2 + m 2 )\ . 
~ x(1 _ x) H*) + ln (! " ^oar)] - -p$(x) In (^^^4 "J + °( e )' 

where xo is given in eq. (5.2). 

From this expansion it can be seen that all the pole parts are contained in the infrared 
singularity operator 

lg]^,Q 2 ,rn Q ,0,x 2 ^-j^ 

as well as in the splitting kernel Pgg(x). A factor Xq appears explicitly in the mass- 
dependent logarithmic term in the infrared operator 1^ . It is necessary in order to capture 
all poles parts proportional to 5(1 — x) of the integrated antenna, V®(4 g ; 3 9 ,1q). 

Furthermore, we can compare this expanded result for V g g q with the corresponding 
integrated dipoles in [63]. The pole parts of our expression given in eq.(5.43) corresponds 
to the pole part of the sum of two dipoles which have either a massless initial spectator and 
a massive final-state emitter or a massless initial state emitter and a massive final state 
spectator. Full agreement is found providing us with a strong check on our result for V° g q. 

C)Flavour violating antennae 

In addition to the flavour conserving integrated inital-final antennae, we also need to con- 
sider the integration over the massive initial-final antenna phase space of the unintegrated 
initial- final flavour- violating antenna A®. g Q. The phase space parametrisation of the initial- 
final antenna phase required here is the same parametrisation as the one used for integrat- 
ing the initial-final quark-gluon antennae. One uses the parametrisation of the two-to-two 
particle phase space given by d$2(^Q, 0) in eq.(5.27). The integrated flavour violating 
antenna ^q- g Q can be written in term of the phase space measure l2(mQ,0) only. In term 
of this master integral, its integrated form reads, 

(4 7r ) 1 - e e e 7E 



2e(Q 2 + m 2 Q )(Q 2 (l-x) + m 2 Q )(l 



x) 2 
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{Q 4 (l - x) 2 [2(1 + x 2 ) - e(7 - 8x + 6x 2 ) + e 2 (l - 2x) 2 ] (5.44) 
+2m 2 Q Q 2 {\ - x) [2(1 + x 2 ) - e(7 - 9x + 6x 2 ) + e 2 (l - 3x + 2x 2 )] 
[2(1 + x 2 ) - e(7 - lOx + 7x 2 ) + e 2 (l - x) 2 ] } I 2 (m Q , 0) 



Expanding in powers of e we obtain 

/ 2 



^(2,;1 ,3 9 ) = -21^ e ,Q2 )mQ)0)X 2 



^(l-x)-lpW(x) 



4 + f 1 " SI <^ - *) - ^o(x) + 2Pi(x) - £ - 



2 V 12 / 2 4Q 2 (l-x x) 2 

+ 4 ( i-x )(i-x x) + (s'P - " + rb - 111(1 - *»> (5 - 45) 

+^5(1 - x) (5 ln(x§) + 2 In(xg) ln(l - x ) + ln 2 (xg) + ln 2 (l - x )) 

-(1 + x) ln(l - x) - 6( |^ 2) [ln(x) - ln(l - x x)] 

1 m , „ / xl(Q 2 + m 2 )\ 
-- P W(x)\n (^^^4 -) 

All the pole pieces are contained in the following infrared singularity operator 

and the splitting kernel pqq (x) . As for the integrated gluon- initiated quark-gluon antenna 
^ssQ' a f ac t° r x o appears in the mass-dependent logarithmic term in the infrared operator 
It is necessary in order to capture all poles parts proportional to 5(1 — x) of the 
integrated flavour violating antenna, A®(2 q ;3 g , 1q). 

In summary, in this subsection 5.3 we have shown that all integrated massive initial-final 
antennae have their pole parts related either to massive infrared operators or to splitting 
kernels or both. 

6. Check of .4° nfS 

A strong check can be performed on the integrated quark-antiquark antenna -^.qq by 
comparing its expression given in eq.(5.31) and known results from the literature on the 
leading order heavy-quark coefficient functions. 

To compare our results with 7 induced deep inelastic scattering we consider the con- 
traction of the hadronic tensor W^ u with the metric tensor —g^v This corresponds to the 
trace of the hadronic tensor, which in terms of the structure functions T2 and Fl is given 

by, 

-Wji = ~ d -^?L (z, Q 2 ) + ^ F 2 (z, Q 2 ) , (6.1) 



where the structure functions can be expanded in powers of the strong coupling constant. 
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To zeroth order in a s these structure functions are given by the simple parton model 
result 



(0) 



L,q 



(0) 



L,g 



0, 



T. 



(o) 



2,9 



6(l-z) , 



F {0) - 



(6.2) 



We find that the correct normalisation of Wn to be checked against the antenna A°. n ^. is 
as in the massless case [58] given by 

2 



W» = To - - - Ft 

d-2 " 2 d-2^ L - 



such that the following relation should hold 



A - x 



[ Q \l-e)-2m 2 Q ] 



1 - e 



2 - 2e 



(6.3) 



(6.4) 



with and being the leading order heavy quark coefficient functions. 

The factor on the left-hand side of this equation which multiplies the integrated antenna 
corresponds to the massive two-parton antenna Aq-q serving to normalise the three parton 
antenna and given in eq.(3.12). 

The coefficient functions are given up to finite order in e in [85] with a different choice 

o 2 

of variable x. x = = is used instead of our expression for x given in eq.(5.26) by 

Q 2 +2m 2 

x = — 2 , which includes the masses itiq of the final states Q and Q. 

We used this second definition of x in our expression of the integrated antenna A g .QQ 
since this definition of x is required in order to guarantee phase space factorisation in our 
subtraction formalism. Adapting the results given in [85] to this mass-dependent definition 
of x, the heavy quark coefficient functions read, 



[Q 4 (l - 8x + 8x 2 ) - 4Q 2 m 2 Q (-l + 3x + x 2 ) 



{Q 2 + 2m 2 Q f y 

+m A Q {A + 8x)] + [Q 4 (l - 2x + 2x 2 ) + AQ 2 m 2 Q {\ - 3x 2 ) (6.5) 



+m 4 (4 + 8x - 8x 2 )] log I - — - J \+0{e 



1 - v s 



^ (X ' Q2) = w^f\ v [q2{1 ~ x) + 2m Q ] + 2m ^ log (tt£) } + ° {e) - (6 - 6) 



(Q 2 + 2m 2 Q 
To the finite order in e, one can show that 



_Itt(i) , ^ (1) 



[Q\l -2x + 2x 2 ) + ±Q 2 m 2 Q (l - x 2 ) + m£(4 + 8x)] (6.7) 

I -v" 



" (Q 2 + 2m 2 Q f 

+ [Q 4 (l - 2x + 2x 2 ) + AQ 2 m 2 Q + mj(4 + 8x - 8x 2 )] log Q 
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with the variable v defined as in eq. (5.34) since the redefinition of x in [85] does not 
affect the definition of v. The expression given in eq.(6.7) coincides with the 0(e ) term of 
our expanded form for A g .qQ after we undo the normalisation of this antenna function by 
multiplying it by [Q 2 — 2i71q]. The relation given by eq.(6.4) is therefore herewith proven 
at O(e ) giving us a strong check on the integrated antenna A g .qQ itself. 

7. Application to top quark pair production at LHC 

In this section we shall give the colour-ordered real emission contributions for all partonic 
processes contributing to ti and tt + jet production at the LHC present at NLO. Together 
with these, we shall present their corresponding antenna subtraction terms which capture 
all single unresolved (soft, collinear and quasi-collinear) radiation of the real matrix-element 
squared for each partonic process involved. 

The results presented here for the real contributions and subtraction terms for the pro- 
cess pp —> ti +1 jet at NLO are essential ingredients for the computation of the double real 
contributions and their subtraction terms to the production of ti at NNLO. Our subtrac- 
tion terms enable to capture all single unresolved radiation present in those double real 
contributions. Furthermore, concerning the colour decomposition of the real matrix ele- 
ment squared for ti +1 jet production at NLO, in the limit where the heavy quarks present 
in the final state are taken massless, the colour decomposition provided here corresponds to 
the colour decomposition of the real matrix elements to the processes involving a massless 
quark-antiquark pair for two jet production at NNLO. 

For all partonic processes involved, starting from well-known amplitudes given in [19, 
22], the colour decomposition of the real matrix elements squared is presented and where 
possible, checked against results in the literature. Wherever possible, decoupling identities 
are used to reduce the size of the colour ordered decomposition of the real matrix-element 
squared and in some cases to eliminate the interference terms in those. 

The colour ordered subtraction terms are explicitly constructed as sums of terms in- 
volving the product of antenna functions with reduced matrix element squared and jet 
functions, as explained in Section 2. 

Concerning the notation of the matrix elements appearing in the real contributions, 
those matrix elements denoted by M represent colour-ordered amplitudes in which the 
coupling constants and colour factors have been omitted. Furthermore, to explicitly visu- 
alize the colour connection between particles in these colour ordered amplitudes, a double 
semicolon is used in the labeling of the partons contributing to a given matrix elements. 
This double semi colon serves to separate strings of colour connected partons. Partons 
within a pair of double semicolons are colour-connected. 

Since in the antenna framework a parton can only be unresolved with respect to its 
colour-connected neighbours, this notation helps to identify the unresolved limits present 
in a given colour ordered amplitude and therefore helps to construct the corresponding 
subtraction terms. Notationwise, we also denote gluons which are photon-like and only 
couple to quark lines, with the index 7 instead of g, to manifestly separate leading from 
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subleading contributions. In amplitudes where all gluons are photon-like no semicolons are 
used, since the concept of colour connection in not meaningful in those configurations. 

Concerning the notation in the subtraction terms themselves, the reduced matrix- 
element squared present in those are also to be taken without coupling constants and 
coupling factors. In those matrix elements, the crossed momenta are denoted with a hat, 
the remapped final-state momenta are denoted with tildes and the remapped momenta of 
initial state hard radiators are denoted by a bar and a hat 1 . 

7.1 ti production at LHC 

Following the general factorisation formula given in eq.(l.l) for hadronic collision pro- 
cesses, the real NLO contributions to the production of a massive quark-antiquark pair in 
a hadronic collision can be written as 



da 



R 



/ff{?[ 



/?(£i)/?(6)d<T ? ^QQ g + fq(ti)f g (b)da qg ^QQ q (7.1) 



+fq(€l)fg(&)d0qg^QQq 



+ fg(Zl)fg(&)<tegg->QQg \- 



fi(£j) denotes the parton distribution function of parton i, which carries a fraction £j 
of one of the incoming hadron momenta and da are the partonic cross sections. These 
process-dependent partonic cross sections da are, in turn given by 

da qq -^ Q Q g = d<S> 3 (k Q , k Q , k g ;p q , Pq )\M^ QQg \ 2 4 3) {k Q , k Q , k g ), (7.2) 
d& gg^QQq = d ^3(kQ,k Q ,k q ;p q ,pg)\M^ QQq \ 2 4 3 \kQ,k Q ,k q ), (7.3) 
d& gg^QQg = d ^( k Q^ k Q^ k ^Pg^Pg2)\M^ QQg \ 2 J ( 2 3 \k Q ,kQ,k g ), (7.4) 

where the momenta labels in the matrix elements are omitted for the sake of conciseness. 
Each of these process-dependent partonic cross sections involves a massive two-to-three 
parton phase space, d& 3 , the corresponding matrix element squared |M°| 2 and the appro- 

(3) 

priate jet function J 2 related to the selection criteria of 2-jet events. Out of three partons, 
from which two are a Q and Q , an event with two jets having each a heavy quark Q or 
a heavy anti-quark Q in them is formed. In the corresponding subtraction terms, the jet 

(2) 

functions will all be of the type J 2 and will only depend on the hard final-state momenta 
appearing in the reduced matrix-element squared. 

In order to obtain the subtraction terms, the colour decomposition of the real matrix 
elements present in the partonic cross sections given above in eqs. (7.2, 7.3, 7.4) has to be 
performed. However, since colour ordering does not distinguish between initial and final 
state partons, it is sufficient to consider the colour decomposition for the two unphysical 
processes — > QQqqg and — > QQggg, and obtain the matrix elements needed for ti 
production from two initial state partons by appropriate crossings. We shall follow this 
strategy in the following. 



1 This notation for the remapped initial state momenta was already used in [60] 
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7.1.1 Processes derived from — > QQqqg 

At amplitude level, the colour decomposition for this unphysical process — > QQqqg is 



M 5 (l Q ,2 ,3„4 9 -,5 fl ) =g 3 V2 {T a % lM 5 i2 , iz Ml{\ Q ,h g ,4 q --2 Q ,3 q ) 

+(T a ^ t2 5 nM M° 5 (l Q ,4 q ; ; 2q, 5 g , 3 q ) (7.5) 
-^{T^) ilM 8 i3M Ml^Q, 5 g , 2 Q - ; 3 q , 4 q ) 

— jj-(T a6 )i 3ji4 5i 1 j 2 M5(lQ, 2q; ; 3 q , 5 g , 4 q ) 

Squaring, and using the photon decoupling identities stating that 

M° 5 (l Q ,2 Q ,3 q , 4 q , 5 7 ) = M° 5 (1q, 5 g , 4 q ; ; 2 Q , 3 q ) + Mg(l , 4 q ; ; 2q, 5 9 , 3,) 

= M° 5 (1q, 5 g , 2 Q ; ; 3„ 4 ? -) + M°(l , 2^; ; 3„ 5 ff , 4 ? ) (7.6) 



gives 



|M 5 (1q,2q,3 ? ,4,-,5 5 )| 2 = g\N 2 c - 1) 

iV c (|Mg(l , 5 fl ,4 ? -; ; 2<g, 3 9 )| 2 + |Mjj(l Q) 4 5 ; ; 2 Q , 5 9 , 3 g )| 2 ) 

+ ^ 5 fl , 2 ; ; 3,, 4,-)| 2 + |M°(l g , 2 Q ; ; 3„ 5 g , 4 q )\ 2 (7.7) 



-2|^(1q,2q,3 <? ,4,-,5 7 )| 2 ) 



where in M®(1q, 2q, 3 q , 4 q , 5 7 ) the parton with momentum label 5 is a [7(1) photon-like 
gluon that only couples to the quark lines. 

Upon the crossing of 3 q and 4 q to the initial state, an expression for the squared matrix 
element \M^_^^Q g \ 2 in terms of colour-ordered amplitudes is obtained by replacing 3 q and 
4 q by 3 q and 4 q in the expression given in eq.(7.7). 

This, together with the phase space and the jet function as given in eq.(7.2) leads to 
the corresponding real emission differential cross section due to this partonic process in a 
colour ordered way. 

Analogously, the crossing of 4 q and 5 g in eq.(7.7) together with the phase space and the 
jet function in eq.(7.3) gives the colour ordered real emission corrections to the production 
of a top-antitop pair from the partonic process qg — > QQq. The corresponding subtraction 
terms for these two crossings are given below. 

After the crossing of 3 q and 4 q in eq.(7.7), we find that the subtraction term for the 
partonic process qq — > QQg is 



&O qq ^QQ g = 9 & { N l ~ l)d$ 3 (^lQ, k 2 Q, k 5g ] p 4q , p 3 q) 



xlN r 



AU4 q ;l Q ,5 g )\Ml((15) Q ,2 Q ,3 q r4 q )\ 2 4 2 \k r5 ,k 2 ) 
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+^(3 g -;2Q,5,)|^(lQ,(25)Q,3 g -,4 g )| 2 J 2 (2) (/c 1 ,fc 5 g) 



N r . 



AUl Q ,^,2 Q )\Mi((15) Q ,(25) Q ,3 q A q )\ 2 4 2) (k r5 ,k T5 



(7.8) 



+A»(4 q , 3 q ; 5 g )\Mi(l Q ,2 Q ,3 q , i q )\ 2 jf\k u k 2 ) 



Only A- Type quark-antiquark antennae are required in all three configurations. In the 
intial-final configuration, both flavour conserving and flavour violating massive antennae 
are needed. 

If and 5 g are crossed instead, the subtraction term for qg — > QQq will be constructed 
with massive flavour-conserving quark-antiquark A- Type and quark-gluon E-Type anten- 
nae in both final-final and initial-final configurations in addition to massless initial-initial 
A-type antennae. It reads, 



d& qg ^QQq = 9 ( N c ~ l)d<£>3(kiQ,k 2 - Q ,k 3q ;p 4q ,p 5g ) 



Ai(5 g] 1 Q , 2 Q ) [ \Mi((12) Q , 5 Q ,3 q , i q )f + \M^Q, (12) Q , 3 q , 4 q )\ 2 jffe k 3 ) 



+A° 3 (4 q , 5 g ; 3 q )\Mi(l Q , 2q, t q , 5 q )\ 2 J 2 (2 } (ki, h) 
1 



+^3 (4,; 3„ 1q) ( \M U A ((13) Q , 5 g , l g , 2«g)| 2 + \Ml((13) Q ,4 g , 5 g ,2 Q )\ 2 ) jf\k Tv k 2 



+-E»(4 q ; 3 q , 2 Q ) [\Ml(l Q ,5 g ,i g , (23)q)| 2 + \M U 4 (1 Q A 9 , 5 fl , (23)q)| 2 J jffa, k~) 
iUlfe 1 Q , 2 Q ) (\M° 4 ((12) q , 5q, 3„ 4,)| 2 + |A4°(Bq, (1~2)q, 3„ l q )\ 2 ) J 2 (2) fe k 3 ) 



N r . 



+A° 3 (4 q , 5 g ; 3 q )\Ml(l Q , 2 Q , l q , 5 q )\ 2 J {2) (hM) 
+±£°(4,; 3 q , 1q)|M2((1~3)q, t 7 , 5 7 , 2q)| 2 J 2 (2 > (k%, k 2 ) 



(7.9) 



+^ 3 u (4 g ; 3 q , 2q)\MI(1q, 4 7 , 5 7 , (23) )| 2 J 2 (2 > (fci, fc 



2:-! ' 



7.1.2 Partonic process gg — > QQg 

The colour decomposition for the unphysical process — > QQggg is 

M 5 °(1q,2 ,3„4„5,) = (gV2f £ [T a ^T a ^Ml{l Q , i g , j g , k g , 2 Q ). (7.10) 

(ij',fc)6P(3,4,5) 

Squaring and crossing gluons 4 S and 5 3 to the initial state gives 
|M 5 °(l Q ,2Q,3 s ,4 s ,5 s )| 2 =/(iV 2 -l) 



x 



E 

(i,i)eP(4,5) L 



Ni (\MUlQ,3 g ,i g J g ,2 Q )\ 2 + \MUlQ,i 9 ,3 9 ,j 9 ,2 Q )\ 



+\M° 5 (l Q ;i g ,j 9 ,3 g ,2 Q )\ 2 ) 
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-(|^0(lQ,3„^,i 7 ,2Q)| 2 + |^(l Q ,? s ,3 3 ,i 7 ,2Q) 
+ \M° 5 (l Q ,i g ,j g ,3 y ,2 Q y 2 ' 



(7.11) 



+ 



N 2 + l 



|^(l Q ,3 fl ,4 fl ,5 fl ,2 )p 



We have used the following photon decoupling identity, 

M 5 {l Q ,ig,jg, kg, 2 Q ) + M 5 (l Q ,ig, kg,jg, 2q) + M° 5 (1 Q , kg,ig,jg, 2 Q ) 

where gluon k is a U(l) boson decoupled from the other gluons, and 

M° 5 {lQ,3 g ,4g,5g,2 Q ) = M° 5 (l Q ,ig,jg,kg,2 Q ), 

(i,j,fc)eP(3,4,5) 



(7.12) 



(7.13) 



where all gluons are photon-like. 

The colour decomposition of the matrix element squared with two gluons in the initial 
state given by eq.(7.11) is obtained as follows: After squaring the amplitude for the process 
— > QQggg and expanding the sum of permutations concerning the three final state gluons, 
the crossing of two gluons 4 g and 5 g is performed. As the final state gluon labelled 3 g is 
fixed, the terms are then regrouped and a sum over the permutations for the two initial 
state gluons 4 9 and 5 g only is performed. 

Only colour ordered squared matrix elements are involved in eq.(7.11) for which unre- 
solved radiation can be captured by a single antenna function. Massive flavour-conserving 
quark-antiquark A- Type antennae, massive quark-gluon D-type antennae together with 
massless initial-initial gluon-gluon F-Type antennae are required. 

The subtraction term reads, 

d& gg^QQg = ^ ( N c ~ l)d$ 3 (^lQ, k 2 Q, k 3g ;p 4 g,P5g) 

^3(^;lQ,2Q)(|Al2((r2)Q,3 p ,i 5 ,fQ)| 2 + |^(?Q,3 3 ,i s ,(l~2)Q)| 2 

MMl((l2) Q rjg,3gt Q )\ 2 + \Ml(i Q ,jg,3g(l2) Q ^ 

+Z?3^^;3 9 ,lQ)|A4°((r3)Q,i 9 ,j 9 ,2Q)| 2 jf ) (fc T 3,A; 2 ) 
+Z?3^i,;3 9 ,2Q)|A40(lQ,i s ,I 9 ,(2~3)Q)|V 2 (2) (/c 1 ,fe23) 

+F°(i g Jg-,3g)\Ml(l Q X,lAQ)\ 2 4 2) (h,k~ 2 ) 



{ £ 


N 2 


L (i,j)eP(4,5) 





-AUlQ,3g,2Q)\M^(l3) Q , i gJg,(23)Q)\ 2 4 2) (k r3 ,k 2li ) 

-^ 9 ;1q,2q)(|A4°((^ 

MMl((12) Q ,j g ,3 g t Q )\ 2 + \Ml(iQ,j g ,3g(12) Q )\ 2 



(7.14) 
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-D)i(i 9 -,3 g A Q )\Mlm)Q,hJy^Q)\ 2 4 2 \k^,k 2 ) 
-Z?3^^;3 s ,2 )|^(l Q ,j- 7 ,I 7 ,(2%)| 2 J 2 (2) (fc 1 ,A;23) 

+^Al(i g] l Q ,2 Q ){\Ml((12) Q ,3^rh,tQ)\ 2 
MMl(lQ,3^j^(l2) Q )\ 2 )4 2 \k r2 ,k 3 ) 



+ 



Ni 



^ j ^(lg, 3 g , 2 Q )\Mi((13) Q , 4 7 , 5 7 , (23) )| 2 J< 2) fcg) 



7.1.3 Consistency check 

We have performed a powerful check on all subtraction terms required for ti production 
given above. For each of those subtraction terms, which is a sum of terms multiplied by 
colour factors proportional to N c , we have verified that it gives the correct non-colour 
ordered collinear or/and quasi-collinear behaviour. We have checked this feature for all 
collinear and quasi-collinear limits present in all subtraction terms presented above. We 
have verified that each subtraction term obey 



a\\b 



9 2 C 



Pab(z) 
Sab 



x \M° m \ 2 x d$ m , 



(7.15) 



where C = Ca, Cf, Tr is the corresponding Casimir, is the non- colour ordered reduced 
matrix element and P a b(z) stands for a massless or massive spitting function governing the 
particular collinear or quasi-collinear limit as defined in Section 4. 

Let us now check all collinear and quasi-collinear limits of the subtraction term ^ a gg ^.QQ g 
given in eq.(7.14) in this way. This subtraction term has the following collinear and quasi- 
collinear limits: 1q||3 s ,3 9 ||4 9 and 1q||4 s . 

The following relation between the non-colour ordered and colour ordered amplitudes 
squared for the process gg — > QQ is needed. It reads, 

/ N 2 

\M° 4 (i Q ,2 Q ; lg ,j g )\ 2 ' ■ 



Nr. 



N^[\M u 4 (,l Q ,i g ,j g ,2 Q ) 



+ \Ml(l Q ,j g ,i g ,2 Q )\ 2 ) -\M 4 (l Q ,l ( ,j~n 2 Q)f 



(7.16) 



When the final-state quasi-collinear limit 1q||3 9 is taken in the subtraction term given 
in eq.(7.14), we obtain 



da 



1qI|3 9 



X 



gg^QQg 

2 



N 



N 



'Ml > 



S13 
2 , I aaO, 



Nt \MI{{1 + 3) , 4 9 , 5 g , 2q)\* + \Mt((l + 3) Q , 5 g , 4 g , 2 Q 



|^((1 + 3) q ,4 7 ,5 7 ,2q)|^ 



(2) - _ 

^2{k(i + z)Q,k 2 Q;p Ag ,p bg )r 2 J (k 1+3 ,k 2 ) (7.17) 



= g 2 C F 



Pqg^Q(z,fJ-, 



'Ml ' 



S13 



|M 4 ((l + 3) Q ,2 (3 ,4 9 ,5 s )| i 
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xd$ 2 (% +3 )Q, k 2 Q'>P4g>P5g) J 2 2 \ k i+3, h), 

where we have used eq.(7.16) with the momentum of the massive final state quark 1q 
being given by (ki + ks). In the limit 1q||3 s , the contributing terms in eq.(7.14) are 
proportional to the antenna A^(lQ,3 g ,2Q) multiplied by a reduced colour-ordered matrix 
elements squared involving the remapped momenta k,f 3 \ and fe( 2 ~3) ^ n ^ ne ■'■qII^s limit, those 
remapped momenta k/f 3 \ and fc( 2 ~3) are respectively given by (k\ + ks) and k 2 . 

Similarly, in the same subtraction term given in eq.(7.14) we can see that in the 3 9 ||4 9 
limit we obtain 

d& s 3 9 \\\ 3 ^ - I P gg ^ G {z) 



gg^QQg y jy2 



c 



S34 



N' c \M U 4 (1 Q , (4 - 3) g , 5 g , 2 Q )\* + \Mt(l Q ,5 g , (4 - 3) g , 2 Q ) 



d$ 2 {h Q , k 2 Q ; P(4_ 3 ) fl , P5g) jf ] {h , k 2 ) (7.18) 



-|A^(1q,(4-3) 7 ,5 7 ,2q) 

= g 2 C A ^ \M2(lQ, 2 Q , (4^3) g , 5 g )\ 2 

S34 

(2) 

xd$ 2 (felQ, fc 2 Q;P( 4 -3) 3 ,P5 ff )^2 (fcl, fo)- 

We have used eq.(7.16) with incoming momenta (p4 — ks) and p^. In the 3 g \\4 g limit, the 
contributing terms in eq.(7.14) are proportional either to D®(4 g ; 3 g , 1q) or to F3 (4 S , 5 9 ; 3 g ). 
In the first case, the remapped momenta (13) and 4 present in the reduced matrix elements 
multiplying by D®(4 g ; 3 g , 1q) are given by k\ and (p4 — ks) respectively. In the second case, 
the remapped momenta 4 and 5 present in the reduced matrix-element squared multiplying 
F®(4 g ,5 g ;3 g ), are given by (p± — ks) and ps respectively. 
Finally, for the 1q||4 s limit, we see that 

6 (N 2 - 1\ 2 P q ^ G {z,l4g) 

d °gg^QQg ~ ► 9 [—^ ) — 

N 2 ( \MK(4^1)q, 5 g , 3 g , 2 Q )\ 2 + |M2((4^1)q, 3 g , 5 g , 2 Q ) 

(2)/ 



-|^((4-1)q,3 7 ,5 7 ,2 

Pqq^G{z, fJ%g) , . , / 



d$2(fc 2 Q,fe3 S ;P(4-l)Q>P53)^2 (*<2>*&) (7.19) 



= g , Cp ^ -^ |M 4 "((4 - 1) , 2 , 3 9 , 5 9 )r 

X d$ 2 (fe 2 Q , &3 9 5 P(4-l)<5 > Phg)^ (*2 , *&)■ 

Again, the relation between the full matrix element squared and the partial amplitudes 
given by eq.(7.16) has been used, in this case the final state quark denoted by 1q and the 
gluon denoted by 4 g have been crossed to the initial state. 

Since colour decomposition does not distinguish between initial and final state coloured 
particles, this crossings can be safely done with the relation still being true. The terms 
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in eq.(7.14) that contribute to this limit are those involving ^3(4^; 1q, 2q) and the cor- 
responding remapped momenta appearing in the reduced squared amplitude multiplying 
this antenna become k^ ~^ ^2 and p| — > (j>4 — k±) . 

These same verifications have been performed on all collinear and quasi-collinear limits 
of all subtraction terms listed in this section providing us with strong consistency check on 
our results for the subtraction terms required to compute the cross section for pp — > ti at 
NLO. 

7.2 ti + jet production at LHC 

The real emission correction to the production of a massive quark-antiquark pair in asso- 
ciation with a jet in a hadronic collision is given by, 

+fq(t,l)fg(&) d Vq g ^QQq g + fq(£,l)fq(&)da q g^ Q Q qq - 

+mi)m2)da qq _ QQqq + mi)m2)da qq _ QQqq ] (7.20) 

+ lU(&)fq(&)da qq ^ Q Q qlq , + fq(£i)fq>(&)da qql ^QQ qql 

q+q' 

+ fq{t,q)fq l {&)dd qq ,^ Q Q qql + fq(£l)fq> (&)da qq ,^QQ qq ,] 
+/< 7 (£l)/< 7 (6) (dVgg->QQgg+d& gg ^ Q Q qq ) j, 

with the partonic cross sections given by 

dtqq^QQgg = d$ 4 (k Q , k Q , k gi , k g2 ; p q , p q ) \ M^ QQg g \ 2 jf } (k Q , k Q , k gi ,k g2 ) (7.21) 
d& qg^QQqg = d MkQ,kQ,k g ,k q ;p q ,p g )\M^ QQq f4 4 \kQ,kQ,k q ,k g ) (7.22) 

d°qq^QQq>q> = d<5> 4 (k Q , k Q , kg, , kg, \ Pq,Pq) \ M^ QQg ,., | 2 J< 4) (k Q , k Q , , (7.23) 

d^^QQ^, = d^Q,*^,^, k q r )Pq ,p q ,)\M Q qq ,^ QQqq ,\ 2 {k Q ,kQ,k q ,k q ,) (7.24) 
d ^99->QQ99 = ^(kg,^,^!,*^;^-^)!^ (7.25) 
d& gg^QQqq = d ®4( k Q> k Q> fe <?> k v Pgn Pg2)\ M gg ^ Q Q q q\ 2 jt\ k Q> k Q^ h)- ( 7 - 26 ) 

The jet functions appearing in the partonic cross section are all of the type and 
correspond to the selection criteria of 3-jet events. Out of four partons, from which two are 
a QQ pair, an event with three jets is build. From these 3-jets, one jet is made of massless 

partons and two other jets have each a heavy quark Q or a heavy antiquark Q in them. 

(3) 

In the corresponding subtraction terms, the jet functions are of the type Jg and those 
depend only on the hard final-state momenta appearing in the reduced matrix-element 
squared. 

In order to obtain the colour decomposition of the matrix element squared for the 
partonic processes defined above and given in eq.(7.26), we will follow the same strategy 
as for ti production. We will use the colour decomposition of unphysical processes and 
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then consider appropriate crossings. For ti +1 jet production the unphysical processes to 
be considered are — > QQqqq'q', — > QQqqgg and — > QQgggg- 

However, the presence of one additional parton in the final state introduces a few dif- 
ficulties. In the first place, the number of partial amplitudes as well as the number of 
unresolved limits to subtract increases. Also, identical quark flavour contributions must be 
considered. But most importantly, in the contributions related to the partonic processes 
derived from — > QQqqgg and — > QQgggg, the colour decomposition of the partonic 
amplitudes squared leads to interferences between partial amplitudes with different colour 
orderings that cannot be removed using any decoupling identities. Those require subtrac- 
tion in an uncommon way in the antenna formalism, which we shall explain below. 

7.2.1 Interference terms 

The subtraction of infrared singularities in a colour ordered squared amplitude is easy 
within the antenna formalism: a suitable antenna function multiplied by a reduced squared 
matrix element with its momenta properly remapped accomplishes the task. However, this 
is not the case for interferences between partial amplitudes with different coulour orderings. 
Those interference terms which are most generally of the form 

M° n+1 (...,a,s,b,...)M° n+1 (...,c,s,d,..y 

with gluon s colour connected to partons a and b in one amplitude and colour-connected 
to partons c and d in the other amplitude, lead to soft singularities when integrated over 
the phase space. Those singular behaviours cannot be straightforwardly subtracted with 
just one suitable antenna function multiplied by a reduced squared matrix element with 
remapped momenta. 

In order to understand how to subtract the soft singularities in these interferences 
terms, we appeal to the factorisation properties of colour ordered matrix elements at the 
amplitude level. Quite generally, when a gluon s with helicity A becomes soft between 
partons a and b, the colour ordered amplitude factorises [87, 88] as 

M° n+1 (..., a, s x ,b, ...) e x ( Ps ) ■ [J a ( Ps ) - J b ( Ps )} M° n (..., a, b, ...), (7.27) 

where e\(p s ) is the polarisation vector of the soft gluon, and Ja(p s ) is the soft gluon current, 
given by 

J£(p.) = / a . (7.28) 

V2p a -Ps 

Eq. (7.27) holds in d dimensions and is absolutely general: it does not depend on the 
identity of partons a and b (they can be either gluons, massive quarks or massless quarks) 
nor on their helicities. Therefore, summing over the helicities of the soft gluon, we find 
that the interference of two partial amplitudes with different colour orderings factorises in 
the soft limit, k s — > 0, as 

Yl • M n+l(-' fl > S A, b, ...)M° +1 (-.., C, 8\, d, ...)t 

A=± 
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k I^> ^(PsV-MJ [JaiPs) ~ J b (p S )Wc(p S ) ~ Jd(Ps)]u (7.29) 

xM n (-,a,b,...)M n (-,c, d,...) f . 
Since colour ordered amplitudes are gauge invariant, we can replace 

£ e£(p.)£ A (p.) -> -«T (7.30) 

A=± 

and thus obtain the following limit for the interference term: 

M° n+1 (...,a,s,b,...)M° n+1 (...,c,s,d,..y (7.31) 



k s -*0 I Sad She S ac Sbd 



^jM n (...,a,b,...)M° n (...,c,d,..y. 



^SasSds Sb s S cs S as S cs Sbs^ds 

In this soft limit, we see that the interference term given by 

Mn+i(-,a, s, b, ...)Mn+i (-, c, s, d, ...) t 

factorises into a difference of four (half) eikonal factors times the interference of the reduced 
(i.e where the soft gluon is absent) interference term given by 

M° n (...,a,b,...)M° n (...,c,d,..y. 

Each of these half eikonal factors can be reproduced by a single antenna function. The 
soft singularities of an interference term like .M° +1 (..., a, s, b, ...)7W^ +1 (..., c, s, d, for the 
case where a, b, c and d are all different can therefore be subtracted with 

^X 3 (a, s, d)M° njl (..., as, bs, cs, ds, ...) f 

+\x°(b, s, c)A< 2 (.., as, bs, ...K, 2 (.., cs, ds, ...)t (7.32) 

--X$(a, s, c)M° i3 (.., as, 6s, ...)7W° 3 (..., 6s, ds, ...) t 

s, d)M° nA {..., as, bs, ...)M° >4 (..., c~s, da, ...)*. 

In this equation, each of the antenna functions denoted as X® can be either final-final, 
initial-final or initial-initial, depending on whether the partons a, b, c and d are in the 
initial or final state, and they can be either massive or massless depending on whether the 
hard radiators are massive or massless. The absence of mass terms in eq.(7.31) is respected 
by eq.(7.32) even when one of the hard radiators (and therefore all the antenna functions 
involving that parton) is massive: In the soft limit, although a massive antenna yields a 
massive soft eikonal factor including explicitly mass terms, as seen in Section 4, in the 
combination of massive antennae required in eq.(7.32) those mass terms cancel amongst 
each other. 

In eq.(7.32), the different labels on the reduced matrix elements denoted by \ do not 
mean that the matrix elements themselves are different, but that the momentum mapping 
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is, in principle, different for each term. The remapping in each case is done in accordance 
with the type of antenna function involved. 

A very important feature about this way of treating the interference terms is that 
not only soft singularities are subtracted with eq.(7.32), as it was originally intended, but 
also no "extra" collinear singularities are introduced by this subtraction term itself. The 
collinear limits of gluon s collinear with partons a, b, c and d are correctly dealt with. 
Indeed, if a and b are different from c and d, the interference term does not contain any 
collinear singularities. It has square-root singularities that do not give rise to e poles upon 
integration. In this case, the subtraction terms in eq.(7.32) do not contain any collinear 
singularities either: all the collinear limits introduced by the first two terms, are subtracted 
by the last two. And, in the case where two partons are equal in eq.(7.29), for example, 
when a and c are the same, the interference term develops a collinear singularity when a\\s. 
In the limit where s becomes soft, setting c = a in eq.(7.29) gives 

M° n+1 (..., a, s, b, ...)M° n+1 (..., a, s, d, ...)t (7.33) 

fes-H3 / S a d S a b S^d m V\ » A Q i , U/l 0/ j \t 

— > + 2~ M n {...,a,b,...)M n {...,c,d, ...)'. 

\ s asSds SasSbs S^gSds s as / 

The corresponding subtraction terms are 

l>X$(a, s, d)M n>l (..., as, bs, ...yw° ;1 (..., as, ds, ...) f 

+^Xl(a, s, 6)M°, 2 (.., as, bs, ...)A< 2 (..., as, ds, ...)* (7.34) 

-\xl{b, s, d)M° n ^(..., as, bs, ...)M° n ^(..., as, ds, ...) f , 

from where it can be seen that these combination of subtraction terms also subtracts the 
a\\s (quasi)-collinear limit, in addition to the soft s limit. 

Taking this construction of the subtraction terms required for the interference terms of 
massive amplitudes into account, let us now present our results for the real contributions 
and their corresponding subtraction terms for all partonic processes involved in tt+jet 
production at NLO. We start with processes involving only quarks and derived from — > 
QQqqq'q', which is the easiest case since only colour-ordered matrix-element squared are 
involved. 

7.2.2 Processes derived from — > QQqqq'q' 

We choose to separate the colour decompositions of the real contributions related to pro- 
cesses derived from — > QQqqq'q with and without identical quarks explicitly. 

The colour decomposition for the unphysical partonic process — > QQqqq'q' hi the 
non-indentical flavour case (q / q') is given by, 

M &{^Q,2Q,3 q Aq,^ q ' ,Qq') = g 4 [5i ltii 5i 3t i 6 5i 5t i 2 Ml(lQ,4q; ;3 q ,6g' ;5g/,2g) 

+^il,i 6 ^3,i 2 ^5,*4-^6(lQi j 3g, 2q; ; h q >, 4g) 
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-^iuie^iABM-^eO-Q^ 6 g r ' ' 3 <?' 4 ?; ; 2 q) ( 7 - 35 ) 

~~ jY~^i,i2^3,j6^5,M-^6(lQ; 2q; ; 3 g , 65'; ; 5 g /,4^) 

+^2 ^i,j 2 ^3,i4^5,«6- M 6( 1 Q' 2 o; ; 3 <?> ; V' 6 «0 • 

Squaring, and using the fact that the partial amplitudes satisfy 

M° 6 {l Q ,2 Q ;3 q ,4 q ;5 q ,,6 q ,) 
= M° 6 (l Q Ag, 5 3„ 6 r ; ; 5,/, 2$) + M° 6 (1q, 6 r ; ; 3 5 , 2^; ; 5^-) (7.36) 

= i (K(lQ>%; ; 3„ 2q; 5,/, 6 r ) + 6 r ; ; 3„4 ? -; 5 q ,, 2 Q ) + M° 6 (1 Q , 2 Q ; ; 3„ 6 ? ; 5^,4^)) 

gives 

|M 6 °(1q,2q,3 9 ,4 ? -,5^6 9 -0| 2 = /(iV c 2 - 1) 
|tV c ^°(1 q ,4 9 -; ; 3„ 6 r ; ; 5,/, 2 Q )\ 2 + |Mg(l Q ,6 r ; ; 3„ 2^; ; 5 q ,, 4 g -)| 2 ) 

+^(|K(lQ,4 $ ;;3 g ,2Q;;5 g ^6 g -0| 2 + |K(lQ,6 g - S ;3 g ,4 9 -;;5^,2Q)| 2 (7.37) 

2q; ; 3„6 r ; ; 5^,4 5 )| 2 - 3\M° 6 (1q, 2 q ; ; 3^,4,-; ; 5 q > , 6 q >)\^ } 

which represents the colour decomposition for squared amplitude of the unphysical process 
— > QQqqq'q' (with q and q' of different flavour). It has no interference terms. Depending 
on which partons are crossed in eq.(7.37) a sum over final state quark flavours might need 
to be performed. When a quark-antiquark pair of the same flavour is crossed to the initial 
state resulting in the the colour ordered squared amplitude for the process qq — > QQq'q', 
the result must be multiplied 2 by (Np — 1). 

The identical flavour case is obtained taking the difference of two colour-ordered non- 
identical cases. At the amplitude level, we use the following relation 

M 6 (l Q ,2Q,3 g ,4 g -,5 g ,6 g -) = M 6 °(lQ,2 o ,3 g ,4 5 ,5 g ,,6 5 - M°(1 Q ,2 Q , 3 q ,6 q ,5 qf ,4f ), (7.38) 

which, upon squaring, gives the colour decomposition of the matrix-element squared for 
the identical flavour case, 

|M 6 °(1 Q , 2 Q ,3 q ,4 q , 5„ 6 q )\ 2 = g 8 (N 2 - 1) 

x ^N c ^\Ml(l Q ,4q, 5 3„ 6 ? ; ; 5 q >, 2 Q )\ 2 + \M° 6 (l Q ,6 q ; ; 3„4^; ; 5 q >, 2 Q )\ 2 

+lK(lQ,6 9 -,;;3 9 ,2Q;;5 9 ^4 g -)| 2 + |K(lQ,4 g -,;;3 (? ,2Q;;5 (? ,,6 g -)| 2 ) 

+2Re(A^0(l Q , A ql - ; 3„ 2 Q ; ; 5,/, 6 9 -)M°(l Q , 4,-; ; 3,, 2 Q ; ; 5,/, 6 r )t) 

2 Note that the result needs to be multiplied by (Nf — 1) and not Nf since we are explicitly separating 
the identical flavour contributions from the non-identical ones. 
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+2Re(K(l Q , 4 r ; ; 3,, 6,-; ; 5,/, 2 q )M%{1q, 4 g -; ; 3„ 6 r ; ; 5,,, 2g)t) 
+2Re(.M°(l , 6 q ; ; 3„ 2^; ; 5 q ,, 4 f )M°(lg, 6 r ; ; 3„ 2^; ; 5^,4 ? )t) 
+2Re(Mg(l g , 6 5 ; ; 3„ 4 r ; ; 5,/, 2q)M{|(1q, 6 r ; ; 3„ 4 ? -; ; 5,,, 2$)*) 
-2Re(Mg(l g , 2q; ; 3„ 6 5 ; ; 5,/, V )Mg(l Q , 2 Q ; ; 3„ 4,-; ; 5,,, 6 r ) f ) 



(7.39) 



+ ^(^l^6(lQ,4 g -;;3 g ,2Q;;5 (? ,,6 g -0| 2 + |^(lQ,6 5 ;;3 g ,2Q;;5 (? ,,4 g -0| 2 

6 f ; ; 3 9 ,4 9 -; ; 5,,, 2 Q )\ 2 + \M° 6 (l Q A q >; ; 3„ 6 q ; ; 5 q ,, 2 Q )\ 2 
+\M° 6 (1q, 2 Q - ; 3„ 6 r ; ; 5^,4 ? )| 2 + |A<g(l , 2^; ; 3,,%; ; 5 q >, 6 q )\ 2 



-3|Mg(l , 2 Q ; ; 3„4 ? -; ; 5^,6^)| 2 " 3|Mg(lg, 2 Q ; ; 3„ 6 ? ; ; 5 q ,, 4 q ,)\ 2 



2 (fflte(-Mg(lQ, 2q; ; 3„ 6 ? -; ; 5 ? >, %)Mg(l , 2^; ; 3„ 4 ? ; ; 5,/, 6 ? -, ) T ) 
-2Re(Mg(l Q , 2 ; ; 3„ 6 q ; ; 5 9 >, 4^)Mg(l , 2^; ; 3„ 6 r ; ; 5^,4 ? )t) 
-2Re(A^°(l Q , 2 ; ; 3„ 4 ql ; ; 5,,, 6,-)Mg(l , 2 ; ; 3„ 4,-; ; 5,,, 6g,)V 



The interference terms present in eq.(7.39) lead only to finite contributions when inte- 
grated over the phase space. Indeed, they lead to square root singularities which do not 
need subtraction. 

Furthermore, apart from those interference terms, the identical-quark flavour contribu- 
tions involve all six colour-ordered matrix element squared appearing in the non-identical 
flavour case plus additional ones. In all cases, identical or non- identical quark and in all 
crossings required, the only singular behaviours which have to be captured by the anten- 
nae in the subtraction terms are final-final and initial-final collinear singularities between 
identical-flavour massless quark-(anti)quarks. Due to the presence of a massive radiator 
in the final state, only one antenna, the quark-gluon massive antenna E® is required in 
final-final and initial-final configurations. 

After the crossing of 3 q and 4 q to the initial state, in the non-identical matrix element 
squared for the process — > QQqqq'q' given in eq.(7.37) we obtain the colour ordered 
matrix element squared for qq — > QQq'q'. Its subtraction term is, 




d$4(£:iQ , k 2 Q , k 5q > , k &ql ; P3q,P4q) 



+\M° 5 ((15)q, (56) 9 , 4 9 ; ; 3 q , 2q)| 2 ) J? } (%, k 2 , k~) 



+E 3 (2 Q ,5 Q/ ,6 q ,) (\M° 5 (1q, 4 ? ; ; 3 q , (56) 9 , (25) )| 2 



+\M° 5 (1 Q , (56) 9 , 4 9 ; ; 3 q , (25) )| 2 ) jf (k u k Tb , k~) 



+ W C [^(iQ^VAO (\M° 5 ((15) Q ,2 Q ;;3 q , (56)„4,)| 2 
+\Ml((15) Q ,(56) g ,2 Q ;;3 q ,4 q )\ 2 



(7.40) 
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-2|M°((15) Q , 2q, 3 q , 4 q , (56) 7 )| 2 ) jf fc 2 , fc ; 

+£ 3 °(2q, 5,,, 6 r ) (|Mg(l Q , (25) ; ; 3 ? -, (5~6) g ,4 ? )| 2 
+|M°(l g ,(56) 9 ,(25) ;;3 ? -,4 ? )| 2 

-2|X0(l Q ,(25)Q,3 g -,4 g ,(56) 7 )| 2 ) jf (fci.fc^fc 



If, instead, 4g and 6^/ are crossed in eq.(7.37), we obtain the squared matrix element 
for the process qq' — > QQqq'. The corresponding subtraction term reads 



da 



qq'^-QQqq' 



5 8 (iV 2 - 1) 



d$4 (ki Q , k 2 Q , k 3q , k 5q > ; p Aq , Pe q ' ) 



E»(4 q ; 3 q , 1q) (|^((13)q, 6,,; ; 5 q >, 4 g , 2 Q )\* 
+\M° 5 ((13) Q , tg, V; ; 5^,2q)| 2 ) jf (%, fc 2 , fc 5 ) 



+£ 3 u (4 g ; 3 9 , 2g) |Mg(l , 6,,; ; 5 q >, A g , (23)q)|< 



4 a , 6,.; ; 5^,(23)q)|^J J< 3) fc 5 
+^(6^; 5,,, 1q) (|Mg((f5) , 4,; ; 3„ 6„ 2q)| 2 

+|.M°((1~5)q, 6 fl , 4,; ; 3„ 2q)| 2 ) jf >(fejg, k 2 , k 3 ) 
+E$(6 q ,;5 q ,,2 Q ) {\M° 5 (l Q ,4 q ;;3 q ,6 9 ,(25) Q )\ 2 

MMl{i Q XA q v,\A^)Q)\ 2 ) jf\ki,k^M) 



+ 



Nr. 



E»(4 q ; 3 q , 1 Q ) (\Mi((13) Q ,2 Q ; ; 5^, 4 fl , 6 9 0| 2 

+ |^0((r3)Q,l„2 ;;5^,6 g O| 2 

-2\Mi((13) Q , 2 Q , 5 q/ , V, I 7 )| 2 ) Jffe fc 2 , fc 5 ) 

+^3°(4 g ; 3„ 2q) (|^(1q, (2~3)q; ; 5,,, 4 fl , 6,0| 2 



+|M°(l ,4 fl ,(23) ;;5 9 ,,6 9 



-2|^(1q,(23)q,5^,6 ? ,,4 7 )|^ jf (k u k^,k 5 ) 
+E° 3 (6 q ,;5 q ,, 1 Q ) (|Mg((f5)g, 2^; ; 3„ 6„ 4 g )| 2 
+ |X0((r5)Q,6„2 ;;3„4 g )| 2 
-2|M°((1~5) Q , 2q, 3„ 4„ 6 7 )| 2 ) jf (fejg, fc 2 , fc 3 ) 



+j E 3 u (6^;5^,2q) ^|^(1q,(25)q;;3„6„4 9 )| 2 
+\M° 5 (l Q ,6 g ,(25) Q ;;3 q ,l q )\ 2 
+\M° 5 (1q, (2~5)q, 3„ 4,, 6 7 )| 2 ) jf >(*;!, k 3 ) 



(7.41) 
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For the identical flavour-case, the same two crossings, 3 q and 4 q on one hand and the 
crossing of 4 q and 6^ on the other hand need to be considered. As the real matrix-element 
in the identical quark case are given by those in the non-identical quark case plus additional 
terms, similarly the subtraction terms for the identical-quark case are obtained by adding to 
the subtraction terms valid for the non-identical case additional subtraction terms related 
to the new colour ordered matrix-element squared appearing only in the identical-flavour 
case only and given in eq.(7.39). 

If 3 q and 4 q are crossed in the identical flavour squared matrix element given in eq.(7.39), 
the squared matrix element for qq — > QQqq is obtained in a colour ordered way. The 
subtraction term for this process is, 

g 8 (N 2 — 1) 

^qq^QQqq = \ d $<t( k lQ, k 2 Q, k 5q , k Qq ' PZq , P4 q) 



X<N r 



E° 3 (1q, 5„ 6q) (|Mg((15)o, 4 ? ; ; 3 q , (56),, 2 Q )\ 2 

+\M° 5 ((15) Q , (56) s , 4 9 ; ; 3 q , 2 Q )\ 2 ) k 2 , k~) 

+E° 3 (2 Q , 5„ 6,-) (\M° 5 (1 Q , 4,; ; 3 q , (56),, (25) )| 2 

+ |X0(lQ,(56) s ,4 9 ;;3 g -,(25) )| 2 ) jf fcgg) 
+E° 3 (4 q ; 5„ 1q) (\M° 5 ((15)q, 6 q ; ; 3 5 , l g , 2 Q )\ 2 

+\M° 5 ((15) Q , tg, 6 q ; ; 3 q , 2q)| 2 ) jf (%, fc 2 , fc 6 ) 
+£ 3 °(4 f ,; 5„ 2q) (\M° 5 (1q, 6,-; ; 4 g , (25)q)| 2 

4 9 , 6 5 ; ; 3 5 , (25) )| 2 ) jf ^ fc 6 ) 
+E° 3 (3 q ; 6 q , 1 Q ) (\M° 5 ((16)q, 4 q ; ; 5„ 3 fl) 2q)| 2 



1 



+ |.M°((16) , 3 9 , 4 9 ; ; 5„ 2q)| 2 ) jf fc 2 , fc 5 ) 
+£°(3,-; 6,-, 2 ) 4,; ; 5„ 3 fl , (26)q)| 2 

+\M° 5 (1q, 3 g , \\ ; 5„ (26) )| 2 ) Jf (k u k T6 , k 5 ) 
E° 3 (l Q , 5„ 6 ? -) (|Mg((f5) , 2 ; ; 3 ? -, (5~6)„ 4,)| 2 
+|^°((r5)Q,(56) s ,2Q;;3 g -,4 g )| 2 
-2|M°((1~5) 0) 2q, 3 q , 4 q , (56) 7 )| 2 ) J< 3 >(fcjg, fc 2 , fcgg) 

+i?3 (2 , 5„ 6,-) (2~5)q; ; 3 q , (56),, 4,)| 2 

+ |^(lQ,(56) s ,(2'5) ;;3 g -,4 g )| 2 

-2|^°(l Q ,(2~5)Q,3 g -,4 9 ,(56) 7 )| 2 ) 4\k u k Tb ,k~) 
+Ei(4 q ; 5„ 1q) (|^((f5) Q , 2q; ; 1„ 6 5 )| 2 



(7.42) 
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+|A^((15) Q ,4„2Q;;3 g -,6 g -)| 2 
-2|Mg((f5) , 2,3,3,-, 6,-, I 7 )| 2 ) 4 3 \k r5 , k 2 , k 6 ) 
+E° 3 (A q ; 5 q , 2 Q ) (\M° 5 (lQ, (25)q5 5 %, t fl , 6 g -)| 2 
+ |^0(lQ,l„(25) ;;3 g -,6 g -)| 2 
-2|^°(l Q ,(25) ,3 ? -,6 ? -,l 7 )| 2 ) 4 3 \h,k T5 ,k 6 ) 

+E° 3 (3 q ; 6,-, l g ) (|^((f6) Q , 2q; ; 5„ § fl) 4 g )| 2 
+ |^((r6) Q ,§„2Q;;5 (/ ,4 g -)| 2 
-2|Mg((16) 0) 2q, 5„ 4 9 , 3 7 )| 2 ) jffe fc 2 , fe 5 ) 

+£ 3 D (3 ? -; 6,-, 2 ) (2~6)q; ; 5,, 3 S , 4 9 )| 2 

+ |X°(lQ,3„(2'6) ;;5 g ,4 g )| 2 

-2|M°(1q, (2~6)q, 5„ 4 9 , 3 7 )| 2 ) J? k%, fc 5 ) 

The squared matrix element for c/g — >■ QQqq in terms of colour ordered partial ampli- 
tudes is obtained by crossing A q and 6 q in eq.(7.39). In this case, the subtraction term 
reads 

g*{Nt - 1) 



d(7 qq^QQqq 



d^A{kiQ,k 2 Q,k 3q ,k 5q ;p Aq ,p 6q ) 
E$(4 q ; 3„ 1q) (|Mg((l~3) , 6 9 ; ; 5„ 2q)| 2 
+|Mg((13) , 1 5 , 6,; ; 5, , 2q)| 2 ) jffe fc 2 , fe 5 ) 



+£#(4,; 3„ 2q) |M»(1 , 6 9 ; ; 5„ 4 fl , (23)q)|< 



+|A^(1q, 4 s , 6,; ; 5, , (23)q)| 2 J jf^, k%, k 5 ) 
+El(4 q ; 5„ 1q) (|M°((1~5) Q , 6 g ; ; 3„ 4 fl , 2q)| 2 

+\Ml((15) Q , l g , 6 q ; ; 3„ 2q)| 2 ) J^fe fc 2 , fc 3 ) 



+E 3 U (4 ? ; 5, , 2q) (JA^Iq, 6 9 ; ; 3„ 4 fl , (25)q)^ 

+|_M°(l Q ,t s ,6 (? ; 53,, (25)g)| 2 ) Jf^fci.^fcs) 

+E$(6 q ; 5„ 1 Q ) (|Mg((f5)o, 4 g ; ; 3„ 6 fl , 2q)| 2 

+|.M°((1~5)q, 6 g , l q ; ; 3 q , 2q)| 2 ) 4 3 \k r5 , k 2 , k 3 ) 

+E$(6 q ; 3 q , 1 Q ) (\M° 5 ((13) Q ,l q ; ; 5,, 6 fl , 2q)| 2 

+|At°((f3) Q , 6 fl , 4 9 ; ; 5, , 2q)| 2 ) jffe fc 2 , fe 5 ) 



+£ 3 u (6 g ; 5„ 2 ) |^(1q,4,; ; 3„ 6 S , (25)$) 
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+\Ml(i Q ,6 g ,^; ;3„(2~5)g)| 2 ) 4 3) (fci,%,fe 3 ) 
+£ 3 °(6 (? ; 3„ 2q) (\Ml(l Q , 4 9 ; ; 5,, 6 fl , (2~3)g)| 2 

+|Mg(l ,6 ff , 4,; ; 5„ (2~3)q)| 2 ) jf ^fci, fc 5 ) 



+ |A^((13)q,2,,2q;;5 (? ,6 (? )| 2 

-2|.M°((1~3)q, 2q, 5„ 6„ 2 7 )| 2 ) jffe fc, fc 5 ) 
+£ 3 °(4 (? ; 3„ 2q) (|Mg(lg, (2~3) ; ; 5„ 4 fl , 6 g )| 2 

+|At°(l Q ,4 fl ,(23) ;;5 9 ,6 9 )| 2 

-2|Mjj(l g , (2~3)q, 5„ 6„ 2 7 )| 2 ) Jf (fci, fcg, fc 5 ) 
+£ 3 °(4 (? ; 5„ 1q) (|^((1~5)q, 2^; ; 3„ 6 g )| 2 

+|7Wg((15) C3 ,I 5 ,2 ;;3 <? ,6 <? )| 2 

-2|M°((f5) Q , 2 Q , 3„ 6„ 2 7 )| 2 ) jffe fc 2 , fc 3 ) 
+£ 3 °(4 (? ; 5„ 2q) (|Mg(l , (2~5) ; ; 3„ 6 q )\ 2 

+|^(l Q ,I 9 ,(2'5)Q;;3 g ,6 g )| 2 

-2|M°(l g , (2~5)q, 3 ff , 6 q , 2 7 )| 2 ) J® (fcx, fe^g, fc 3 ) 

+£ 3 °(6 (? ; 5„ 1q) (|^((1~5) q , 2q; ; 3 q , 6 g , l q )\ 2 

+ |^((r5)Q,6„2Q;;3„4 g )| 2 

-2|Mjj((f5) , 2q, 3 ff , 4„ 6 7 )| 2 ) jffe fc 2 , fc 3 ) 
+£ 3 °(6 (? ; 5„ 2q) (|Mg(l , (2~5) ; ; 3„ 6 g , l q )\ 2 

+\M° 5 (l Q ,6 g ,(25) Q ;;3 q ,l q )\ 2 

+|>t°(l Q ,(2~5) ,3 9 ,4 ? ,6 7 )| 2 ) Jf^fci.^fcs) 
+£ 3 °(6 (? ; 3„ 1q) (|^((1~3)q, 2^; ; 5„ 6 fl , 4 g )| 2 

+|7Wg((f3) Q ,6 p ,2 ;;5 q ,4 (7 )| 2 

-2\M° 5 ({13)q, 2 Q , 5 q ,4 q , 6 7 )| 2 ) jf (fejg, fc 2 , fc 5 ) 
+£ 3 °(6 (? ; 3„ 2q) (|Mg(lg, (23)^; ; 5„ 6 fl , 4 g )| 2 

+|^(l Q ,6 3 ,(23)Q;;5 9 ,4 g )| 2 

+|At°(l Q ,(23) ,5 9 ,4 9 ,6 7 )| 2 ) 4 3 \h,k^,h) 



E$(4 q ; 3„ 1q) (|Mg((13) Q , 2q; ; 5„ 6,)| 2 



7.2.3 Processes derived from — > QQqqgg 

The colour decomposition for the unphysical process — > QQqqgg is 
M°(l Q ,2 Q ,3 q A q ,5 g ,6 g ) 

= 2g 4 U,MT a ^) llM Mi(l Q ,i 9 ,j 9 ,4g-,-,3 q ,2 Q ) 
(i,j)eP(5,6) 1 

+(T a % lM {T a i )i 3 , i2 M° 6 (l Q , i g , 4,-; ; 3 q , jg, 2 Q ) 

+5 llM (T a 'T^) i3tl2 Ml(l Q Ag-,-,3 q ,i g ,j g ,2 Q ) 

-^(T a %^h iM Ml(l Q ,ig,2 Q] -,3 q ,jgAq) 



N, 



. Siu ^T«T*) iStU M%(lQ,2 Q y,3 q ,i g J g ,4z) 



and the partial amplitudes satisfy 
A^0(1q,2q,3 (? ,4 9 -,5 7 ,6 7 



Yl {M 6 (i Q ,ig,jgAq; ; 3„ 2 Q ) + M° 6 {i Q ,i g Aq, ; 3 q ,jg, 2 Q ) 

(i,j)6P(5,6) 

+M 6 {lQAq,;3q,ig,jg,2 Q )) 

Y {M° 6 (1q, kg, Ig, 2 Q ; ; 3 q , 4 ? -) + Mg(l g , fc fl , 2q; ; 3„ / s , 4 5 ) 



(fc,/)6P(5,6) 



+K(l ,2 ;;3 ? ,A; g ,/ ff ,4 3 -)). 



Squaring, and using the identity eq.(7.45) to rearrange the result, we obtain 

,2_5 8 (iV c 2 -l) 



\M"(l Q ,2 Q ,3 q ,4 q ,5g,6g) 



x E 

(ij)6P(5,6) 
r2 



(7.44) 



(7.45) 



JVc 2 ( lK(lQ, Ws> 5 3„ 2 Q )\ 2 + |Mg(lg, i fl , 4,-; ; 3 q ,j g , 2 Q )f 

+ \Ml{l Q ,A q ;;3 q ,ig,j g ,2 Q y 2 
+\M° 6 (l Q ,ig,jg, 2 Q ; ; 3 q ,A q )\ 2 + \M° 6 (1q, i g , 2 Q ; ; 3, , i fl , 4 9 -)| 5 



+ lK(lQ,2Q;;3 g ,i 9 ,j„4 ? -)p 



+ |7W|](lQ,2Q,3 (? ,4 (7 ,i 7 ,j 7 )| 5 



+N, 



2Re{M 6 (l Q , i g , j g , 4 q ; ; 3, , 2 Q )M° 6 (l Q ,4 q ; ; 3„ i g ,j g , 2$)+) 

+2Re(M°(l Q , i g , j g , A q ; ; 3 q , 2 Q )M° 6 (l Q ,4 q ; ; 3 q ,j g , i g , 2 Q ?) 
+Re(M° 6 (l Q ,i g A q ; ; 3 q ,j g , 2 Q )M° 6 (l Q ,j g A q ; ; 3„ i s , 2q)+) 
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-Re(M° 6 (i Q ,i g ,j g Aq; ; 3„ 2 Q )M 6 {i Q ,jg,i g Aq; ; 3„ 2 Q y) 

-Re(M° 6 (l Q ,4g; ; 3„ i g ,j g , 2q)M° 6 (1q, 4 q ; ; 3„ j 9 , 2 Q )t) 
-2Ite(Mg(l , ; 3„ 2q).M°(1q, i 9 , j 9 , 2 Q ; ; 3„ 4 ? -)t) 

-2Re(M°(l Q , i g , 4 q ; ; 3„ j g , 2 Q )M° 6 (l Q ,i g , j g , 2 Q ; ; 3,, 4 ? -))t) 
-2Re(M° 6 (l Q Aq; ; 3„ 2g)Mg(l Q , i g ,j g , 2 Q ; ; 3„ 4 ? -)t) 
-2Re(M° 6 (l Q ,i g ,j g Aq, ; 3„ 2q)A^°(1q, 2 Q ; ; 3„ j 9 , 4 ? -)t) (7.46) 
-2Re(M° 6 (l Q ,i g Aq, ; 3„j g , 2q)A^°(1q, i g , 2 Q ; ; 3„ j ff ,4 ? -)t) 
-2Re(M 6 (l Q ,i g ,4q; ; 3„i 9 , 2q).M°(1 q , j g , 2 ; ; 3„ i g , 4^) 
-2Re(M°(l Q , 4,-; ; 3„ 2q).M°(1 q , j g , 2 ; ; 3 q ,i g , 4 ? -)t) 
-2Ite(Mg(lQ, 4,-; ; 3„ 2 (5 )Mg(l Q , 2 ; ; 3„ j 9 , 4 ? -))t) 
-2Re(.M°(l Q , 4,-; ; 3„ 2 (5 )Mg(l Q , 2^; ; 3 q ,i g , j g A^) 

-2Re(M° 6 (lQ,ig,4q, ; * q ,3g, ^q)J<(Iq,2 q - ; 3„ j 9 , i 9 , 4 ? -)t) 

+Re(M° 6 (l Q , i g , 2 Q ; ; 3„j ff , 4 q )M° 6 (l Q , j g , 2 Q ; ; 3„ i g , 4 q )^) 
-Re(Mg(l Q , 2 ; ; 3„ 4 q )M° 6 {l Q , j g , i g , 2 Q ; ; 3 q ,4q) f ) 
-Re(M° 6 (l Q , 2 Q ; ■ 3 q ,i g ,jg,4q)M° 6 (lQ, 2 Q ; ; 3 q , j g , i g , 4 g -) t ) 
+2Re(.M°(l Q , 2q; ; 3„ 4,-)M°(l Q , 2 ; ; 3„ i g ,j g A^) 



Besides colour-ordered amplitudes squared this expression includes interference terms 
between amplitudes with different colour orderings which require subtraction. 

After appropriate crossings, for each of the colour-ordered matrix-element squared ap- 
pearing on the right hand side of this equation, an antenna capturing the single unresolved 
radiation is determined. For the soft radiation behaviour present in the interference terms, 
the subtraction terms are constructed with the difference of four antennae as explained 
in Section 7.2.1. For the subtraction terms, we will need all types of massive final-final 
and initial final antennae defined in Section 3 in addition to massless antennae in all three 
(final-final, initial-final and initial-initial) configurations which were defined in [42,57,60]. 

By crossing 3 q and 4 q in eq.(7.46) to the initial state, the matrix element for the process 
qq —7- QQgg is obtained. From the unresolved limits of this squared matrix element we 
obtain the following subtraction term, 



+2Re(M° 6 (l Q ,ig,jg, 2 Q ; ; 3„ 4 q )M° 6 (l Q , 2 Q ; ; 3 q A 9 ,i 9 AqY) 





AT c 2 [^(4 9 ;1q,* 9 )|^^ 

+A° 3 (3 q ; 2 Q , i g ) | Mi (1q , j g , 4 q ; ; % q , (2i) Q ) | 2 jf ( k x , %. , kj ) 

+ 4{l Q ,ig,jg)\Ml{{li) Q ,{ij)gA q V 
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+dg(2Q,i a ,i s )|Alg(lQ,4 g ;;3 ? ,(^) s ,(2O O )|V| 3) (fc 1 ,%.,^0 

+d°(4 9 ;j fl ,i fl )|AlO(lQ,(^0 5 ,4 9 ;;3 9 -,2Q)|Vf(A; 1 ,A; 2 ,^0 

+c$(3 g -;^ 5 )|A<°(lQ^ 

+A§(lg, i fl) 2q) (\M° 5 ((U) q , (2i)g; ; 3 ? , j s , 4„)| 2 

+2Re(^((H) Q , 4^; ; j g , (2z) )M°((ll) Q , j g , 4,; ; 3 5 , (2i) Q ?) 
-2Re(Ml((li) Q , {2i) Q - ■ 3 q ,j g , l q )Ml((li) Q ,l q ; ; 3 q ,j g , (2i)$)t) 
-2Re(M°((ll) Q , (2% ■ 3 q ,j 9 , 4 q )M° 5 ((li) Q , j g , 4 ? ; ; 3 5 , (2l) Q )t)) jf (fcj,, fc s , fc,-) 

-A°(4 g ;l Q ,z g ) (lA^OtCr^Q,^,^;;^^^)! 2 + | ((ri) Q , I Q ; ; 3 g -, 2^) | 2 

+2Re(M° 5 ((U) Q , l q ; ; 3 q ,j g , 2 q )M° 5 ((U)q, 2 q - ; 3 g -, j g , 4,)*) 
+2Re(At°((H) Q , 4,; ; j fl) 2 )M°((ll) Q , j g , 2^; ; 3,-, 4,)t)) (fcj-., fc 2 , %) 

+A°(3 9 -; 1 Q , (2Re(M°((H) Q , 4 g ; ; 3 5 , j g , 2q)M°((H)q, j flJ 2 Q - ; 3 ffj 4,)t) 

+2Re(Mg((fi) , j g , 4,; ; 3 9 -, 2 ( g)Mg((fi) g) 2 Q ; ; j g , 4,)t) 

-2Re(Mg((fi) , j 5 , 4 9 ; ; 2q)M°((H)q, 4 g ; ; 3 q ,j g , 2^)) jf (fcj-., fc 2 , 
+A°(4 ? ; 2q, i fl ) (2Re(A^°(l Q , 4,; ; 3 ? , j fl) (2z)q)M°(1q, (2i) Q ; ; 3 ff , j fl , 4,)t) 

+2Re(.M°(l , j s ,4,; ; 3 q , (2i) Q )M 5 (l Q ,j g , {2i) Q ; ; 3 q , 4,)t) 

-2Re(7W°(l Q ^;; 3^,(2^ 

-A°(3 g -;2Q,z g ) (lyUgClQ^^;^,^^^)^)! 2 + |A^g(lc3 , i 5 , 4 g ; ; 3 9 -, (2i) ) | 2 

+2Re(At°(l Q ^,4,;;3 9 -,(^ 

+2Re(At°(l Q ^,4,;;3 ? -,^ 
+A°(3 9 -, 4 g ; i p ) J fl) % ; 3 q , l q )\ 2 

+2Re(Ml(l Q , l q ; ; 3 q ,j g , 2 Q )M° 5 0-Q,jg, ^ 5 2g) f ) 
-2Re(.M°(l Q , 4,; ; 3 9 - J g , 2 )M°(1 Q J g , % ; 4 9 )t) 
-2Re(.M°(l ~j g ,l q \ ; 2q)^(1q, J 9 , % ; 3,-, 4,)*)) J< 3) (*i, fc, k) 
+d° 3 (l Q ,i g ,j g ) (\M° 5 ((U) Q , (ij) g , 2 Q ; ■ 3 q , 4 g )| 2 

+2Re(M°((ll) Q , 4 g ; ; 3 5 , 2q)^°((H)q, 2^; ; 3 ? , 4*) f ) 
-2Re(M°((ll) Q , 4,; ; 3 q , (ij) g , 2 Q )M° 5 ((U) Q , (ij) g , 2 Q ;^ q , 4 g )t) 
-2Re(Mg((li) 0) 4,; ; 3 g -, 2 g ).M°((H) Q , 2^; ; 3 9 -, 4,)t) 
-2Re(M°((ll) Q , 4,; ; 3 9 -, 2 g ).M°((H) Q , 2q; ; 4,)t)) jf , fc 2 , fc..) 

+rf°(2Q, Ws) (\M° 5 (1q, (ij) g , (2i) Q ; ■ 3 q , l q )\ 2 
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+2Re(M° 5 (l Q , (ij) g ,4 q ^3 q , (2i) Q )M° 5 (l Q , (2% ; 3 q , (ij) g , 4 g )t) 
-2Re(_M°(l Q , 4 g ; ; 3 9 -, (2i) )X°(l Q , (2% ; 3 ? -, 4,)t) 
-2Re(.M°(l Q , 4 g ; ; 3 ? -, (2i) )^(l Q , (y) fl , (2i) Q ; ; % 4 g )t) 
-2Re(M°(l , 4,; ; 3 q , (2z) )M°(l Q , (2z)q; ; 3,-, 4,)t)) jf fcg, fc.) 

+d°(4,;i ff ,j fl ) (|X0(l Q ,2Q;;3 g -,(i7-) 9 ,3 g )| 2 

+2Re(M°(l , 4,; ; 3,-, 2q)^(1q, (ij) g , 2 Q - ; 3 5 , 4,)t) 

-2Re(Ml(l Q ,l q ; ; 2q)^(1q, 2 q ; ; 2 9 ) f ) 

-2Re(M° 5 (l Q , (ij) g , 4 ? ; ; 3 ? -, 2q).M°(1q, 2^; ; 3 ? , 4,) f ) 

-2Re(M°(l Q , l 9 ; ; 3 ? -, 2q).M°(1q, (ij) g , 2 Q ; ; 3 5 , 4,)+)) 4 3 \h,k 2 , k~) 
+d° 3 (3 q ; i g ,j g ) (\M° 5 (1q, 2 Q - ; 3,-, (ij) g ,l q )\ 2 

+2Re(.M°(l Q , (ij) g ,4 q ; ; 3 q , 2 Q )M° 5 (l Q , (ij) g , 2 Q ; ; 3 q , 4,)*) 
-2Re(_M°(l Q , 4 q ; ; 3 9 -, 2 <5 )A<°(1q ) {ij) g , 2 Q ; ; 3,-, 4,)t) 
-2Re(_M°(l Q , 4 g ; ; 3 ? -, 2 < g)Mjj(l , 2 _. . 4,)t) 
-2Re(.M°(l Q , ^%,4 ? ; ; 3 ? -, 2q)M°(1 q , 2^; ; % 4,)*)) jf ^ k 2 , k~) 

AUlq, i g , 2q) (2\M° 5 ((U)q, (2i) Q , 3 q , 4 g , j 7 )| 2 - \M%{{U) Q , (2i) Q - 3 q , j g , l q 



1 



Vl 2 



-\M° 5 ((li) Q ,j g , (2i)g; ; 3 5 , 4 g )| 2 ) jf (fe r ., fc 2l , fc,-) 



+A§(3,-, 4 g ; i g ) (2\M° 5 (1q, ~2q, 3,-, 4„ j 7 )| 2 - \M° 5 (1q,2 q ; ; 3 ? -, j 9 , 1, 
-|^(lQ,i fl ,2Q;;3,-,4 9 )| 2 ) jf\h,~k 2 ,k) 

;3g,4g)^)jg ^ ( 1 , fc^- , j ) 



Q 

-^(3 g -;l ,z g )2Re(Atg((H) Q ,2 ;;3,-,i g ,4,)^((H) Q ,i g ,2g 
-^(4,;2 ,z g )2Re(^(lQ,(20 ;;3,-,i g ,l,)^(lQ,i g ,(20Q 



When 4 g and 6 g are crossed to the initial state in eq.(7.46), the squared matrix element 
for qg — > QQqg in terms of colour ordered subamplitudes is obtained. The subtraction 
term in this case reads, 



d& qg^QQqg = 9 & ( N c ~ l)d$ 4 (&iQ, k 2Q , k 3q , h g ;p 4q , P(ig ) 



N 2 

1 v c 



Ai(3 q , 5 g , 2 Q )\Mi(l Q ,6 g , 4,; ; (35),, (25)q)| (h, ftgg, fcgg) 
+A°(4,; 1 Q , 5 fl )|M°((l~5) Q , l ? ; ; 3,, 6 fl , 2q)| 2 jf \k%, k 2 , fe 3 ) 
+^§(6 fl ; 1q, 2q) (|Mg((f2) Q) 5 9 , 4 9 ; ; 3„ 6q)| 2 + |7Wg((12) Q , 4 g ; ; 3 q , 5 g , t Q )f 
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+\M° 5 $ Q , 5 g , l q ; ; 3„ (1~2)q)| 2 + \M° 5 $q, \\ 5 3„ 5 g , (1~2)q)| 2 ) J 3 (3) (%, fc 3) k 5 ) 
+A° 3 (4 q , 6 g ; 3 q ) (\M° 5 (1 Q ; 6 ?) 5 fl , %)| 2 + |A^(1q, 5 fl , I,; ; 6,-, ~2 Q )\ 2 ) jf\hMM) 
+D° 3 (6 g ; 5 g , 1 Q )\M° 5 ((15) Q , %,4 q ; ; 3„ 2q)| 2 jf fc 2 , fc 3 ) 
+D 3 °(6 S ; 5 fl , 2 Q )\M° 5 (l Q ,4 q ; ; 3„ 6 g , (2~5) )| 2 jf fc 3 ) 
+D 3 °(6 9 ; 5 fl , 3,)|Mg(l„ 4,; ; (35), ,B ff , 2q)| 2 jf \k u k 2 , k~) 
+D° 3 (4 q , 6 g ; 5 g )\M° 5 (l Q , 6 g ,l q ; ; 3„ 2q)| 2 jf (fci, fc 2 , fcs) 
+^ 3 °(4 ? ; 3„ 1q) (|Mg((f3) Q , 5 9 , 6 ff , l 9 , 2q)| 2 + \M° 5 {(13)q, 5 fl , 6 ff , 2q)| 2 
+ |^g((ll) Q ,l 9 ,5 9 ,6 9 ,2 ( 3)| 2 + |^((r3)Q,6 9 ,5 9 ,l 9 ,2Q)| 2 
+|.M°((1~3) Q , 6 9 , 5 9 , 2q)| 2 + |Mg((f3) , l 9 , 6 9 , 5 g , 2 )| 2 ) jf k 2 , k 5 ) 

+^ 3 °(4 ? ; 3„ 2 ) (|Mg(l Q) 5 9 , 6 s , l s , (23)q)| 2 + \M° 5 {1q, 5 g ,l g , 6 g , (23) <5 )| 2 

+|Mg(l , 6 ff , 5 9 , (23)q)| 2 + \M° 5 {l Q ,Zg, 5 g , 6 g , (2~3)q)| 2 

6 S , 4 fl , 5 ff , (23) )| 2 + \M° 5 (1q, l g , 6 g , 5 g , (23) )| 2 ) jf } (fci, fc 5 ) 

5 g , 2q) (|Mg((f5)(3, (2~5)q; ; 3 g , 6 g , l q )\ 2 

+2Re(M°((f5) Q , 4 5 ; ; 3„ 6 a , (2~5)q)M°((1~5)q, 6 9j 4 g ; ; 3„ (25) <3 )t) 
-2Re(Mg((f5) , 4,; ; 3„ 6 g , (2~5) ).M°((1~5)q, (2~5)q; ; 3 g , 6 g , 4,)t) 

-2ReCM°((l~5) Q) 6 g , 4 g ; ; 3„ (25) < g)Mg((r5) , (2~5) ; ; 3 g , 6 9 , 4 g )t)) jf fcg, fc 3 ) 
+A°(1 , 5 g , 3,) (2Re(Mg((f5) , 4 ? ; ; (35),, 6 fl , 2 <5 )M°((f5) Q , 6 fl , 2^; ; (3~5) 9 , 4,)t) 

+2Re(^°((f5)Q, 6 fl , 4 9 ; ; (3~5)„ 2 )^((1~5) Q , 2^; ; (3~5)„ 6 g , 4,)t) 

-2Re(M°((l~5) Q , 4,; ; (3~5)„ 6 3 , 2 ).M°((1~5)q, 6 9 , 4,; ; (35),, 2g)t)) jffe k 2 , k~) 
-A° 3 (3 q , 5 g , 2q) 4 3 ; ; (35),, 6„ (2~5) )| 2 

+|K(lQ,6 ff ,4 ? ;;(3~5), L (2~5)^| 2 

+2Re(A^°(l Q , 6 9 , 4 g ; ; (35),, (25) Q )M° 5 (1 Q , (25) ; ; (35),, 6 g , 4,)t) 
+2Re(X°(l Q , 6 ff , 4,; ; (35),, (25) <5 )^g(l , 6 ff , (2~5) ; ; (35),, 4,)t)) jf (k ± , k%, k~) 
-A° 3 (A q ; 1 Q , 5 g ) (|Mg((f5) Q , 6 S ,4 ? ; ; 3„ 2q)| 2 
+ |^°((f5)Q,l (? ;;3 (? ,6 9 ,2 ( 3)| 2 

+2Re(At°((f5) Q , 4,; ; 3„ 6 9 , 2 )At°((f5) Q , 2^; ; 3„ 6 9 , 4,)t) 
+2Re(X^((15)Q,Iq; ; 3„ 6 g , 2 Q )M° 5 ((15) Q , 6 g , 2 Q ; ; 3„ 3,)t)) (fejg, fc 2) fc 3 ) 

+A° 3 (A q ; 2 Q , 5 g ) (2Re(M° 5 (l Q , l q ; ; 3„ 6 9 , (2~5)q)A^°(1q, (2~5) ; ; 3„ 6 ff , 4,)t) 
+2Re(A^°(l Q , 6 9 , l q ; ; 3„ (25)q)A^°(1q, 6„ (2~5) ; ; 3„ 3,)t) 
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-2R*(M°(1 Q , l q - ; 3 ? , 6 g , (2~5) )>t°(l Q , 6 g , 4,; ; 3„ (25)<g)t)) J^(fei, fegg, fc 3 ) 

+^(4 g ; 3„ 5 fl ) 6 g , 2q; (35),, I,)) 2 

+2Re(M° 5 (l Q , l q ; ; (3~5) ? , 6 5 , 2 Q )M° 5 (1 Q , 6 g , l q ; ; (35),, 2 Q ?) 
-2Re(M° 5 (l Q , l q ; ; (3~5) ? , 6 g , 2 Q )M° 5 (l Q ,6 g , 2 Q ; ; (35),, 4,)t) 
-2Re(Mg(l , 6 g , t q ; ; (35),, 2 Q )M° 5 (l Q ,6 g , 2 Q ; ; (35),, t,)t)) 4 3 \h,k 2 , k~) 

+^s(6 9 ; lQ, 2q) (|^((1~2)q, 6q; ; 3„ 5 ff , 4 g )| 2 + \M° 5 ((12) Q , 5 g , 6 Q ; ; 3„ 4 g )| 2 
+ |A^O(6 , (1~2) ; ; 3„ 5 g , l q )\ 2 + \M° 5 $ Q , 5 g , (12) Q ; ; 3„ 4 g )| 2 
-\M° 5 ((12) Q , l q - ; 3„ 5p, 6 Q )| 2 - |Mg((f2) , 5 g , 4,; ; 3^ 6 Q )| 2 
4 ? ; ; 3 q , 5 g , (12) Q )\ 2 - \M° 5 $Q, 5g, ^ \ 3„ (1~2)q)| 2 
-2|M°((1~2) , 6q, 3„ 4„ 5 7 )| 2 - 2\M° 5 (t Q , (12) Q , 3 q ,4 q , 5 7 )| 2 ) J^fe fc 3 , fc 5 ) 

+A§(4,, 6 9 ; 3,) (|Mg(l , 5 ff , % ; 6 ? -, l 9 )| 2 + |M°(1 Q , % ; 6 q , 5 g , l q )\ 2 

-\M° 5 (l Q , 5 g , 4 ? ; ; 6,-, 2q)| 2 - |Mg(l , 4,; ; 6 5 , 5 9 , %)| 2 

-2|Mg(l , 2q, 6 ? , 4„ 5 7 )| 2 ) Ji(h, ~k 2 , h) 
+D° 3 (6 g ; 5 g , 1 Q ) (\M° 5 ((15) Q , 6 g , 2 Q ; ; 3„ 4^)1 2 

+2Re(M° 5 ((lh) Q ,4 q ; ; 3„ 6 g , 2 Q )M° 5 {(15) Q , 2q; ; 3„ 6 g , 4,)t) 
-2Re(M°((l~5) Q , 4 g ; ; 3„ 6 g , 2q).M°((1~5)q, 6 g , 2^; ; 3„ 4,)t) 
-2R*(Mjj((f5) , 6 ff , 4 g ; ; 3„ 2q).A^((1~5)q, 6 ff , 2 Q ; ; 3„ 4,)t) 
-2Re(.M°((l~5) Q , 6 g , l q - ; 3„ 2 )^((1~5) Q , 2 ; ; 3„ 6 g ,4 q )+)) jffe fc 2 , fc 3 ) 
+£ 3 °(6 9 ; 5 9 , 2 ) 6 g , (2~5) ; ; 3„ 4,)| 2 

+2Re(M° 5 (l Q ,6 g ,4 q ; ; 3„ (2~5) )^(1 Q , (2~5) ; ; 3„ 6 ff , 4,)t) 
-2Re(M°(l Q , 4 g ; ; 3„ 6 9 , (2~5)q)^(1q, (25) q ; ; 3„ 6 9 , 4,)t) 
-2Re(Ml(l Q , 4 q ; ; 3„ 6 ff , (25) Q )M° 5 (1 Q , 6 g , (25) Q ; ; 3„ 4 g )t) 
-2Re(^°(l Q , 6 ff , 4 g ; ; 3„ (2~5)q)^(1q, 6 9 , (2~5) ; ; 3„ 4,)t)) jf ^fci, fc 3 ) 
+J D 3 (6 9 ; 5 g , 3 q ) [\M° 5 (1 Q , 2 q ; ; (35),, 6 g , l q )\ 2 

+2Re(M° 5 (l q , 6 g ,4 q ; ; (35),, 2q)^(1q, 6 ff , 2^; ; (35),, 4,)t) 
-2Re(A^°(l g , 4 g ; ; (3~5) 3 , 6 g , 2 Q )M° 5 (1 Q , 2 Q ; ; (35),, 6 g , 4,)t) 
-2Re(A^°(l g , 4 g ; ; (35),, 6 fl , 2q)^(1q, B 9 , 2^; ; (3~5)„ 4 g )t) 
-2Re(A^°(l g , 6 9 , 4 g ; ; (35),, 2q)^(1q, 2 ; ; (35),, 6 g , 4,)t)) jf fc 2 , fegg) 
+J D 3 °(4 g , 6 9 ; 5 9 ) ; 3 g , 6 g , l q )\ 2 
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+2Re(M°(l Q , 4 q ; ; 3„ 6 9 , 2 Q )M° 5 (l Q ,6 g , % ; 3„ 4,)*) 
-2Re(A^0(l Q , l g ; ; 3„ 6 a , 2q)A*°(1q, % ; 3„ 6 9 , t,)t) 
-2Re(M°(l Q , 6 g , 4 q ; ; 3„ 2q)^(1q, % ; 3„ 6 ff , 4,)*) 



-2Re(Mg(l Q , 6 g , 4 q ; ; 3„ 2 Q )M u 5 (l Q ,6 g , 2 Q ; ; 3„ 4,)T) J jWfa, k 2 , k 3 ) 



-E»(4 q ; 3„ 1 Q ) (|A^((13)q, 5 9 , 6 9 , 4 7 , 2q)| 2 + |A^((13) , 6 9 , 5 g , 4 7 , 2 Q )\ 2 
+|Mg((f3) Q , 5 9 , t 9 , 6 7 , 2q)| 2 + |Mg((l~3) , t ff , 5 g , 6 7 , 2q)| 2 
+|M°((1~3) Q , 6 9 , 4 ff , 5 7 , 2q| 2 + |M°((1~3) Q , 4 fl) 6 9 , 5 7 , 2q| 2 ) jf \k^, k 2 , k 5 ) 



-E»(4 q ; 3„ 2q) (|Mg(l , 5 9 , 6 fl , 4 7 , (23)g)| 2 + |^(1q, 6 g , 5 g , 4 7 , (23)q)| 2 
+|Mg(l , 5 ff , 4 fl , 6 7 , (23) )| 2 + |Mjj(lg, 4 fl , 5 ff , 6 7 , (2~3)g)| 2 
+|Mg(l , 6 g , tg, 5 7 , (23) )| 2 + \M° 5 (1 Q , l g , 6 g , 5 7 , (23) )| 2 ) jf (fci,^, fc 5 ) 



+ 



A°(l , 5 9 , 2q) (2|A<°((15)q, (25)g, 3„ 4„ 6 7 )| 2 



■|Mg(15) Q , (25) ; ; 3„ 6 fl , 4 ? )| 2 - |A^(15)q, 6 ff , (25) ; ; 3„ 4 9 )| 2 ) jf fo^, fc 3 ) 



+^(3 ? ,5 ff ,2Q)2Re(At°(l ,(2~5)Q 
-^(lg,5 p ,3,)2Ite(Mg((r5)o,2 < 5 
+A§(4,;l ,5 p )2R*(Mg((f5)Q,2Q 
-A°(4 g ;2Q,5 fl )2Re(M°(l Q ,(2~5)Q 



(35),, 6 9 , 4,)M°(1q, 6 g , (25) ; ; (35),, 4^) jffa, kg, fegg) 
(35),, 6 fl , 4 q )M° 5 ((15) Q , 6 g , 2 Q - ; (3~5)„ 4,)t) jf fc 2 , fegg) 
3 q , 6 g , 4 q )M° 5 ((15) Q , 6 g , 2 Q - ; 3 q , 4,)t) jf } (fcjg, fc 2 , fc 3 ) 
3„ 6 fl , 4 9 )M°(1 , 6 9 , (2~5)q; ; 3„ t,)t) jf } fc 3 ) 



+A«(4, ; 3„ 5 fl ) 2|^(1 Q , 2$, (35) 9 , 4 q , 6, 



6 5 , 2q; (35),, 4 g )| 2 - \M° 5 (1q, 2 q ; ; (35) 9 , 6 fl , 4,)| 2 ) jf \k u k 2 ,k^) 



+-A»(6 g ; 1 Q , 2q) (2|Mg((12) , 6 Q , 3„ 4„ 5 7 )f + 2\M u 5 {6q{12) q , 3 q ,4 q , 5 7 )| 2 
-|Mg((f2)0, 5 9 , 6 Q ; ; 3„ 4,)| 2 - \M° 5 (6 Q , 5 g , (1~2) ; ; 3„ 4 g )| 2 
-|Mg((f2) Q , 6 Q ; ; 3„ 5 9 , 4,)| 2 - |A^°(6q, (1~2) ; ; 3„ 5 9 , 4,)| 2 ) (feg, fc 3 , fc 5 ) 

+^(4,, 6 9 ; 3,) (2|Mg(l Q , 2 , 6 ? -, 4 q , 5 7 )| 2 

5 ff , 2 ; ; 6 q , l q )\ 2 - \M° 5 (1 Q , % ; 6 ? -, 5 9 , l 9 )| 2 ) j£\hM,h) 
+^ 3 °(4 ? ; 3,, lg)|Mg((13) , 2g, 4 7 , 5 7 , 6 7 )| 2 jf^, fc 2 , fc 5 ) 
+^ 3 °(4 ? ; 3„ 2 )|^(1q, (2~5)q, 1 7 , 5 7 , 6 7 )| 2 jf (fci, h) 

Finally, when gluons 5 g and 6 g are crossed to the initial state in eq.(7.46), the squared 
matrix element for gg — > QQqq is obtained. This partonic process only contains collinear 
and quasi-collinear limits but no soft limits as no gluons are present in the final state. 
Using the decoupling identities given in eqs.(7.6, 7.12, 7.13), the reduced matrix elements 
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multiplying each antenna function can be rewritten in a fairly compact form with no 
interference terms left. After doing this simplification the subtraction term for this process 
reads, 



d& gg^QQgg = 9 Nf ( N c ~ l)d$ 4 (fciQ, k 2 Q, k 3q , k 4q ;p 5 g,P6g) 

*< £ 

(i,j)eP(5,6) 



+ 



^°(* fl ;l Q ,2 )(|At^ 

+\M° 5 ((i2) Q ,j g ,Ag;;S q ; iQ )\ 2 + \MlCi Q 3g^y^(i2) Q )\ 2 )j { i\k T2 MM) 

+ iA°(i s ;3 g ,4 (7 )(|>(0(lQ,! g ;;(34) 9 ,i,,2Q)| 2 + (34) 5 ; 2q)| 2 

+\M° 5 (1q, j g ,i q ; 5 (34),, 2 )| 2 + \M° 5 (l Q ,j g , (34),-; ;i q , 2 Q )\ 2 ) jf\k u k 2 ,k^) 
+^3°(lQ,3 ? ,4 ? -) (|A<g((f3) J ff J ff ,(3l^ 

MMi{{n) Q ,{M) g ;t g ;j g ,2 Q \ 2 ) j { i\k Tv k 2 ,k~) 

+ ^2q,3,,4 5 )(|A^(1 q ,^ 

MMl(l Q AM) g ; ig ,j g A23) Q )\ 2 )4 3 \k 1 ,k T3 ,k r4 )_ 

±A° 3 (i g ;l Q ,2 Q ) (|^((r2) Q ,fQ;;3 g ,i,,4 g -)| 2 + |^(!Q,(r2) ;;3,,j 9 ,4 g -)| 2 

+ |^((r2)Q,j 9 ,iQ;;3 g ,4 5 )| 2 + |^(!Q,i,,(r2)Q;;3 (? ,4 g -)| 2 
-|^0((r2)Q,4 5 ;;3 9 ,i 9 ,f Q )| 2 -|^(!Q,4 5 ;;3 g ,i 9 ,(r2) )| 2 (7.47) 
-|^((r2)Q,j 3 ,4 g -;;3 g ,f Q )| 2 .-|^(iQ,i s 4 g -;;3 g ,(r2)Q)| 2 
-2|^0((r2)Q,f Q ,3 g ,4 g -,i 7 )| 2 -2|^(f ,(r2)Q,3 9 ,4 g -,i 7 )| 2 ) jffe fc 3 , & 4 ) 

+ l -Al(t g - 3 ? ,4,-) (|Mg(l , 2q; ; (34),, j,,!,)! 2 + |M°(1 Q , 2^; ;i q ,j g , (34) 5 )| 2 
+|Mg(lg J fl , 2 ; ; (34)„i,)| 2 + \M° 5 (l Q ,j g , 2 Q ;;i q , (34) 5 )| 2 
-\M° 5 (l Q ,i q ; ; (34) 9 Jp, 2 )| 2 - \M° 5 (1 Q , (M) q ;;i q ,j 9 , 2 Q )\ 2 
-\M° 5 (l Q ,j g ,i q ; ; (34),, 2 )| 2 - \M° 5 (l Q ,j g , (34),-; ; 2q)| 2 
-2|^0(lQ,2Q,(3l) g ,i g ,^)| 2 -2|^°(lQ,2 ,l 5 ,(3l)ff,j 7 

-^ 3 °(1q, 3„ 4 5 ) (|-M°((f3) , i fl , j 9 , (34) 7 , 2q)| 2 + \MU(&) Q ,i g ,(M) g ,h,2 Q )\ 2 

+ \M° 5 ((13) Q , (34) 9 , i fl , j 7 , 2q)| 2 ) jf (%, fe, fcg) 
-^3°(2Q,3 g ,4 5 ) (|A^(l Q ,i 9 J 9 ,(3l) 7 ,^ 



+|M^(l0,(34) fl ,i p J 7 ,(23) (3 )| :4 )4 3) (A;i,*53,fc 
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2JV? 



A° 3 (i g ;l Q ,2 Q ) (|^°((12) ,! Q ;;3,J g ,4,-)| 2 + |>f°(!g,(12) ;;3 ? ,i ff ,4^ 



+ |^0((12)Q,j 9 ,i Q ;;3 9 ,4 (? )| 2 + |^(iQ,i,,(12) (3 ;;3 g ,4 g -)| : 
-2|^((r2)Q,lQ,3 g ,4 g -,i 7 )| 2 -2|^(i ,(r2)Q,3 g ,4 9 -,j 7 )| 2 ) jf\k^k 3 M) 
+A° 3 (i g ;3 q ,4g) (|.Mg(l Q ,2Q; ; (34) g , i 5 , l Q ) | 2 + 2^; (34) g -)| 2 

+\M° 5 (l Q ,j g , 2 Q ; ; (34) g ,I g )| 2 + \M° 5 (l Q ,j g , 2$; (34),-)| 2 
-2|^(lQ,2Q,(34)„i g ,i 7 )| 2 -2|^(lQ,2Q,f g -, (34) 9 -,j 7 )| 2 ) 4 3) (h,k 2 ,k 

j - (^ 3 ^l g ,3 g ,4 9 -)|>t°((ll)Q,(3l) 7 ,^J 7 ,2Q)| 2 jf ) (^,A ;2 ,A ; 3 2 ) 
+^3^2Q,3 g ,4^|^(l Q ,(3l) 7 ,i 7 ,i 7 ,(23) )| 2 jf (h,k T3 ,k r4 ) 



N. 



For all three crossings considered, we have performed the consistency check explained in 
Section 7.1.3. We have checked that in all collinear or quasi-collinear limits contained in 
these subtraction terms, the sum of the terms contributing to a given limit collapses to 
the appropriate colour factor multiplying a given massless or massive splitting function, 
corresponding to the given collinear splitting, times the appropriate non-colour ordered 
matrix-element squared. 

7.2.4 Partonic process gg — > QQgg 

The colour decomposition for the unphysical process — > QQgggg is 

M 6 °(1q,2q,3 3 ,4„5„6,) = 4/ £ (T a ^T a "T a % tl2 M(i,j,k,l), (7.48) 

(iJ,fc,0eP(3,4,5,6) 

where we used A4(i,j, k,l) = A4q(1q, i g , j g , k g , l g , 2q) for simplification of the formulae. 
Squaring gives 

K(1 q ,2q,3„4„5„6,)| 2 = 5 ' ; X ; ~ 1 



E 

(iJ,fc,0GP(3,4,5,6) 



N 2 



N?\M(iJ,kJ)\ 2 -N*\M(iJ,k;;l)\ 2 + 1 -^\M(i,j;;k,l)\ 2 



N*Re^(M(j, i, /, k) +M(j, I, i, k) +M(j, I, k, i) 
+M(k,i,l,j) +M(k,j,l,i) +M(l,i,k,j) 
+M(l,j,i,k)+M(l,k,j,i)) x M{i,j,k,l) ] 
(Ar 4 _ 3AT2 _ 1 



(7.49) 



■\M(i,j,k,l)\< 



4! 

where when the gluon labelled with I, when the gluons labelled with I and with k or when 
all four gluons decouple we have respectively, 



M(i,j,k;;l) = Mq(1q, i g ,j g , k g , Z 7 , 2q) 



(7.50) 
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= M{i,j,k,l) + M{i,j,l,k) + M{i,l,j,k) + M{l,i,j,k) 

M(i,j;;k,l) = M%{1q, i g , j g , k 7 , Z 7 , 2q) 

= M(i, j, k, I) + M(i, j, I, k) + M(i, k, j, I) +M{i, l,j, k) (7.51) 
+M(i,k,l,j)+M{i,l,k,j)+M(k,i,j,l)+M(l,i,j,k) 
+M{k,i,l,j)+M{l,i,k,j)+M{k,l,i,j)+M{l,k,i,j) 

Ai(3,4,5,6) = K(1q,3 7 ,4 7 ,5 7 ,6 7 ,2 ) (7.52) 

Yl Ms(lQ,ig,jg,k g ,lg,2Q) 

(j,i,fc,0eP(3,4,5,6) 

This colour-ordered decomposition for the squared amplitude for gg — > QQgg has also been 
derived in [87]. Our result is in complete agreement with it. 

Crossing gluons 5 9 and 6 g in eq.(7.49) gives the squared matrix element for gg — > QQgg 
reading, 

|M 6 °(l Q ,2 ,3 g ,4 g ,5 g ,6 g )| 2 = g^=^ 
x < > \n® (\M(i,j, k, l)\ 2 + \M(i,k, j,l)\ 2 + \M(i,k,l,j)\ 2 



-N 4 



(*J)6P(3,4),(fc,i)6P(5,6) L V 

+\M(k,i,ji)\ 2 + \M(k,i,i,j)\ 2 + \M(k,i,i,j)\ 2 

-^ c 4 (|^(i,i,fc;;0| 2 + |^(i,^j;;0| 2 + |^(i,fc,/;;i)| 2 
MM(k,z,jy,l)\ 2 + \MCk,iJy,j)\ 2 + \M(kJ,i;;j)\ 2 
Re(^(M(j,i,l,k) + M(j,l,i,k) + M(j,l,k,i) + M{k,i,l,j) 

+M(k,j,i,i)+M(i,i,k,j)+M(i,j,i,k)+M(l,k,j,i)) x M(i,j,k,l)^ 
+Ref(M(k,i,l,j) +M(k,l,i,j) +M(k,l,j,i) + M{j,i,l,k) 

+M(j, k, 1, i) + M{1, k) +M{1, k, + M{l,j, k, i)) x M{i, k,j, Z) f 
+Re((M(k,i, + M(k,j,i,l) +M(k,j,l,i) + M(l,i,j,k) 

+M(l,k,j,i)+M(j,i,l,k) + M(j,k,i,l) + M{j,l,k,i)) x M{i,k,l,j) ] 
+Re(^(M(i, k, l,j) + M(i, 1, k,j) + M(i, l,j, k) + M(j, k, 1, i) 

+M(j,i,l,k) + M(l,k,j,i) + M(l,i,k,j) + M(l,j,i,k)) x M(k,i,j,l)^ 
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+N 2 



+Re\(M(i,k,j,l) + M(i,j,k,l) + M(i,j,l,k)+M(l,k,j,i) 

+M(l, i, j, k) + M(j, k, I, i) + M(j, i, k, 1) +M(j, I, i, k)) x M(k, i, I, j) f 
+Re ^ (M (/, k, j, i) + M(l, j,k,i)+M (/, j,i,k)+M{i, k, j, 1) 

+M(iJj,k) + M(j,k,i,i) + M(j,i,k,i)+M(j,i,l,k)) x M(k,l,i,tf 
\M(i,j;;k,l)\ 2 + \M(i,k;;j,l)\ 2 + \M(k,i;; j,l)\ 2 + \\M(k,l;;i,k)\ 



2 

(iV c 4 -3iV 2 - 1), -,. . t f. l2 
+— — ^ '-\M(i,j,k,l)\ 2 

(7.53) 

The expression above is obtained as follows. We start by expanding the sum over permu- 
tations for all four gluons present in the final state in the colour decomposition for the 
unphysical process — > QQgggg , given in eq.(7.49). Then, gluons 5 g and 6 g are crossed to 
the initial state. Finally, the terms are recombined in a double sum over permutations for 
the two initial state and the two final state gluons. The corresponding subtraction term 
reads, 

d& gg^QQgg = ~ 1 ) d$ 4(^iQ, k 2 Q, k 3g , h g ;p 5g ,p 6g ) 

(i,j)€P(3,4),(fc,0eP(5,6) 
r3 



|^3(^;lQ,2 ) (|^°((12) Q ,i 5 ,^,^,Z Q )| 2 + (12)q)| 2 

+|7wO((12)q,^,^,j 5 ,Zq)| 2 + iTWgcf^,^,^,^, (fs:)^)! 2 

+ \M° 5 ((12) Q , kg,i g ,jg, l Q )\ 2 + \M 5 (!q, kg,i g ,jg, ( 12 ) ^ ) | 2 ) jf (fcjg, k 3 , k 4 ) 
+ 4(l Q ,i g ,j g )\M 5 ((Ti) Q , {ij)g, kg, tg, 2 Q ) | 2 jf ( fcg , fc 2 , fc ~- ) 

+rf3(2 <5 ,^,J 9 )|A^°(l Q ,fc 5 ,Z 5 ,^) 3 , (^Oq)! 2 ^ 3 ^^,^.,^.) 
+^3^; iff, 1q) (l>tg((fi)o, t g ,jg, k g , 2 Q )\ 2 

+\M° 5 ((Ti) Q , t g , k g ,j g , 2 Q )) \ 2 4 3 Hk Ti , k 2 , kj) 

+D° 3 (l g ;i g , 2 Q ) (\M° 5 (l Q ,j g , k g , t g , (2*)q)| 2 

+ \M 5 (l Q ,kg,jg,tg,(^) Q )\ 2 ) jf , %• , ^ ) 

+fi(.l 9 -,igJg) (\M° 5 (1q, (ij) g , t g , Jc g , 2 Q )\ 2 + \M° 5 (l Q , l g , (ij) g , k g , 2 Q )\ 2 

+|^g(l Q , ^, 4, 2 )| 2 + |^fg(l Q , ^, 4, 2 )| 2 ) ^J 3) (An, fc 2 , ^) 

+F^(k g ,l g ;i g ) (|AtO(l Q ,l g ,f g ,^,%)| 2 + |Atg(lQ,jvJ g ,f g ,2 a )| 2 ) J^XhMM) 
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A° 3 (l Q ,i g ,2 Q ) (\M%((li) Q J g ,k g ,l g ,(2i) Q )\ 2 - (l/2)\Ml((li) Q J 1 ,k 1 ,l 1 ,(2i) Q )\ 2 
+|A^g((fi)Q,A ;ff ,^,Z 9 ,(2i) ( 5)| 2 + |A^garO Q ,fc p ,4,i^(2i) )| 2 

+MM° 5 ((Ti)Q,jg, kg, lg, (2i) Q )M° 5 ((Ti)Q,j g , lg, kg, (2i) Q ?) 

+MM° 5 ((Ti)Q,kgJ g J g ,(2i) Q )M 5 ^ 

+Re(>t°((k2Qi g j;,i ff ,(2z2Q)^((lOo>4>^>i^( 2 OQ) t ) 
+2Re(M 5 ((li) Q ,jg, kg, lg, {2i) Q )Ml{{li) Q , k g ,j g , l g , {2i) Q ?) 

+ 2Re(M 5 ((li) Q ,jg, kg, lg, (2i ) Q ) M I ( ( U ) Q , lg, kg,jg, (2i) Q )^) 

+ 2Re(M 5 ((li) Q , kg,jg, lg, (2z)q)^((H) Q , kg, lg,jg, (2*)^)) jf (fcg, , , ) 

lQ, 2 ) (|^rg((12) Q , ^, ^, Zq)| 2 + |^°(f , (12)q)| 2 

+|^rg((i2) Q , ^, f Q )| 2 + |^°(f , ^, (f2) )| 2 

MMl{{l2) Q ,k g ,i g ,j g j Q )\ 2 + \Ml{t Q ,k^ 
+|yW°((12) Q ,^,j 9 ,^,Z Q )| 2 + |yWg(f ,z 9 ,j^,^, (12) )| 2 

+i^((r2)Q, i9 ,^,i 7 j Q )i 2 + i^(fQ,i 9 ,^,i 7 ,(r2) ( 3)i 2 
+iA^((ii>)Q,^, V j 7 jQ)i 2 + ^ 
-(i/2)i^((r2)Q,i 7 ,j 7 ,fc 7 ,f Q )i 2 
-(i/2)|xO(f ),i 7 ,i 7 ,fc 7 ,(r2)Q)i 2 ) 4 3 \k^,ki,kj) 

+4(i Q ,i g ,j g ) (i^((ri) Q ,fe%,fe 9 ,/ 7 ,2Q)i 2 

+Re(At°((H) Q , fcg, i fl , 2q)^((H)q, ^ k g , 2g)t) 
-Ite(Mg((li)o, 2 (5 )M°((li) , ^ 2g)t) 

-Ite(Mg((li)o, 4, (u)^ 2q)^°((1z2q, f 9 , Mi?%, 2g)t) 
+2Re(M°((1i)Q, l g , (ij) g , 2 Q )M° 5 ((li) Q , l g , (ij) g , k g , 2 Q )t) 
-2Re(>t°((H) Q , 4, 2q)A^°((H)q, 2Q) 1 ") 

-2Re(M°((H) Q , 2q)M° 5 ((U) q , l g , k g , {ij) g , 2 Q ^j) jf\k^, k 2 , k~) 

+4{2 Q ,i g ,j g ) (|^(1q,^,^%,/ 7 ,(2z) )| 2 

+Re{Ml{l Q ,k g: {ij) g ,l g ,{2i^ 

-Re(M° 5 (l Q , (ij) g , k g , lg, (2i ) q ) M % ( 1 q , (2*)^) 
-Re(M°(l Q , fc g J ff , (22q)^^(1q, f 9 , Mu%, ^q) 1 ) 
+2Re(^°(l Q , (ij) g , k gi l g , (2i) Q )M° 5 (l Q ,i g , (ij) g , k g , (2i) Q ?) 
-2Re(M°(l Q , fc fl) (ij) g , (2i) Q )M° 5 (l Q ,i g , (ij) g , k g , (2i) <5 )t) 

-2Re(>t°(l Q ,to%i 9 J 9 ,^ 
+Z?3°(/ 9 ; i 9 , 1 Q ) (|Mg((fi) , U 9 , fcy, 2 )| 2 + |Mg((fi) 0) l g , k g ,j y , 2q)| 2 



+2Re(M° 5 ((Ti) Q ,j g , T g , k g , 2 q )M° 5 ((U)q, k g ,j g , T g , 2$)t) 
+2Re(M° 5 ((ii) Q ,j g , t g , kg, 2 Q )M° 5 ((Ti) Q , k g , l g ,j g , 2 Q ?) 

+ 2Re(M 5 ((U) Q ,jg, kg, Ig, 2 Q ) ((U) Q , kg, tg,jg, 2 Q f) 

-2Re(M 5 ((U) Q ,jg, tg, kg, 2 Q ) ((li) Q , tg, kg,jg, 2^) 

-2Re(M° 5 ((ii) Q ,jg, kg, tg, 2 Q )Ml((Ti) Q , t g , k g ,j g , 2 <3 )t) 

-2Re(M 5 ((li) Q , tg, jg, kg, 2 Q ) M^U) Q , tg, kg,jg, 2^) 

-2Re(^°((ll) ,i fl ,fc fl ,f ff ,2 a )^((H) g i g ,i g ,f g ,2Q)t) 
-2Re(M° 5 ((li) Q , lg,j g , kg, 2 Q )M° 5 ((li) Q , k g ,j g , l g , 2g)t) 

-2Re(M 5 ({li) Q ,tg,jg, kg, 2 Q ) M^U) Q , kg, lg,j g , 2^)) | 2 (fcg, fc 2 , fcj ) 

+^3(^^^2q) (|A^i(lQ,i ff ,4,^,(^) )| 2 + |A^g(lQ,fc 5 ,f ff ,j T ,(2i) (3 )| 2 
+2Re(M°(l , f fl , fc fl) (2z)q)A4°(1q, f g , j g , (2*)$)*) 
+2Re(7Wg(l Q , j s , f fl , kg, (2i) Q )Ml(l Q , k g , l g ,j g , (2i) Q )^) 
+2Re(AtO(l Q ,f fl ,i fl ,fc fl ,(2l)Q)^0(l Q i g ,f g ,j g ,(2z)Q)t) 
-2Re(^°(l ,i fl i fl ,f fl ,(2l)Q)>lO(lQj g ,^,j g ,(2l)Q)t) 
-2Re(^g(lQ,4,i s ,fc s ,(2^) )^0(lQ,4,^,i 9 ,(2^)Q)t) 
-2Re(Mg(lo,^, f fl , kg, (2i) Q )M° 5 (l Q , k g ,j g , l g , (2*)$)*) 

-2Re(Ml(l Q ,jg,kg,tg,&)Q)Ml(l Q ,kg,jg,tg,&) Q )^ 
-2Re(M 5 (l Q , tg,jg, kg, (2i ) Q )M jj ( 1 Q , kg,jg, lg, (2l) ) 

-2Re(M° 5 (l Q ,j g , kg, tg, (2i) Q )M° 5 (l Q , k g , t g ,j g , (2i) Q ?j) jf (h,k^, k 3 ) 
+fs(lg\ igJg) (\M° 5 (1q, (ij) g , l g , K 2 Q )\ 2 + \M° 5 (l Q ,l g , (ij) g , fc 7 , 2 Q )\ 2 

+ 4Re(Ml(l Q , {ij)g, kg, tg, 2 Q )Ml(l Q ,tg, kg, {ij)g, 2 Q )^) 

-4Re(M° 5 (l Q , {ij) g , tg, kg, 2 Q )M° 5 (1 Q , k g , {ij) g , l g , 2$)t) 

-4Re(M 5 (lQ, tg, {ij)g, kg, 2 Q )M%{l Q , kg, {ij)g, lg, 2 Q )^) 
-4Re(M°(l Q , (ij)g, tg, kg, 2q)^(1 Q , kg, tg, (ij)g, 2^) 

-4Re(.M°(l Q , f fl , kg, 2 Q )M° 5 (1 Q , k g , t g , (ij) g , 2$)t)) jf } fc 2 , fc..) 

+ 2Re(Ml(l Q , kg,jg, tg, 2q)A^(1 Q , tg,jg, kg, 2^) 

-4Re(M° 5 (l Q ,j g , kg, t g , 2 Q )M%{l Q , kg, tg,j g , 2$)t) 
-ARe(M° 5 (l Q ,jg, %, t g , 2 Q )M° 5 (1 Q , t g , k g ,j g , 2$)t)) jf } (k, k, k) 
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1 



A° 3 (l Q ,i g ,2 Q ) (|A^(0((H)q,J 9 ,^,^,(20q)| 2 + |A^g((H)Q,^,i 5 ,^,(2i)Q)| 5 
+\M° 5 ((U) Q ,k g ,i g ,j 1 ,(2i) Q )\ 2 

-(3/2)|>t°((ll) Q ,i 7 ,fc 7 ,f 7 ,(2z) )| 2 ) jf (k^k^kj) 



+^i(^;lQ,2 ) (|A^g((r2) Q ,^,^,^,f Q )| 2 + |A^°(f Q ,^,^,fc T ,(12) )| 2 

MMl((l2) Q ,i g rk g J 7 M 2 + \Ml(t Q ,i g rk g J 7 ,(^)Q)\ 2 
+ |A^0((12) c3 ,^,i, ? ,^,fQ)| 2 + |7Wg(f a ,^,^,jV, (12)q)| 2 ) ^ 3) (A;^, ^, 

+l4(lQ,igJg)\M 5 ((U) Q , ky, l y , 2 Q )\ 2 J? } (fc s , k 2 , k fj ) 

+Z?°(/ 9 ;^,lQ)|A^°((ll)Q,i 7 ,fc 7 ,f 7 ,2Q)|Vf ) (fc ll ,A; 2 ,/c i ) 
+Z?0(/,;^,2Q)|^(lQ,j 7 ,fc 7 J 7 ,(2l)Q)|Vf ) (/c 1 ,^,fe J ) 

i g , 2 Q )\M° 5 ((U) Q ,h, ky, / 7 , (20 )| 2 Jf (fe s , fc,-) 
1 Q , 2q) (|Mg((f2) , i 7 , j 7 , L,, h)\ 2 
+|A^5(fg, iy, j 7 , k 7 , (12)g)| 2 ^ Jg ^ (fcj^, ki, kj 

For this subtraction term, we have also checked that in all its collinear and quasi- 
collinear limits, it reduces to the product of a Casimir factor multiplied by a splitting 
function corresponding to the given limit and the appropriate non-colour ordered matrix- 
element squared as explained in Section 7.4.1, providing us with a powerful check on the 
correctness of our result. 
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8. Conclusions 

We have presented the extension of the antenna formalism required for the calculation of 
hadronic processes involving massive final states in association with jets at the NLO level. 
The construction of massive subtraction terms with all its mass-dependent ingredients is 
presented in all required configurations (final-final, initial-final, initial- initial) . The un- 
known massive antenna functions are derived, their limiting behaviour is presented and 
those are finally integrated over a factorised form of the massive phase space. Besides 
the massive extension of the flavour-conserving antennae, new massive flavour-violating 
antennae were derived in unintegrated and integrated forms. One of the integrated mas- 
sive initial-final antenna, •A^qq-, can be directly related to the well-known heavy quark 
coefficient function. A special section is dedicated to this comparison and full agreement 
is found. In Section 5, when all antennae are integrated over the appropriate factorised 
massive phase space, we showed that we can capture all poles of the massive integrated 
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antennae in universal factors. Those poles are related either to massive colour ordered 1^ ■ - 
type operators or well-known splitting kernels P(ij){x) associated to initial-final massless 
collinear singularities. The colour-ordered massive I^-type operators are here presented 
for the first time. As a first application of our massive extension of the antenna formal- 
ism, we constructed the colour ordered real contributions and subtraction terms for the 
production of a top quark pair and for the production of a top quark pair and a jet at 
NLO. In the second case, the presence of interference terms in the colour decomposition 
of the real matrix element squared renders the construction of the subtraction terms more 
involved. The treatment of those terms is explained in detail in the paper in Section 7.2.1. 
All colour-ordered subtraction terms constructed have been checked in all soft, collinear 
and quasi-collinear limits of the real matrix element squared. Furthermore, for each sub- 
traction term, it has been verified that in all collinear and/or quasi-collinear limits present 
in it, the sum of all terms contributing to a given limit add up to reproduce the product 
of the required splitting function with the corresponding Casimir factor multiplied with 
the non colour-ordered matrix element squared. This check is explicitly derived for the 
subtraction term related to the process gg — > QQg contributing to the process pp — > ti at 
NLO. This verification provides us with an extremely powerful test on the correctness of 
our results for the subtraction terms for ti and ti +1 jet production at NLO. 

The results presented in this paper represent a substantial step towards the calculation 
of the NNLO corrections to top quark pair production within the antenna subtraction 
formalism. The decomposition of the real matrix elements into colour-ordered amplitudes 
squared and the identification of the different leading and subleading colour structures, 
including the treatment of interference terms, described in Section 7.2.1, will allow the 
application of the NNLO antenna subtraction method to compute the double real radiation 
contributions to pp — > ti at NNLO. The NLO antenna subtraction terms for pp — > ti +jet, 
provided in this paper, are already part of the NNLO corrections to ti production at hadron 
colliders. 
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